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Preface 


S ome of the primary objectives of the SeaWiFS Project, as stated in Volume 1 of the SeaWiFS Prelaunch 
Technical Report Series (Hooker, et al. 1992) are to achieve a radiometric accuracy to within 5% absolute 
and 1% relative, water- leaving radiances to within 5% absolute, and chlorophyll a concentration to within 35% 
over the range of 0.05-50 mgm -3 . These goals are challenging, because they require very accurate estimation of 
a large number of variables, both radiometric and biological. The SeaWiFS Project, in collaboration with the 
SIMBIOS Project, continues to address calibration, measurement protocol, and data analysis issues associated 
with reducing the uncertainties in the key variables required to achieve the stated measurement goals. This 
overall effort has now been underway for a decade and continues. The work completed prior to the launch of 
SeaWiFS is documented in Hooker and McClain (2000). All 14 volumes of the SeaWiFS Postlaunch Technical 
Report Series have been dedicated to calibration and validation topics (these include two index volumes). 

Volume 15 is the report of the second Data Analysis Round Robin executed in March 2000 (DARR-00), which 
focused on variations in spectral upwelling radiance, downwelling irradiance, and diffuse attenuation coeffi- 
cients resulting from differences in data processing methods. The first Data Analysis Round Robin (DARR-94) 
demonstrated that different analysis techniques can produce reasonably similar results with so-called good data 
(Siegel et al. 1995); however, with noisy data, the results were inconsistent. DARR-94 did not recommend 
specific techniques as being clearly superior to others. DARR-00 represents a more detailed intercomparison, 
but does not include the same suite of analysis methods. The participants in DARR-00 were the Goddard Space 
Flight Center, the Joint Research Center, and Satlantic, Inc. All three groups have extensive experience and 
have worked on a variety of field programs together. Each has developed its own analysis approach making 
an intercomparison a natural extension of their ongoing collaboration. As might be expected, the results show 
differences in the derived products and in some cases, the differences are substantial. These results are of interest 
to the bio-optics community because many groups have also developed their own analysis methodologies and 
should be aware of the considerations outlined in this study. 

These analysis round robins were conducted on a voluntary basis with the philosophy and purpose of addressing 
difficult processing problems so as to deliver better data products. The DARR-00 group has dealt with a number 
of perplexing issues and has documented their approach and findings in an open and straightforward manner. 
Their efforts and honest appraisals are laudable, and advance the state of the art a step closer to meeting our 
data accuracy requirements. 


Greenbelt, Maryland 
May 2001 


— C. R. McClain 


Results of the Second SeaWiFS Data Analysis Round Robin, March 2000 (DARR-00) 


Table of Contents 

Prologue 1 

1. The Second SeaWiFS Ocean Optics DARR (DARR-00) 4 

1.1 Introduction 4 

1.2 Instrumentation 5 

1.2.1 The AMT Instruments 6 

1.2.2 The AAOT Instruments 7 

1.3 The DARR-00 Database 8 

1.3.1 The GSFC Data 8 

1.3.2 The JRC Data 9 

1.3.3 The Data Products 9 

1.4 Data Analysis Methods 9 

1.5 Results-Standard (VI) 16 

1.5.1 L u (0~) 17 

1.5.2 E d ( O') 19 

1.5.3 K l and K d 20 

1.5.4 E d ( 0 + ) 21 

1.5.5 R ts 21 

1.5.6 [Tw] n 23 

1.5.7 i? u (0 _ ) and Q n ( 0") 25 

1.5.8 Ancillary Variables 25 

1.5.9 VI Summary 27 

1.6 Results-Nonstandard 29 

1.6.1 V2 Summary 29 

1.6.2 V3 Summary 31 

1.6.3 V4 Summary 37 

1.6.4 V5 Summary 37 

1.7 Discussion 40 

1.8 Conclusions 45 

2. The GSFC Data Processing System 46 

2.1 Introduction 46 

2.2 Data Collection 47 

2.3 Data Format 47 

2.4 Processing Steps 47 

2.4.1 Data Conversion 47 

2.4.2 Data Visualization 49 

2.4.3 In- Water Data Products 49 

2.4.4 Above- Water Data Products 50 

2.4.5 Ancillary Data Products 50 

2.5 Summary 50 

3. The JRC Data Processing System 52 

3.1 Introduction 52 

3.2 Data Requirements 53 

3.3 Data Processing Steps 53 

3.3.1 Calibration and Formatting 53 

3.3.2 Filtering and Processing 55 

3.4 Summary 56 


IV 


S. Hooker, G. Zibordi, J-F. Berthon, D. D’Alimonte, S. Maritorena, S. McLean, and J. Sildam 


Table of Contents ( cont .) 

4. The ProSoft Optical Data Processor 57 

4.1 Introduction 57 

4.2 Data Processing 57 

4.2.1 Level -1 Data Processing 58 

4.2.2 Level-2 Data Processing 61 

4.2.3 Level-3 Data Processing 62 

4.3 Summary 65 

Appendix A 66 

Appendix B 66 

Glossary 66 

Symbols 67 

References 69 

The SeaWiFS Postlaunch Technical Report Series 71 


V 



S. Hooker, G. Zibordi, J-F. Berthon, D. D’Alimonte, S. Maritorena, S. McLean, and J. Sildam 


Abstract 

The accurate determination of upper ocean apparent optical properties (AOPs) is essential for the vicarious 
calibration of the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) instrument and the validation of the derived 
data products. To evaluate the importance of data analysis methods upon derived AOP values, the Second Data 
Analysis Round Robin (DARR-00) activity was planned during the latter half of 1999 and executed during March 
2000. The focus of the study was the intercomparison of several standard AOP pareimeters: a) the upwelled 
radiance immediately below the sea surface, L u (0“,A); b) the downward irradiance immediately below the 
sea surface, Ed{ 0“, A); c) the diffuse attenuation coefficients from the upwelling radiance and the downward 
irradiance profiles, Ki{ A) and Kd( A), respectively; d) the incident solar irradiance immediately above the sea 
surface, Ed(0 + , A); e) the remote sensing reflectance, R rs (A); f) the normalized water-leaving radiance, [Lw(A)] n ; 
g) the upward irradiance immediately below the sea surface, £ u (0~), which is used with the upwelled radiance 
to derive the nadir Q-factor immediately below the sea surface, Q n (0~,A); and h) ancillary parameters like 
the solar zenith angle, 0 , and the total chlorophyll a concentration, CV a , derived from the optical data through 
statistical algorithms. In the results reported here, different methodologies from three research groups were 
applied to an identical set of 40 multispectral casts in order to evaluate the degree to which differences in data 
analysis methods influence AOP estimation, and whether any general improvements can be made. The overall 
results of DARR-00 are presented in Chapter 1 and the individual methods used by the three groups and their 
data processors are presented in Chapters 2-4. 


PROLOGUE 

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) 
Project is tasked with executing a program to acquire a 
global ocean color data set, validate and monitor its ac- 
curacy and quality, process the radiometric data into geo- 
physical units using a set of atmospheric and bio-optical 
algorithms, and distribute the final products to the scien- 
tific community through the Goddard Space Flight Cen- 
ter (GSFC) Distributed Active Archive Center (DAAC). 
The long-standing objective of the SeaWiFS Project is 
to produce water-leaving radiances to within 5% absolute 
(Hooker and Esaias 1993). 

The accurate determination of upper ocean apparent 
optical properties (AOPs) is essential for the vicarious cal- 
ibration of an ocean color sensor like SeaWiFS and the val- 
idation of the derived data products, because the sea-truth 
measurements are the reference data to which the satellite 
observations are compared (McClain et al. 1998). The only 
economically feasible approach for minimizing spatial bi- 
ases for a global data set is to maximize the acquisition of 
in situ measurements by soliciting data from the oceano- 
graphic community at large. The SeaWiFS Bio-Optical 
Archive and Storage System (SeaBASS) is the database 
maintained by the Project for these data and SeaWiFS 
sponsored field activities (Hooker et al. 1994). 

The uncertainties associated with in situ AOP mea- 
surements have various sources, such as, the deployment 
and measurement protocols used in the field, the absolute 
calibration of the radiometers, the environmental condi- 
tions encountered during data collection, the conversion 
of the light signals to geophysical units in a data process- 
ing scheme, and the stability of the radiometers in the 


harsh environment they are subjected to during trans- 
port and use. In recent years, progress has been made 
in estimating the magnitude of some of these uncertain- 
ties and in defining procedures for minimizing them. For 
the SeaWiFS Project, the first step in the process of con- 
trolling sources of uncertainty was to convene a workshop 
to draft the SeaWiFS Ocean Optics Protocols (SOOP). 
The SOOP adheres to the Joint Global Ocean Flux Study 
(JGOFS) sampling procedures (Joint Global Ocean Flux 
Study 1991) and defines the standards for optical measure- 
ments to be used in SeaWiFS radiometric validation and 
algorithm development (Mueller and Austin 1992). The 
SOOP is periodically updated as deficiencies are identified 
and outstanding issues are resolved (Mueller and Austin 
1995, and Mueller 2000). 

The follow-on inquiries into controlling sources of un- 
certainty investigated a variety of topics. The SeaWiFS In- 
tercalibration Round-Robin Experiment (SIRREX) activ- 
ity demonstrated that the uncertainties in the traceability 
between the spectral irradiance of calibration lamps were 
approximately 1.0%, and the intercomparisons of sphere 
radiance was approximately 1.5% in absolute spectral ra- 
diance and 0.3% in stability (Mueller et al. 1996). The first 
SeaWiFS Data Analysis Round Robin (DARR-94) showed 
differences in commonly used data processing methods for 
determining L u (0”,A) and £^(0 + ,A) were about 3-4% of 
the aggregate mean estimate (Siegel et al. 1995). Hooker 
and Aiken (1998) made estimates of radiometer stability 
using the SeaWiFS Quality Monitor (SQM), a portable 
and stable light source, and showed the stability of their 
radiance and irradiance sensors in the field during a 36- 
day deployment was on average to within 1.0% (although 
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some channels occasionally performed much worse). More 
recently, Hooker and Maritorena (2000) quantified differ- 
ences in the in-water methods and techniques employed 
for making radiometric measurements in support of Sea- 
WiFS validation and demonstrated a total uncertainty in 
the measurement of in-water AOPs at approximately the 
3% level. 

As part of the SeaWiFS Project calibration and val- 
idation activities, the SeaWiFS field program conducted 
specific experiments to investigate these issues. The ex- 
periments took place during several Atlantic Meridional 
Transect (AMT) cruises on board the Royal Research Ship 
James Clark Ross (JCR) between England and the Falk- 
land Islands. The odd-numbered, southbound cruises sam- 
pled the boreal autumn and austral spring; while the even- 
numbered, northbound cruises sampled the boreal spring 
and austral autumn (Aiken and Hooker 1997). 

The accuracy of any AOP determination is a function 
of the quality of the observational measurement, the data 
acquisition methodology, and the data processing method 
employed. The former includes the quality of the optical 
calibration and the radiometric stability of the instruments 
while they are being used in the field. To minimize observa- 
tional uncertainties, the SeaWiFS Project has sponsored a 
variety of multidisciplinary workshops to outline the obser- 
vations and sampling protocols required for bio-optical al- 
gorithm development (Mueller and Austin 1992 and 1995). 
One of the consequences of the workshops was the estab- 
lishment of the series of SIRREX activities to demonstrate 
and advance the state of the art for calibrating the instru- 
ments used in field activities. 

Although the SeaBASS architecture allows for some 
quality control (Hooker et al. 1994), it is based primar- 
ily on resolving obvious clerical errors in the reporting of 
where and when data acquisition activities took place — it 
does not attempt to quantify differences in the data analy- 
sis methods employed. The latter is, in part, a function of 
how the individual software packages deal with data filter- 
ing, binning, smoothing, and removing possible artifacts 
from changes in surface illumination, ship shadow, or re- 
flections and wave focusing. 

The focus of this intercomparison study is the estima- 
tion of a variety of optical parameters derived from ver- 
tical profiles of the downward irradiance, Ed(z.X), the 
upwelled (nadir) radiance, L u (z, A), and the upward ir- 
radiance, E u (z, A). In the results reported here, different 
methodologies from three international groups are applied 
to the aforementioned multispectral profiles to evaluate the 
degree to which data analysis methods influence AOP es- 
timation, and whether any general improvements can be 
made. In addition to GSFC, there were three groups in- 
volved: 

a) The University of California at Santa Barbara 
(UCSB) Institute for Computational Earth System 
Science (ICESS), 


b) The Joint Research Centre (JRC) Marine Environ- 
ment Unit of the Space Applications Institute, and 

c) Satlantic, Inc. (Halifax, Canada). 

The overall results of the second Data Analysis Round 
Robin (DARR-00) are presented in Chapter 1 and the indi- 
vidual methods of the three groups are presented in Chap- 
ters 2-4, respectively. The participants in the activity are 
given in Appendix A. A summary of the material presented 
in each chapter is given below. 

1. The Second SeaWiFS Ocean Optics DARR 

(DARR-00) 

The focus of this round-robin study was the estimation 
of a variety of commonly used data products derived from 
in-water optical measurements. Eleven parameters impor- 
tant to bio-optical analyses were determined and compared 
from three different in-water optical data processors from 
the JRC, GSFC, and Satlantic, Inc. The parameters were 
calculated for a data set covering a large range of total 
chlorophyll a concentration (0.08-2.43 mg m -3 ). All three 
processors were intercompared using 40 optical profiles; the 
JRC and GSFC processors were further intercompared us- 
ing an additional 10 casts (the larger data set increased the 
amount of data in very clear waters, thereby extending the 
lowest total chlorophyll a concentration to 0.027 mgm' 3 ). 
In terms of overall spectral averages, many of the JRC and 
GSFC results intercompared to within 2.5%, but none of 
the Sa/tlantic results intercompared with the other proces- 
sors at this level. Band-ratio averages, however, frequently 
intercompared to within 2.5% for all processor combina- 
tions, even when the overall spectral averages did not. 
These results suggest a database constructed with pro- 
cessed data from a wide source of contributors will have 
substantially higher uncertainties than a database con- 
structed with raw data which is processed with a single 
processor, although band ratios regularly provide reduced 
uncertainties with respect to individual spectral uncertain- 
ties. 

2. The GSFC Data Processing System 

The GSFC data processing system was designed to cal- 
culate apparent optical properties from underwater profiles 
of upward and downward irradiance, as well as upwelling 
radiance, from a wide variety of open ocean regimes. The 
full processing involves three steps: a) the calibration of 
the data, b) the establishment of the extrapolation inter- 
val, and c) the derivation of the final data products. In the 
first step, special files containing dark voltages are used to 
remove detector offsets and apply a pressure tare to the 
depth sensor. The second step is done interactively while 
viewing all the optical and ancillary (sensor tilts, fluores- 
cence, etc.) fields. The final step is usually executed in 
batch mode using a file of extrapolation limits created in 
the second step. All of the code is written in the Interac- 
tive Data Language (IDL) programming environment from 
Research Systems, Inc. (Boulder, Colorado). 
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3. The JRC Data Processing System 

The JRC Data Processor was developed for the analy- 
sis of field data collected in coastal and marine regions in 
support of ocean color calibration and validation activi- 
ties. The code specifically supports the calibration, for- 
matting, visualization, and processing of in-water optical 
profiles of L u (z, A), Ed{z, A), and E u (z, A) collected with 
the Satlantic seven-channel OCR-200 and OCT200 sen- 
sor series. The user can define a variety of options which 
permit the selection of the extrapolation interval for the 
computation of the near-surface values, and the removal 
of outliers through statistical filtering. Relevant features 
include the possibility to minimize measurement perturba- 
tions such as instrument self shading, bottom effects, and 
superstructure shading (although the latter is restricted to 
measurements performed at the Acqua Alta Oceanographic 
Tower (AAOT) site in the northern Adriatic Sea). Each 
step of the processing is supported by the graphic presen- 
tation of intermediate results. Output products are the 
subsurface values L u (0",A), 2^(0“ , A), and JF n (0”,A); the 
diffuse attenuation coefficients Kd{ A), ), and K u ( A); 

plus the derived products R( 0“, A), R rs ( A), [Lw(A)] n , and 
Qn( O', A). 


4. The ProSoft Optical Data Processor 

The ProSoft data processor was developed for the analy- 
sis of optical data collected with instruments manufac- 
tured by Satlantic, Inc., including profiling and fixed-depth 
(buoy) systems. The data processing steps are organized 
into different levels which follow, as much as possible, the 
procedures given in the Ocean Optics Protocols for Sea- 
WiFS Validation, Revision 1 (Mueller and Austin 1995). 
The processing involves three steps wherein the level -1 raw 
data (usually binary) are converted to calibrated values 
(engineering units) using one of several options for dark 
current correction, the level -2 data are averaged and cor- 
rected for acquisition and performance problems (if any), 
and the level -3 data are used to compute the final data 
products, e.g., R rs , [Lw(A)] n , photosynthetically avail- 
able radiation (PAR), etc. Processing can be clone manu- 
ally (and interactively) by the user, guided by the ProSoft 
processing steps, or in a fully automated batch mode. A 
notable difference between ProSoft and other data proces- 
sors is the diffuse attenuation coefficient can be estimated 
using a traditional estimation of the slope of the optical 
variables within a user-defined extrapolation interval or us- 
ing a new methodology based on multiresolution wavelet 
analysis. 
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Chapter 1 


The Second SeaWiFS Ocean Optics DARR 

(DARR-00) 

Stanford B. Hooker 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 

Giuseppe Zibordi 

Marine Environment Unit , Space Applications Institute 
Joint Research Centre of the European Commission 
Ispra, Italy 

Stephane Maritorena 

UCSB Institute for Computational Earth System Science 
Santa Barbara, California 

Abstract 


The focus of this round-robin study was the estimation of a variety of commonly used data products derived 
from in-water optical measurements. Eleven parameters important to bio-optical analyses were determined 
and compared from three different in-water optical data processors from the JRC, the GSFC, and Satlantic, 
Inc. The parameters were calculated for a data set covering a large range of total chlorophyll a concentration 
(0.08-2.43 mg m” 3 ). All three processors w r ere intercompared using 40 optical profiles; the JRC and GSFC 
processors were further intercompared using an additional 10 casts (the larger data set increased the amount of 
data in very clear waters, thereby extending the lowest total chlorophyll a concentration to 0.027 mgm" 3 ). In 
terms of overall spectral averages, many of the JRC and GSFC results intercompared to within 2.5%, but none 
of the Satlantic results intercompared with the other processors at this level. Band-ratio averages, however, 
frequently intercompared to within 2.5% for all processor combinations, even when the overall spectral averages 
did not. These results suggest a database constructed with processed data from a wide source of contributors 
will have substantially higher uncertainties than a database constructed with raw data which is processed with a 
single processor, although band ratios regularly provide reduced uncertainties with respect to individual spectral 
uncertainties. 


1.1 INTRODUCTION 

The SeaWiFS calibration and validation plan (Hooker 
and McClain 2000) relies on radiometric measurements 
made at sea by a diverse community of investigators. One 
of the long-standing objectives of the SeaWiFS Project is 
to produce spectral water-leaving radiances, Lvv(A), within 
an uncertainty of 5% (Hooker and Esaias 1993), and the 
sea-truth measurements are the reference data to which the 
satellite observations are compared (McClain et al. 1998). 
The accuracy of the field measurements are, therefore, of 
crucial importance. 

If a total 5% uncertainty level is to be maintained for a 
vicarious calibration exercise (remote plus in situ instru- 
mentation), approximately half of the uncertainty budget, 
i.e., 2.5% (actually if quadrature sums are used, the ground 


truth component is closer to 3.5%), is available for the 
ground truth component. Hooker and Maritorena (2000) 
investigated whether or not current deployment and pro- 
cessing schemes meet the 2.5% objective, and found several 
systems almost did. Their study considered four sources 
of uncertainty: calibration, in situ stability, data collec- 
tion, and data processing. They estimated data processing, 
which was restricted to a single data processing package, 
contributed approximately 2% uncertainty to the total un- 
certainty. 

The primary objective of DARR-00 was to use multiple 
data processors to investigate whether or not a data pro- 
cessing uncertainty below the 3-4% result determined for 
DARR-94 is achievable. Two uncertainty thresholds were 
considered for this inquiry: 2.5% as a hoped for minimum 
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(because it represents almost all of the ground truth un- 
certainty budget, but is below the DARR-94 results), and 
1.0% as a needed goal (because it permits some expansion 
in the other components of the total uncertainty budget). 

1.2 INSTRUMENTATION 

As part of the SeaWiFS Project calibration and vali- 
dation activities, the SeaWiFS field program participated 
in AMT cruises on board the JCR and in deployments to 
the AAOT site. The AMT Program exploits the passage 
of the JCR as it transits more than 100° of latitude and 
50° of longitude in the North and South Atlantic Oceans 
between the UK and the Falkland Islands in support of 
British Antarctic expeditions. In September, the JCR sails 
from the UK, and the following April it makes the return 
trip with the vast majority of the sampling taking place in 
Case-1 waters (Aiken et al. 2000). 

The AAOT is located in the northern Adriatic Sea 
(12.51°E,45.31°N) approximately 15 km east of the city of 
Venice (Italy). The site has been visited on a monthly basis 
since the fall of 1995 in support of ocean color calibration 
and validation activities. Because of its near-coastal loca- 
tion, the site may be characterized by Case-1 or Case-2 
conditions. 

The AOP instruments used on AMT and AAOT de- 
ployments included the SeaWiFS Optical Profiling System 
(SeaOPS) and the Low-Cost NASA Environmental Sam- 
pling System (LoCNESS) for the former, and the Wire- 
Stabilized Profiling Environmental Radiometer (WiSPER) 
and a variant of the miniature NASA Environmental Sam- 
pling System (miniNESS) for the latter. SeaOPS and WiS- 
PER are deployed using a winch and crane, whereas, LoC- 
NESS and miniNESS are floated away from the sampling 
platform and deployed by hand. 

Although there are many differences between optical 
instruments deployed with winch systems versus free-fall 
units, the primary differences are related to ship-induced 
perturbations, wave motion, and the time required to per- 
form a cast. Optical instruments deployed with winch sys- 
tems potentially suffer from more disadvantages than ad- 
vantages: 

1. Cranes have a limited reach, so platform perturba- 
tions can be a problem. 

2. For a ship, the vessel is not decoupled from the 
ocean surface, so roll and pitch can cause measure- 
ment problems (particularly for irradiance sensors). 

3. Winches and cranes are sophisticated systems, so 
there is a continuing vulnerability associated with 
breakdowns, especially with the hydraulic subsys- 
tems (this vulnerability is much reduced in the shal- 
low coastal environment where small electric winch- 
es are the most frequently used). 

4. At sea, winches and cranes are stowed during pas- 
sage, so they require relatively lengthy preparation 


time for operations to begin and end, which means 
a quick cast during optimal sky conditions is fre- 
quently difficult to achieve. 

5. The instruments are not submerged until they are 
far from the side of the platform, so there is very 
little chance of them being damaged by wave action. 

6. Many winches have low wire speeds, so the cycle 
rate for a complete cast is relatively long (in shallow 
water this is not the case). A long cycle rate means 
cloud contamination during a cast is likely, and the 
cast frequently must be temporarily halted to allow 
the cloud to pass, which adds to the time required 
(of course, a low wire speed means high vertical 
sampling resolution). 

Tethered (or free-fall) systems, in comparison, have more 
strengths than weaknesses: 

A. The profiler (and sometimes the reference) can be 
deployed away from the platform clear of any plat- 
form shading or reflections to the in situ light field. 

B. The profiler is not subject to wave action, but it 
must be properly trimmed to ensure minimal tilts 
during descent. 

C. A free-floating reference is not decoupled from sur- 
face motion, but engineering solutions and deploy- 
ment practices can be adopted to reduce this effect. 

D. There is a direct cable connection between the in- 
struments and the data acquisition units, so there is 
no complicated (hydraulic) machinery or electrical 
(slip ring) connection which can require long repair 
times in the event of a failure. 

E. The cable is usually a lighter material than steel 
(e.g., Kevlar® f), and is more easily damaged than 
the standard hydrowire used with winch systems. 

F. The instruments are usually hand lowered close to 
the side of the ship, so they are vulnerable to dam- 
age by wave action. 

G. Deployment and recovery can be accomplished with 
only two scientists which ensures a rapid cycle rate. 

The latter is particularly important, because it means casts 
can be executed in between cloud passage and more casts 
can be done in a particular unit of time. It also means 
station scheduling can be kept informal with the ship being 
stopped only when the illumination conditions are optimal. 

All of the radiometers, including any spares, were man- 
ufactured by Satlantic, Inc. (Halifax, Canada). The basic 
subcomponents of the light systems were the power and 
telemetry units, called DATA-100s, the series 200 Ocean 
Color Irradiance (OCT200) sensors, and the series 200 
Ocean Color Radiance (OCR-200) sensors. The radiome- 
ters have seven channels, which were chosen to correspond 


f Kevlar is a registered trademark of E.I. du Pont de Nemours 
and Company, Wilmington, Delaware. 
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with the SeaWiFS instrument wavelengths and bandwidths 
(Hooker et al. 1993). 


1.2.1 The AMT Instruments 


SeaOPS is composed of two sets of instruments. The 
in-water component is composed of a downward-pointing 
radiance sensor (OCR-200), which measures L u (z, A), an 
upward-pointing irradiance sensor (OCT200), which mea- 
sures Ed(z. A), and a downward-pointing irradiance sensor 
(OCT200), which measures E u (z. A). Internal tilt sensors 
quantify the vertical orientation (<£>) of the profiler as it 
is lowered through the water. The above- water compo- 
nent. an OCT200. is fitted to a pole on a mast and mea- 
sures the incident solar irradiance above the ocean surface. 
E c i( 0 + , A). A custom-built T-shaped frame is used to carry 
SeaOPS. which is lowered and raised through the water 
column bv a winch: data are collected during the lowering 
and raising of the frame as down and up casts, respectively. 

The positioning of the equipment on the SeaOPS frame 
was developed with a geometry that ensured the light sen- 
sors were in close proximity to one another while prevent- 
ing the radiance sensor from viewing any part of the sup- 
port (Fig. 1). The narrow geometry of the frame was de- 
signed to provide a minimal optical cross section. The field 
of view of the irradiance sensor was only influenced by the 
7 mm winch wire. 
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Fig. 1 . A detailed layout of the SeaOPS com- 
ponents showing the dimensional relationships be- 
tween the pieces of equipment. 


For most stations, the sun was kept on the same side from 
which SeaOPS was deployed except if adverse sea or wind 
conditions prevented this. The crane used had about a 
10 m reach over the side of the ship. 

The LoCNESS profiler is not a new instrument per se, 
but instead is built up from the SeaOPS components. Sea- 
OPS has two pairs of internal tilt sensors, one pair for 
when it is oriented horizontally and one pair for when it is 
oriented vertically for the LoCNESS configuration. In its 
prototype configuration, LoCNESS measured two compo- 
nents of the light field (Aiken et al. 1998); the DATA-100 
and the two light sensors were connected in line using ex- 
tension brackets with an OCR-200 at the nose (pointing 
downwards to measure L u ) and an OCT200 at the tail 
(pointing upwards to measure Ed)- The addition of weight 
to the nose bracket and buoyant (foam) fins to the tail 
bracket produced a balanced package that fell through the 
water column with minimum tilts (less than 2°). 

The Three-Headed Optical Recorder (THOR) option 
allowed LoCNESS to measure the same three components 
of the in- water light field as SeaOPS (Fig. 2). In all con- 
figuration options, the power and telemetry cable to the 
profiler extended through the field of view of the irradi- 
ance sensor, but the small diameter of the cable (7 mm) 
minimized any negative effects on the measured light field. 
The irradiance incident at the sea surface was measured 
with a sensor installed on a mast situated above the ship’s 
deck (often called the reference or deck cell measurement). 



Careful attention was paid to the balance of the instru- 
ment frame, even though SeaOPS had tilt and roll sensors. 
At the start of each cruise, the frame was trimmed with 
lead weights in air, accounting for the in- water weights of 
the sensors: final trim checks were carried out in situ dur- 
ing the first (test) station. The typical lowering and rais- 
ing speed of the winch was approximately 0.20-0.25 ms -1 . 


Fig. 2. A schematic of the LoCNESS profiler. 

The principle advantage of LoCNESS in comparison to 
conventional free-fall profilers is its cost and flexibility; it 
can be assembled from relatively low cost components and 
it can be quickly reconfigured, because the radiometers 
used are not integral to the design. For example, the two 
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Table 1 . Channel numbers (A*) and center wavelengths (in nanometers) for the radiometers used with the 
AMT sampling systems. 


Sensor 

SeaOPS (SO) 

LoCNESS (LN) 

Type 

Ai 

a 2 

A 3 

A 4 

A 5 

A 6 

A 7 

Ai 

a 2 

a 3 

A 4 

A 5 

^6 

A 7 

L u (z) 

411.0 

442.8 

489.8 

509.7 

555.0 

664.8 

682.7 

411.6 

442.7 

489.9 

510.3 

554.2 

665.3 

682.6 

Ed{z) 

411.5 

442.5 

489.3 

509.6 

555.4 

665.7 

683.2 

411.3 

442.5 

489.3 

509.1 

554.8 

666.0 

682.9 

E u {z) 

411.4 

442.7 

490.0 

509.3 

554.3 

665.9 

682.4 

411.6 

443.2 

490.8 

510.2 

554.9 

665.5 

682.4 

E d ( 0 + ) 

411.9 

443.0 

489.8 

511.0 

555.5 

665.2 

683.7 

411.9 

443.0 

489.8 

511.0 

555.5 

665.2 

683.7 


different nose sensors can be replaced with two identical 
sensors so they can be intercompared in situ. 

A summary of the center wavelengths for the SeaOPS 
and LoCNESS radiometers is given in Table 1. All of the 
instruments have very similar center wavelengths, and a 
common solar reference was used for the two systems. 

1.2.2 The AAOT Instruments 

The WiSPER system is permanently installed on the 
AAOT and is operated from a 7 m platform extension on 
the southeastern side of the second level. A custom-built 
profiling rig and the positioning of the equipment on the rig 
were developed with a geometry that ensures all radiome- 
ters do not view any part of the mechanical supports. The 
radiometers are mounted on a 1 m extension boom which 
displaces them from the bulk of the underwater instrumen- 
tation and puts them approximately 7.5 m from the nearest 
tower leg. Two taught wires anchored between the tower 
and the sea bottom prevent the movement of the rig out 
of the vertical plane defined by the wires. 

The narrow geometry of the WiSPER frame was de- 
signed to provide a minimal optical cross section. The field 
of view of the irradiance sensor is obstructed by the power 
and telemetry cable as well as the stabilization wires, but 
all of these have very small cross sections. Careful atten- 
tion was paid to the rigidity and stability of the rig, so 
there is no need for tilt or roll sensors. 

WiSPER uses the same kind of optical sensors (the 7- 
channel OCR-200 and OCT200 series of instruments) and 
makes the same measurements as SeaOPS: one OCT200 to 
measure Ed{z, A), one OCR-200 to measure L u (z, A), and 
one OCT200 to measure E u (z,\). One difference between 
SeaOPS and WiSPER, is the solar irradiance sensor for the 
latter is regularly occulted at the end of a cast to determine 
the diffuse (or indirect) solar irradiance, ^(0 + ,A). 

A special kind of DATA- 100, the MOR-005, provides 
the analog- to-digital (A/D) conversion and telemetry ca- 
pability for the WiSPER instruments. The equipment 
is powered directly from 12 V lead- acid batteries which 
are stored and kept charged on the tower. WiSPER is 
raised and lowered using a small electrical winch. The 
typical lowering and raising speed of the winch is approx- 
imately 0.08ms" 1 . Visual sea- and sky-state observations 
are recorded during all casts. A schematic of the WiSPER 
system is shown in Fig. 3. 



Fig. 3. The measurements made with the WiS- 
PER system. The in-water sensors are shown in a 
distorted configuration, so each can be seen clearly. 
Their actual physical layout is similar to SeaOPS. 

The miniNESS profiler makes the same measurements 
as WiSPER, except it is deployed as a tethered, free-fall 
package (and has internal tilt sensors). It is a variant of 
LoCNESS and is built with the same modular components: 
a DATAT00 (with 16-bit A/D converters) for power and 
telemetry, and 7-channel OCR-200 and OCT200 sensors. 
The main difference between the LoCNESS and miniNESS 
profilers is the former has all of the light sensors mounted 
at the ends of the rocket-shaped profiler, whereas the latter 
has the E d and L u sensors mounted on the fins and E u on 
the nose. An adapter is used to mount the E u sensor on 
the nose (a flared metal cage can be used to protect this 
sensor against accidental bottom impacts). 

Putting light sensors on the fins destabilizes the profiler 
(although, careful trimming of the profiler has regularly 
produced tilts less than 2°), and it makes the L u sensor 
more susceptible to shading. This problem was minimized 
by choosing where the mechanical termination was with 
respect to the sensors and the sun. In general, the two 
sensor fins, which are 180° apart, will align perpendicular 
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Table 2. Channel numbers (A*) and center wavelengths (in nanometers) for the AAOT sampling systems. 


Sensor 

WiSPER (WP) 

miniNESS (MN) 

Type 

Ai 

A 2 

A 3 

a 4 

A 5 

a 6 

A 7 

Ai 

a 2 

A 3 

a 4 

A 5 

a 6 

A 7 

Lu{z) 

412.3 

442.8 

490.5 

510.8 

554.9 

665.8 

683.9 

412.5 

442.2 

490.9 

510.3 

554.5 

665.4 

684.0 

E d (z) 

412.4 

443.5 

490.6 

509.1 

555.9 

665.4 

682.1 

412.3 

442.1 

490.5 

510.3 

554.5 

665.7 

683.8 

E u (z) 

412.5 

442.2 

490.7 

509.8 

554.7 

664.8 

683.2 

412.4 

443.5 

490.8 

509.9 

554.7 

664.9 

683.2 

411.5 

442.8 

490.9 

510.4 

554.3 

664.8 

682.7 

412.5 

442.2 

490.7 

509.8 

554.7 

664.8 

683.2 

E d ( 0 + ) 

411.5 

442.8 

489.9 

510.3 

554.5 

664.8 

683.2 

411.5 

442.8 

489.9 

510.3 

554.5 

664.8 

683.2 

412.1 

442.9 

490.2 

509.3 

555.2 

665.1 

683.0 

412.1 

442.9 

490.2 

509.3 

555.2 

665.1 

683.0 


to the mechanical termination when the cable is pulled in 
to bring the profiler to the surface (before a profile). To 
minimize L u shading, all that is required is to choose which 
of the other two fins should be used for the mechanical 
termination, so the L u sensor aligns towards the sun. A 
schematic of the miniNESS profiler is shown in Fig. 4. 



Fig. 4. The sensors and corresponding measure- 
ments of the JRC miniNESS profiler. 


A summary of the center wavelengths for the WiSPER 
and miniNESS radiometers is given in Table 2, all of which 
are very similar. Multiple entries for a particular sensor 
type represent alternative configurations. For the solar ref- 
erence measurement, both sampling systems always used 
the same reference. 

1.3 THE DARR-00 DATABASE 

The DARR-00 data processing groups participated as 
volunteers, so to minimize the work associated with ingest- 
ing a potentially wide variety of formats, one data format 
was chosen: tab-delimited files complying with the Amer- 
ican Standard Code for Information Interchange (ASCII), 
i.e., standard text files with the entries separated by tabs. 


Ten optical casts with supporting calibration and ancillary 
documentation from a diverse a set of bio-optical regimes, 
in terms of chlorophyll a concentration, were chosen for 
four different instrument systems (for a total of 40 casts). 
The types of instruments were classified into two broad 
categories based on the deployment method used, free fall 
or winch and crane (20 casts for each), so the importance 
(if any) of deployment methods could be determined. 

To further simplify the comparisons, the instruments 
involved were restricted to the OCR-200 and OCT 200 Sat- 
lantic series of light sensors. Restricting the data to ASCII 
files from one instrument series limited the data sources to 
the JRC and GSFC groups. The GSFC group submitted 
10 casts each from SeaOPS and LoCNESS. while the JRC 
group submitted 10 casts each from WiSPER and mini- 
NESS. The instruments were also classified according to 
the primary sampling environment involved, deep ocean or 
shallow coastal (20 casts for each), which correspond in a 
general sense to Case-1 or Case-2 conditions, respect ively. 

Although the submitted 40 casts covered a wide range 
of conditions, they did not adequately represent very clear 
waters, so an additional 10 casts from very low chloro- 
phyll a conditions were added to the data set (Sect. 1.5). 
Because the participating groups were volunteering their 
time and effort, not all of them could accommodate the 
extra burden of 10 more casts, so the only groups that 
analyzed these extra data were the JRC and GSFC. As 
will be shown subsequently, this had a minimal effect on 
the conclusions derived from the different results. 

1.3.1 The GSFC Data 

The GSFC data are filed using a simple naming scheme: 
MCCHDNNN.EEE, 

where the encoding letters are defined as 

M A one-letter code for each major campaign (A for 
AMT); 

CC A two-digit sequential campaign number: 

H A single letter indicating the hosting instrument for 
the measurement which is based on the serial num- 
ber (S/N) of the DATA- 100 involved (0 for SeaOPS 
S/N 4, K for LoCNESS S/N 11, and L for LoCNESS 
S/N 4); 
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D A single letter indicating the deployment direction 
(D for down and U for up); 

NNN A three-digit sequential cast number; and 

EEE The file extension (SHO for 6 Hz in-water optical 
data, and SHM for 6 Hz solar reference data). 

An inventory of all the data used for the DARR-00 activity 
is presented in Appendix B. 

1.3.2 The JRC Data 

The JRC data are also organized using a simple algo- 
rithm to set the file names: 

ICCSNL.EEE, 

the encoding for which is defined as 

I A one-letter code for each instrument (W for WiS- 
PER and P for the miniNESS profiler; 

CC A two-digit sequential campaign number; 

S A single letter indicating the data type (S for station 
and D for dark); 

N A one-digit station number; 

L A sequential letter incremented for each cast (A for 
1, B for 2, etc.); and 

EEE The file extension (MOR for WiSPER, OCP for mini- 
NESS, and MVD for the solar reference). 

1.3.3 The Data Products 

Water-leaving radiance is a primary variable in the cal- 
ibration and validation process, not only because it is a 
first-order quantity (i.e., it appears explicitly in the radia- 
tive transfer equation governing the remote sensing obser- 
vations), but also because it is the central physical quantity 
for bio-optical studies in the upper ocean. Many variables 
used for the latter (e.g., R rs (A)) are derived directly from 
Lw{ A) and are the input data for the algorithms used to 
invert ocean color observations into chlorophyll a concen- 
trations. 

Many in- water techniques for calculating Lw{X) are 
based on formulations well expressed by the Smith and 
Baker (1984) method. Variations are a consequence of 
the procedures (and platforms) used to acquire the data, 
and how the in-water data are propagated to the surface. 
The L u (\,z) near-surface profile is usually used to com- 
pute the diffuse attenuation coefficient, Kl(\,Zq), as the 
local (around the depth z 0 ) slope of ln[L u (A, z)]. Ki(X , zq) 
is used to extrapolate the upwelled radiance through the 
upper layer to null depth. Lw{ A), is then obtained by 
properly propagating the near-surface upwelling radiance, 
L u (A,0~), through the air-sea interface. 

The data products compared are those commonly as- 
sociated with bio-optical analyses: 

1. The upwelled radiance immediately below the sea 
surface, L u ( 0",A); 


2. The downward ir radiance immediately below the 
sea surface, E d ( 0", A); 

3. The diffuse attenuation coefficients from the up- 
welling radiance and the downward irradiance pro- 
files, Ki(X) and K d {\ ), respectively; 

4. The incident solar irradiance immediately above the 
sea surface, E d { 0 + , A); 

5. The remote sensing reflectance, R rs ( A); 

6. The normalized water-leaving radiance, [Lw(A)] n ; 

7. The upward irradiance immediately below the sea 
surface, E u (0~), which is used with the upwelled 
radiance to derive the nadir Q-factor immediately 
below the sea surface, Q n (0', A); and 

8. Ancillary parameters like the solar zenith angle, 0, 
and the total chlorophyll a concentration, CV a , de- 
rived from the optical data. 

The latter is usually determined using one of the standard 
ocean color (OC) algorithms (O’Reilly et al. 1998). For the 
DARR-00 activity, the OC2v2 algorithm is intercompared, 
because all of the processors used this algorithm. 

Plots of calibrated optical measurements for selected 
casts from the JRC data set are shown in Figs. 5-7. A 
subset of wavelengths are shown to maintain clarity of pre- 
sentation. The first file, W56S7A from WiSPER (Fig. 5), 
was selected to show the influence of multiple stratification 
on the in-water light fields. The second file, P60S4R from 
miniNESS (Fig. 6), is an example of the significance of sur- 
face effects, i.e., wave focusing, which is most clearly seen 
in the E d (z, A) data. The third file, W56S1A from WiSPER 
(Fig. 7), shows an ideal case with no in- water extra sources 
of variance and very stable solar illumination. 

Plots of calibrated optical measurements for selected 
casts from the GSFC data set are shown in Figs. 8-10. 
Again, a subset of wavelengths are shown to maintain clar- 
ity of presentation. The first file, A07KD025 (Fig. 8), corre- 
sponds to high chlorophyll a conditions (Cr a — 1.18 mgm" 3 ) 
wherein several in-water parameters that characterize the 
water column are confined to the upper layer: the depth 
of the chlorophyll maximum (DCM) as estimated from the 
fluorometer data is 37 m, the mixed layer depth (MLD) is 
29 m, and the one-percent light level (OLL), based on the 
decay of the surface PAR value, is 38 m. The second file, 
A07KD040 (Fig. 9), is an example of a cast in very clear 
water (Cr a =0.08 mgm -3 ) with typical deep ocean water 
properties: the DCM is 139 m, the MLD is 121 m, and the 
OLL is 124m. The last file, A060D033 (Fig. 10), is pre- 
sented to show a case with abnormally large vessel motion 
in eutrophic waters (C Ta =2.43 mg m -3 ) during otherwise 
very good conditions. 

1.4 DATA ANALYSIS METHODS 

In this study, no one data processing system is assumed 
to be more correct than another, so an unbiased parame- 
ter is needed to compare the various methods. The aver- 
age unbiased percent difference (UPD) between two data 
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10-6 io-5 io-4 io-3 io-2 io-l 10° 0.9 1.1 1.3 1.5 


E u (z, A) [W rrr 2 nrir 1 ] Ed{ 0 + ,A) [W rrr 2 nrrr 1 ] 

Fig. 5. The in-water and solar reference data for W56S7A.M0R and W56S7A.MVD, respectively: a) L u (z. A), b) 
Ed(z, A), c) E u (z 1 A), and d) Ed{ 0 + . A). shown as a function of the corresponding z value rather than t. 
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Fig. 6. The in-water and solar reference data collected for P60S4R.OCP and P60S4R.MVD, respectively: a) 
L u (z, A), b) E d (z, A), c) E u {z, A), and d) E d { 0\\). 
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Fig. 9. The in-water and solar reference data collected for A07KD040.SHO and A07KD040 . SHM, respectively 
a) L u (z, A), b) E d (z, A), c) E u (z, A), and d) E d (0 + ,X). 
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Fig. 10. The in-water and solar reference data collected for A060D033.SH0 and A060D033.SHM, respectively: 
a) L u (z, A), b) E d (z, A), c) E u (z,\) t and d) E d ( 0 + ,A). 
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products and was computed as: 


^ S (A) 


200 |g /l (A,ti) - g B ( A,tj)| 
N mfr{ S' 4 (-Mi) + 


where t{ is the time associated with the sample (the i index 
selects a sample), N m is the number of measurements be- 
ing compared, and 5 represents the data product, L u (0~), 
Ed{ 0"), Kd, etc. The A and B codes identify the proces- 
sors involved: J for JRC, G for GSFC, and S for Satlantic. 
The UPD between the upwelled radiances (at null depth) 
estimated with the JRC and GSFC processors, for exam- 
ple, is ?/^ g (A) and represents the uncertainty between the 
two processors in estimating L n (0“,A). 

Spectral averages of the UPD are calculated by includ- 
ing a summation in (1) over the requisite number of spec- 
tral bands, N 


The processing of the data does not always end with the 
calculation of L w { A). The WiSPER and miniNESS data, 
for example, are collected in the proximity of a large struc- 
ture sited in shallow water, so there are other processing 
aspects which need to be considered: 

1. Instrument self-shading, which affects all in- water 
instruments, particularly in the red domain where 
water absorption is elevated; 

2. Bottom reflection, which can be important in shal- 
low water; and 

3. Platform shading, which can effect any instrument 
deployed in close proximity to a large structure (e.g., 
a research vessel or offshore tower). 

Correction terms for these effects are generally not applied, 
because they require accurate knowledge of several atmos- 
pheric and marine optical parameters, and they may re- 
quire extensive computations. 


' t s B = p) 

Although several schemes can be devised for spectral aver- 
ages, the ones considered here are as follows: the average 
of the blue-green wavelengths, Tj^(Ai- 5 ), the average of 
the red wavelengths, (Xq-j), and the average of all the 
wavelengths. ^ s (Ai- 7 ). 

The primary interest in this study are spectral esti- 
mates of water-leaving radiance, but the most frequently 
used ocean color algorithms, OC2v2 (O’Reilly et al. 1998), 
uses band ratios, e.g., the 490 and 555 nm bands. Conse- 
quently, another important parameter for the above-water 
analyses presented here is the UPD for a band ratio: 


(Afc/i) 


200 y, l^V^U) - <B B (\ k n,tj)\ 


where 93 is the band ratio of a particular data product for 
wavelengths k and l (the data product for wavelength k 
is divided by the data product for wavelength /). For this 
study, the only band ratio presented is between 490 and 
555 nm. For example, the average band ratio results for 
L u ( 490) and L u (555) for the JRC and GSFC processors is 
denoted as (A 3 /s)- 

Most results are presented in tables as a function of 
data products and method codes. In these cases, a short- 
hand notation is used for spectral averages and band ratios: 

Ti- 7 The average of all (seven) wavelengths, 412-683 nm; 

Ti _5 The average of the (five) wavelengths in the blue- 
green part of the spectrum, 412-555 nm; 

T 6 -7 The average of the (two) wavelengths in the red part 
of the spectrum, 665-683 nm; and 

4>3/5 The average of the band-ratio results calculated for 
490 and 555 nm (bands three and five, respectively). 


1.5 RESULTS-STANDARD (VI) 

The first results presented are the comparisons derived 
from the processing of data using standard options for each 
processor, i.e., those procedures normally employed by the 
processing groups. These results are also referred to as 

the version 1 (VI) results, and subsequent changes to the 
processors are identified by increasing version numbers. 

In addition to the instrument types (SO, LN, WP, and 
MN), the data sets can be naturally divided into a vari- 
ety of categories depending on the sampling location, the 
deployment methodology, and the total chlorophyll a con- 
centration. The former two are encoded as follows: 

DO For the deep ocean data (SO and LN); 

SC For the shallow coastal data (WP and MN); 

WC For the winch and crane data (SO and WP); and 
FF For the free-fall data (LN and MN). 

The deep ocean data are all from Case-1 conditions, and 
the shallow coastal data are from Case-2 conditions or from 
water near the threshold between Case-1 and Case-2 in 
terms of the Loisel and Morel (1998) classification scheme. 

The Cr a (in milligrams per cubic meter) is used to sepa- 
rate the data into three concentration (or trophic) regimes, 
which are encoded as follows: 

ET For the eutrophic regime ( Cr a >1); 

MT For the mesotrophic regime (0.1 < Cr a < 1); and 
OT For the oligotrophic regime ( Cr a < 0.1). 

These groups are a mixture of all instrument types, except 
all the oligotrophic data come from the LN data set. 

The various combinations of processor intercomparisons 
are identified by two-letter codes: 

JG For the JRC and GSFC processors: 

JS For the JRC and Satlantic processors; and 
GS For the GSFC and Satlantic processors. 
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As will be shown below, the JRC and GSFC processors 
are architecturally similar. Unlike the Satlantic processor, 
they are also maintained by individual investigators with 
narrowly defined user groups, so they can be quickly mod- 
ified without negatively effecting a large community. 

The contribution of each intrument type to the concen- 
tration regimes is presented in Table 3, which is divided 
into two groups given by the subtotal and total partitions. 
The former represents the data used for intercomparing all 
the processors, whereas the latter represents the additional 
data beyond the original 40 casts to ensure a reasonable 
number of samples in each concentration regime (which 
were only processed by the JRC and GSFC groups). 

Table 3. The number of casts from each instru- 
ment type within each concentration regime. Also 
shown are the sampling location and deployment 
method codes for each instrument, as well as the 
groups that processed the indicated casts. 


Instrument 

Groups 

ET 

MT 

OT 

SO (DO and WC) 

J G S 

8 

2 

0 

LN (DO and FF) 

J G S 

1 

6 

3 

WP (SC and WC) 

J G S 

3 

7 

0 

MN (SC and FF) 

J G S 

3 

7 

0 

Original (40-Cast) Subtotal 

15 

22 

3 

SO (DO and WC) 

J G 

0 

0 

3 

LN (DO and FF) 

J G 

0 

0 

7 

Extended (50-Cast) Total 

15 

22 

13 


The distribution of the instrument types within the var- 
ious categories shown in Table 3 will produce biases in 
some of the results. For example, the majority of the shal- 
low coastal data are in the mesotrophic regime, which also 
means the majority of these data are for stable references. 
These biases are not necessarily detrimental to the compar- 
ison process, because they allow a quantitative assessment 
of the differences between the instrument types, sampling 
locations, and deployment methods. It is important to re- 
member these biases, however, because they do mitigate 
the application of the results to broader generalizations. 

The various processing options that can influence the 
final determination of Lvv(A) can have a far-reaching effect 
on bio-optical analyses. Five different processing options 
are considered in the presentation of the results: 

1. Tilt filtering (this is usually a simple threshold — if 
the tilt of the instrument exceeds a particular value, 
the data are not used); 

2. E d (0 + ) normalization (to account for the influence 
of changing illumination conditions during the pro- 
file, some processors normalize the in-water light 
data by E d (0 + , t)/E d (0 + , t 0 ), where t 0 is the start 
time of the profile); 

3. The type of linear fit used in the extrapolation pro- 
cess (simple least-squares regressions are common, 


but some processors use routines wherein any data 
with a large variance are weighted the least); 

4. Whether or not data are statistically filtered, i.e., 
any point more than a specified number of standard 
deviations (na) from the linear regression is consid- 
ered a so-called outlier and is ignored (this is a more 
extreme form of weighted fitting); and 

5. The limits of the extrapolation interval, i.e., the up- 
per and lower depths. 

For the three processors, the standard options used for 
each are as follows: 

J Tilt filtering (5° threshold), E d ( 0 + ) normalization, 
3cr filter for outliers, least-squares (unweighted) lin- 
ear fit, and extrapolation depths chosen to ensure 
all wavelengths are well parameterized; 

G No tilt filtering, no E d ( 0 + ) normalization, not fil- 
tered for outliers, least-squares (unweighted) linear 
fit, and extrapolation depths chosen so the blue- 
green wavelengths are optimally parameterized; and 

S Tilt filtering (5° threshold), no E d { 0 + ) normaliza- 
tion, filtered for outliers, least-squares (unweighted) 
linear fit, and near-surface extrapolation based on 
binned data. 

All of the results presented in this section are based on 
these standard processing options. 

1 . 5.1 L u ( 0 ") 

All the DARR-00 processors calculated L\y( A) from 
profiles of upwelled spectral radiance, L u (z, A), so for this 
study, the most important parameter is arguably L u ( 0“, A), 
i.e., the null depth estimation of the upwelling radiance. 
Unfortunately, with the instruments and deployment meth- 
odologies currently available, it is not practical to measure 
L u ( 0", A) precisely at an infinitesimal depth below the sur- 
face. A profile of L u (z, A) must, therefore, be measured 
over the upper few optical depths with sufficient accu- 
racy to determine Kl(z , A) for L u {z, A), and to propagate 
L u (z , A) to the surface over a near-surface depth interval. 

The percent differences (or uncertainties) between the 
three data processors in the calculation of L u ( 0“, A) is pre- 
sented in Table 4f. There are several aspects of the inter- 
comparisons which are worth noting: 

a) The smallest uncertainties are associated with the 
JRC and GSFC intercomparisons (the Satlantic in- 
tercomparisons are, on average, a little more than 
twice as large); 


| In all of the results tables, the intercomparisons of the data 
processors based on the four instrument types (SO, LN, WP, 
and MN), are for the original 40-cast data set (usually the 
top three boxes of information in the table, but sometimes 
only the top box). The JG comparisons of different sampling 
locations (DO and SC), deployment methods (WC and FF), 
and concentration regimes (OT, MT, and ET), are for the 
extended 50-cast data set, and are always the bottom two 
boxes of information. 
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Table 4. A summary of the ip. 4L and 4> values for the computation of L u ( 0 , A) 


Comparison 

-0(412) 

ip( 443) 

-0(490) 

-0(510) 

-0(555) 

^0(665) 

V>(683) 

*1-5 

$6-7 

*1-7 

*3/5 

so 

JG 

(40) 

3.7 

4.0 

3.2 

3.0 

2.6 

26.4 

11.7 

3.3 

19.0 

7.8 

1.0 

LN 

JG 

(40) 

3.6 

3.0 

2.6 

3.0 

4.0 

23.9 

22.9 

3.2 

23.4 

9.0 

2.0 

WP 

JG 

(40) 

1.0 

1.0 

0.9 

0.9 

0.9 

2.4 

3.2 

0.9 

2.8 

1.5 

0.2 

MN 

JG 

(40) 

2.6 

2.4 

1.9 

1.9 

2.0 

5.5 

7.3 

2.2 

6.4 

3.4 

0.4 

Averag 

e 

2.7 

2.6 

2.1 

2.2 

2.4 

14.6 

11.3 

2.4 

12.9 

5.4 

0.9 

so 

JS 

(40) 

12,8 

24.0 

9.7 

9.3 

8.5 

48.7 

31.3 

12.9 

41.3 

20.6 

2.0 

LN 

JS 

(40) 

7.0 

6.4 

5.4 

5.6 

6.9 

43.1 

42.2 

6.3 

41.0 

16.1 

6.6 

WP 

JS 

(40) 

1.6 

1.3 

1.0 

1.1 

1.2 

3.7 

4.9 

1.3 

4.3 

2.1 

0.2 

MN 

JS 

(40) 

6.9 

5.7 

4.0 

3.9 

3.8 

10.8 

10.2 

4.9 

10.5 

6.5 

0.3 

Averag 

*e 

7.1 

9.3 

5.1 

5.0 

5.1 

.26.6 

22.1 

6.3 

24.3 

11.3 

2.3 

so 

GS 

(40) 

9.7 

21.4 

7.8 

7.4 

7.2 

31.0 

26.8 

10.7 

28.9 

15.9 

1.6 

LN 

GS 

(40) 

8.3 

7.5 

6.6 

6.9 

8.9 

46.8 

45.0 

7.6 

44.3 

18.0 

7.5 

WP 

GS 

(40) 

2.0 

1.8 

1.5 

1.5 

1.6 

4.0 

5.0 

1.7 

4.5 

2.5 

0.3 

MN 

GS 

(40) 

7.5 

6.0 

4.0 

3.8 

3.7 

10.8 

10.1 

5.0 

10.4 

6.6 

0.7 

Average 

6.9 

9.2 

5.0 

4.9 

5.3 

23.1 

21.7 

6.3 

22.0 

10.7 

2.5 

DO 

JG 

(50) 

2.9 

2.8 

2.5 

2.6 

3.0 

25.9 

21.9 

2.7 

23.9 

8.8 

1.5 

sc 

JG 

(50) 

1.8 

1.7 

1.4 

1.4 

1.4 

4.0 

5.2 

1.5 

4.6 

2.4 

0.3 

wc 

JG 

(50) 

2.1 

2.3 

1.8 

1.8 

1.7 

18.5 

12.7 

1.9 

15.6 

5.8 

0.6 

FF 

JG 

(50) 

2.7 

2.5 

2.2 

2.3 

3.0 

16.0 

17.4 

2.5 

16.7 

6.6 

1.4 

OT 

JG 

(50) 

1.6 

1.5 

1.6 

1.8 

2.7 

25.5 

28.7 

1.8 

27.1 

9.1 

1.8 

MT 

JG 

(50) 

2.1 

2.0 

1.8 

1.9 

2.4 

13.4 

12.7 

2.0 

13.0 

5.2 

0.8 

ET 

JG 

(50) 

3.7 

3.7 

2.7 

2.6 

2.2 

15.5 

7.4 

3.0 

11.4 

5.4 

0.7 


Table 5. A summary of the ip. 4/, and 4> values for the computation of Ed ( 0 , A). 


Comparison 

-0(4 12) 

CO 

-0(490) 

0(510) 

0(555) 

0(665) 

0(683) 

*1-5 

$6-7 

*1-7 

$3/5 

SO 

JG (40) 

4.4 

4.8 

4.7 

4.9 

5.5 

8.6 

18.0 

4.9 

13.3 

7.3 

1.6 

LN 

JG (40) 

6.5 

7.9 

5.7 

6.4 

8.3 

20.3 

24.3 

7.0 

22.3 

11.4 

4.3 

WP 

JG (40) 

1.0 

0.8 

0.8 

1.0 

0.9 

1.5 

1.6 

0.9 

1.5 

1.1 

0.7 

MN 

JG (10) 

1.7 

1.5 

2.2 

2.3 

1.8 

2.7 

1.9 

1.9 

2.3 

2.0 

1.1 

Average 

3.4 

3.8 

3.4 

3.7 

4.1 

8.3 

11.5 

3.7 

9.9 

5.4 

1.9 

SO 

JS (40) 

10.3 

12.1 

12.5 

12.2 

12.7 

41.8 

43.7 

12.0 

47.6 

19.5 

2.3 

LN 

JS (40) 

22.2 

20.2 

21.0 

20.9 

22.2 

73.9 

64.8 

21.3 

67.1 

34.1 

7.3 

WP 

JS (40) 

4.0 

3.9 

3.8 

3.5 

5.6 

5.4 

5.2 

4.2 

5.3 

4.5 

3.5 

MN 

JS (40) 

4.7 

4.4 

3.3 

3.8 

3.2 

9.9 

10.2 

3.9 

10.1 

5.7 

1.5 

Average 

10.3 

10.2 

10.1 

10.1 

10.9 

32.8 

31.0 

10.3 

32.5 

15.9 

3.7 

SO 

GS (40) 

7.5 

9.1 

9.8 

9.7 

10.0 

35.5 

34.3 

9.2 

39.9 

15.6 

1.9 

LN 

GS (40) 

26.4 

25.0 

23.9 

23.9 

28.6 

80.7 

76.8 

25.5 

77.5 

40.1 

10.3 

WP 

GS (40) 

3.8 

4.3 

4.0 

3.7 

5.9 

5.8 

5.9 

4.3 

5.8 

4.8 

3.2 

MN 

GS (40) 

5.4 

5.1 

4.2 

4.4 

4.7 

11.2 

10.8 

4.8 

11.0 

6.5 

2.3 

Average 

10.8 

10.9 

10.5 

10.4 

12.3 

33.3 

32.0 

11.0 

33.5 

16.8 

4.4 

DO 

JG (50) 

6.6 

7.5 

7.0 

7.1 

8.1 

15.5 

22.9 

7.3 

19.2 

10.7 

2.5 

SC 

JG (50) 

1.3 

1.2 

1.5 

1.6 

1.3 

2.1 

1.7 

1.4 

1.9 

1.5 

0.9 

WC 

JG (50) 

2.8 

3.0 

3.1 

3.2 

3.3 

5.7 

14.2 

3.1 

10.0 

5.1 

1.1 

FF 

JG (50) 

5.9 

6.7 

6.2 

6.4 

7.2 

14.0 

14.6 

6.5 

14.3 

8.7 

2.5 

OT 

JG (50) 

8.5 

9.6 

9.7 

9.0 

9.7 

16.5 

24.2 

9.3 

20.4 

12.5 

1.5 

MT 

JG (50) 

3.2 

2.8 

3.4 

3.9 

4.1 

8.5 

8.7 

3.5 

8.6 

4.9 

2.0 

ET 

JG (50) 

2.8 

4.2 

2.5 

2.9 

3.5 

7.1 

14.4 

3.2 

10.8 

5.4 

1.9 
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b) All of the intercomparisons have maximum uncer- 
tainties in the red part of the spectrum (665 nm), 
with minimal uncertainties in the blue-green do- 
main (490-510 nm) which increase towards the blue 
(412-443 nm); 

c) The band-ratio results are the smallest uncertainties 
(except for some of the Satlantic intercomparisons); 

d) The shallow coastal data (particulary the WiSPER 
data) have smaller uncertainties than the deep ocean 
data (most notably in the red wavelengths); and 

e) As a function of increasing chlorophyll a concentra- 
tion (OT to MT to ET), there is a general increase 
in uncertainties in the blue-green part of the spec- 
trum, but a decrease in the red wavelengths, and a 
decrease in the band-ratio uncertainties. 

In terms of the desired 2.5% intercomparison objec- 
tive, there are several cases where this threshold is sat- 
isfied. For the Satlantic processor, only the blue-green 
WP uncertainties and some of the band-ratio uncertain- 
ties are below 2.5% — most of the results are, on average, 
significantly above the threshold. For the JRC and GSFC 
processors, all of the WP, and many of the MN uncertain- 
ties, are below 2.5%, as are all of the band-ratio uncertain- 
ties. In addition, many of the shallow coastal plus winch 
and crane average results, as well as the oligotrophic and 
mesotrophic results, are to within 2.5% in the blue-green 
wavelengths, but all the red wavelengths results are above 
the 2.5% threshold. 

1.5.2 E d { 0-) 

Profiles of the downward spectral irradiance, Ed(z, A), 
are required to compute the diffuse attenuation coefficient, 
Kd{z, A), and the irradiance reflectance, R( A). As with 
L u ( 0", A), Ed( 0", A) is usually determined by extrapolation 
over a near-surface depth interval from the profile data, in 
this case, Ej(z, A). The percent differences between the 
three data processors in the calculation of Ed( 0”, A) is pre- 
sented in Table 5. The results partition immediately into 
two groups, irrespective of the processors being compared: 
a) the best results are associated with a stable reference 
(WP, MN, and SC), and b) the worst results are associated 
with an unstabilized reference (SO, LN, and DO). 

The uncertainties associated with the Satlantic proces- 
sor are significantly higher (on average by a factor of 3) 
than the intercomparisons between the JRC and GSFC 
processors, the differences in the red wavelengths are higher 
than the blue and green by more than a factor of 2 (on aver- 
age), and the band-ratio uncertainties are always smaller 
than the individual wavelengths. In terms of the 2.5% 
threshold, the JG uncertainties for the shallow coastal data 
(WP, MN, and SC) are almost always to within 2.5%. 

Regardless of the processors being intercompared, the 
SO and LN (DO) data are always associated with higher 
uncertainties, with the latter being the highest. When 


coupled with the aforementioned shallow coastal data re- 
sults, this produces very good results for the WC and 
SC data categories (the latter are always to within 2.5%), 
and poorer results for the DO and FF categories (which 
are only to within 2.5% for the band-ratio results). The 
uneven distribution of instrument types and deployment 
methods results in a distribution as a function of concen- 
tration regimes which is hard to interpret. The ET data 
have the smallest differences while the OT data have the 
largest, but this is because the former are dominated by 
the SC systems and the latter by the DO systems. 

The overall poorer results obtained with the deep ocean 
data (SO and LN) is a consequence of two factors. First, 
most of these data are in clear water, so wave focusing 
effects can substantially influence where the top of the ex- 
trapolation interval is chosen (the greater vertical resolu- 
tion of the SO data helps reduce this effect). Second, the 
data are from areas with deep mixed layers, so the bottom 
of the extrapolation interval is less constrained than that 
encountered in the shallow coastal environment. Taken 
together, these two factors can result in substantially dif- 
ferent top and bottom depth intervals from different oper- 
ators. If the number of points differ substantially and the 
total possible number of points in the extrapolation inter- 
val is not very large (possible with the high descent rates 
of free-fall instruments), this will effect the extrapolation 
process and, thus, the final Ed{ 0”, A) estimate. 

This point is well quantified by considering the average 
number of points in the extrapolation ranges for the JRC 
and GSFC processing ( Nj and N G , respectively) of the 
SO, LN, WP, and MN data: 

SO Nj = 214 and N G = 345; 

LN Nj = 39 and N G = 68; 

WP Nj — 208 and N G — 216; and 

MN Nj = 20 and N G = 24. 

The percent differences for these data (using the largest 
value as the reference number) are as follows: 38% (SO), 
43% (LN), 4% (WP), and 14% (MN). The ranking of the 
JG results in terms of smallest- to-largest average percent 
differences is WP, MN, SO, and LN — so, the quality of the 
intercomparison is well correlated with how similarly the 
extrapolation intervals are defined. 

Although it is easy to understand why the type of ref- 
erence being used would influence the intercomparison of 
Ed( 0 + ,A) values, it is not so obvious as to why it influ- 
ences the Ej{ 0", A) intercomparisons. The reason is a con- 
sequence of how the two variables are normally used dur- 
ing data processing. For most processors, the Ej( 0‘,A) 
estimate from the in-water profile is properly adjusted for 
the influence of the air-sea interface and then compared 
with the above-water measurement of Ej( 0 + ,A). The two 
should agree to within a few percent, and in most cir- 
cumstances, alternative (but valid) processing options are 
tried until they do. Consequently, there is usually a steer- 
ing effect of the solar reference data on the estimation 
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Table 6. A summary of the 0, $, and $ values for the computation of Kd{ A). 


Comparison 

■0(412) 

■0(443) 

0(490) 

0(510) 

0(555) 

0(665) 

0(683) 

*1-5 

$ 6-7 

$1-7 

$3/5 

SO JG (40) 

3.6 

4.9 

7.5 

7.8 

7.7 

4.1 

7.3 

6.3 

5.7 

6.1 

3.1 

LN JG (40) 

27.7 

32.9 

27.2 

18.0 

25.6 

8.5 

9.7 

34.7 

9.1 

22.8 

21.2 

WP JG (40) 

2.4 

2.4 

3.7 

4.4 

4.8 

1.8 

1.9 

3.5 

1.8 

3.1 

3.4 

MN JG (40) 

4.0 

4.5 

11.6 

15.3 

11.0 

3.0 

2.2 

9.3 

2.6 

7.4 

3.2 

Average 

9.4 

11.2 

12.5 

11.4 

12.3 

4.3 

5.3 

13.5 

4.8 

9.8 

7.7 

SO JS (40) 

7.7 

10.8 

15.1 

14.9 

14.4 

7.8 

9.0 

12.6 

9.8 

11.3 

3.1 

LN JS (40) 

50.3 

48.2 

65.8 

61.1 

40.0 

40.6 

26.2 

56.0 

34.5 

49.8 

45.3 

WP JS (40) 

6.5 

10.8 

16.1 

21.3 

19.4 

6.0 

5.7 

15.3 

5.9 

12.5 

6.5 

MN JS (40) 

4.9 

7.7 

11.5 

15.1 

10.6 

4.8 

3.8 

10.0 

4.3 

8.3 

6.0 

Average 

17.4 

19.4 

27.1 

28.1 

21.1 

14.8 

11.2 

23.5 

13.6 

20.5 

15.2 

SO GS (40) 

5.2 

6.9 

9.5 

9.3 

9.3 

6.2 

7.6 

8.0 

8.1 

7.7 

1.5 

LN GS (40) 

58.7 

51.5 

57.5 

55.3 

42.7 

38.2 

29.5 

56.5 

35.7 

51.3 

52.9 

WP GS (40) 

7.5 

12.1 

17.5 

24.3 

22.5 

7.2 

6.6 

17.3 

6.9 

14.1 

6.2 

MN GS (40) 

7.8 

9.5 

18.8 

23.5 

18.0 

5.4 

4.4 

15.5 

4.9 

12.5 

8.0 

Average 

19.8 

20.0 

25.8 

28.1 

23.1 

14.3 

12.0 

24.3 

13.9 

21.4 

17.1 

DO JG (50) 

29.7 

35.3 

34.4 

23.3 

20.0 

7.3 

9.7 

31.2 

8.5 

23.2 

23.6 

SC JG (50) 

3.2 

3.4 

7.7 

9.8 

7.9 

2.4 

2.0 

6.4 

2.2 

5.2 

3.3 

WC JG (50) 

8.3 

10.8 

13.1 

9.0 

7.4 

3.1 

6.6 

9.7 

4.9 

8.3 

9.8 

FF JG (50) 

28.2 

32.5 

32.7 

25.5 

21.8 

7.2 

6.7 

31.2 

7.0 

22.5 

20.2 

OT JG (50) 

56.9 

67.4 

64.6 

39.3 

24.3 

8.2 

10.9 

50.5 

9.6 

38.8 

44.6 

MT JG (50) 

7.1 

6.5 

11.6 

12.9 

16.2 

4.6 

4.2 

15.3 

4.4 

10.1 

4.8 

ET JG (50) 

2.5 

5.3 

4.5 

6.1 

5.6 

3.9 

6.5 

4.8 

5.2 

4.9 

4.7 


Table 7. A summary of the 0, $, and $ values for the computation of 


Comparison 

0(412) 

0(443) 

0(490) 

0(510) 

0(555) 

0(665) 

0(683) 

*1-5 

$6-7 

$1-7 

$3/5 

so 

JG (40) 

2.8 

3.8 

4.3 

4.2 

3.9 

23.8 

16.1 

3.8 

19.9 

8.4 

1.1 

LN 

JG (40) 

13.4 

12.4 

14.1 

13.0 

12.0 

36.4 

44.0 

13.0 

40.2 

20.8 

7.1 

WP 

JG (40) 

2.3 

2.9 

3.7 

3.9 

4.3 

3.5 

5.9 

3.4 

4.7 

3.8 

0.9 

MN 

JG (40) 

4.4 

5.1 

6.2 

6.7 

7.8 

6.0 

10.1 

6.0 

8.1 

6.6 

1.8 

Average 

5.7 

6.0 

7.1 

7.0 

7.0 

17.4 

19.0 

6.6 

18.2 

9.9 

2.7 

DO 

JG (50) 

9.6 

10.9 

12.2 

10.6 

9.3 

36.3 

39.1 

10.5 

37.7 

18.3 

6.0 

SC 

JG (50) 

3.3 

4.0 

4.9 

5.3 

6.1 

4.7 

8.0 

4.7 

6.4 

5.2 

| 1.4 

WC 

JG (50) 

3.8 

4.7 

5.2 

4.9 

4.4 

21.2 

19.9 

4.6 

20.6 

9.2 

2.0 

FF 

JG (50) 

9.9 

11.0 

12.8 

11.6 

11.0 

25.8 

32.4 

11.3 

29.1 

16.3 

6.0 

i 

OT 

JG (50) 

13.2 

15.6 

17.7 

14.6 

11.7 

45.8 

54.1 

14.6 

50.0 

24.7 

10.6 

MT 

JG (50) 

5.7 

6.2 

7.7 

7.8 

8.5 

17.3 

22.1 

7.2 

19.7 

10.8 

2.5 

ET 

JG (50) 

3.9 

4.3 

4.3 

4.2 

4.0 

13.8 

9.6 

4.2 

11.7 

6.3 

1.0 


of Ed{ 0",A). This point is developed more completely in 
Sect. 1.5.4 when the £<i(0 + , A) intercomparisons are made. 

1.5.3 K l and K ^ 

Diffuse attenuation coefficients are calculated from the 
vertical profiles of the in-water light field. The L u (z, A), 
Ed(z, A), and E u (z, A) data are used to calculate Ki{ A), 
Kd{ A), and K u ( A), respectively, which are needed to calcu- 
late (by linear extrapolation of the log-transformed data) 


the null depth values of the corresponding optical param- 
eters. Of the three attenuation coefficients, only Kd{ A) is 
intercompared for all three processors; Ki{\) is only in- 
tercompared for the JRC and GSFC processors, because 
K l ( A) is not normally reported by the Satlantic processor 
(although it is calculated and used). 

Given the importance of the attenuation coefficients to 
the derivation of the primary optical variables, it is im- 
portant they be estimated as accurately as possible. Sum- 
maries of the intercomparison of the AT(A) and Kl( A) 
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values are presented in Tables 6 and 7, respectively. The 
smallest differences are associated with the JRC and GSFC 
intercomparisons. For the JG comparisons, the winch and 
crane data (WP and SO) are the best, and the free-fall 
data the worst (the LN data are notably poor). For the 
JS and GS comparisons, there is a mixture of best results 
between SO and MN instrument types; the LN results are 
always the worst, and in some cases, significantly so. In 
all cases, the band-ratio results are better than the spec- 
tral averages, but they are not always better than some 
individual wavelengths. 

One quality control parameter that can be applied to 
the Kd{ A) estimates is to compare them to Kd{ A) values 
for pure water, K w ( A), computed from the Pope and Fry 
(1997) absorption data and the Morel (1974) backscatter 
values for pure water. f The intercomparison is made by 
calculating the percentage of casts for each processor and 
the deployment methods (WC and FF) for which K d (X) < 
K w { A) (in the original 40-cast data set), and the results are 
shown in Fig. 11 as a function of wavelength. Noting that 
one cast equates to 5% in each 20-cast category within the 
plot (FF or WC), the WC data satisfy the quality control 
test except for the Satlantic results beyond the blue part 
of the spectrum, for which an increasingly large number of 
casts do not. All of the processors have difficulty with the 
FF data, for which there is a general increase in failures 
from the blue to the red domain. Although not shown in 
the plot, the majority of the failures come from the LN 
data (as inferred in Table 6). 



Fig. 11. The percentage of casts for which Kd( A) < 
K w ( A) as a function of the data processors and de- 
ployment methods. The solid symbols and lines cor- 
respond to the FF category, and the open symbols 
and dashed lines correspond to the WC category. 


f ProSoft posts warnings in the header records of the ouput 
file when Kd < K w . 


1 . 5.4 E d { 0 + ) 

Changing atmospheric conditions, particularly under 
cloud cover, leads directly to variability of the in-water 
light field which must be corrected to obtain accurate es- 
timations of optical properties from irradiance or radiance 
profiles. First-order corrections for this variability can be 
made using above-water (on deck) measurements of so- 
lar irradiance, E d ( 0 + ,A). For the AMT and AAOT mea- 
surements, the emphasis is on collecting data only dur- 
ing clear-sky conditions, so the extra variance from at- 
mospheric variability during optical profiling is minimized. 
Regardless of whether or not the in- water optical data are 
corrected for atmospheric variability on a point-by-point 
basis (Sect. 3. 3. 2. 2), the E d ( 0 + , A) data are needed to pro- 
vide for the normalization of Lw{ A) in the calculation of 
-Rrs(A) and [Lvv^A)]^- 

Although not as notable as the results for E d { 0 , A) 
(Sect. 1.5.2), the E d { 0 + ,A) differences also partition into 
two groups, as shown in Table 8: a) those associated with 
a stable reference (WP, MN, and SC), and b) those asso- 
ciated with an unstabilized reference (SO, LN, and DO). 
Because all of the intercomparisons for the three proces- 
sors are excellent — always to within 2.5% it is easier to 
see this effect in the DO and SC intercomparisons for the 
JG extended data set. For these results, the unstabilized 
reference values are a factor of 5—6 larger with respect to 
the stabilized values (noting that this equates to a differ- 
ence of approximately 1.0%); the band-ratio results are, 
however, to within 0.1% in either case. All of the OT data 
are from unstabilized references, so these results are the 
poorest, with some uncertainties approaching 2.5%. 

As discussed earlier, a useful diagnostic as to the effi- 
cacy of the E d { 0“, A) calculation is the ratio between the 
two with the latter adjusted for the effects of the air-sea 
interface: E d {Q~ , A)/ [0.96 E d {0 + , A)] . A comparison of the 
UPD values for this ratio is presented in Table 9. For 
all processors, the lowest differences occur for the stable 
references, WP and MN (SC), and the largest differences 
occur for the unstabilized references, SO and LN (DO), 
with the latter reaching very high differences for the JS 
and GS comparisons. The variance associated with choos- 
ing the top and bottom extrapolation intervals is included 
in these results, so the true effect of an unstabilized ref- 
erence will not be known until the processing options are 
made more equal (Sects. 1.6.1 and 1.6.2). 

1 . 5.5 R rs 

Derivations of Lw( A) from measurements of L u (z, A) 
made in identical water conditions, but under different 
illumination conditions, will differ. This variability can 
be accounted for, in part, by adjusting Lw{ A) using the 
Ed{ 0 + ,A) measurement to derive the remote sensing re- 
flectance, R ts ( A) = L w (\)/E d (0 + ,\). This normalization 
of the illumination conditions makes R r s(A) the primary 
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Table 8. A summary of the 0. 4c and $ values for the computation of i0/(O + A). 


Comp 

a rison 

0(412) 

0( 443) 

0(490) 

0(510) 

-0(555) 

0(665) 

0(683) 

*1-5 

$ 6-7 

*1-7 

$3/5 

so 

JG (40) 

0.2 

0.2 

0.3 

0.3 

0.3 

0.3 

0.3 

0.2 

0.3 

0.3 

0.0 

LN 

JG (40) 

0.3 

0.7 

0.2 

0.6 

0.4 

0.6 

0.2 

0.4 

0.4 

0.4 

0.2 

VVP 

JG (40) 

0.3 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.0 

MN 

JG (40) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.0 

A verage 

0.2 

0.3 

0.2 

0.3 

0.3 

0.3 

0.2 

0.3 

0.3 

0.3 

0.1 

SO 

JS (40) 

1.0 

1.0 

1.0 

0.9 

1.0 

0.9 

0.9 

1.0 

0.9 

1.0 

0.1 

LN 

JS (40) 

0.5 

0.9 

1.3 

1.4 

1.9 

1.3 

1.2 

1.2 

1.2 

1.2 

0.6 

WP 

JS (40) 

0.7 

0.8 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.0 

MN 

JS (40) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.0 

A verage 

0.6 

0.7 

0.8 

0.8 

0.9 

0.8 

0.7 

0.8 

0.8 

0.8 

0.2 

so 

GS (40) 

0.9 

0.8 

0.8 

0.7 

0.8 

0.7 

0.7 

0.8 

0.7 

0.8 

0.0 

LN 

GS (40) 

0.8 

1.3 

1.2 

1.6 

1.9 

1.5 

1.0 

1.3 

1.3 

1.3 

0.8 

WP 

GS (40) 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.0 

MN 

GS (40) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.2 

0.2 

0.1 

0.2 

0.1 

0.0 

A verage 

0.5 

0.7 

0.6 

0.7 

0.8 

0.7 

0.6 

0.7 

0.6 

0.7 

0.2 

DO 

JG (50) 

1.0 

1.2 

1.1 

1.2 

1.2 

1.3 

1.2 

1.1 

1.2 

1.2 

0.1 

SC 

JG (50) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.0 

wc 

JG (50) 

1.2 

1.3 

1.4 

1.4 

1.5 

1.5 

1.5 

1.3 

1.5 

1.4 

0.1 

FF 

JG (50) 

0.2 

0.4 

0.2 

0.3 

0.3 

0.3 

0.2 

0.3 

0.3 

0.3 

0.1 

OT 

JG (50) 

1.8 

2.0 

2.2 

2.2 

2.3 

2.4 

2.4 

2.1 

2.4 

2.2 

0.1 

MT 

JG (50) 

0.2 

0.4 

0.2 

0.4 

0.3 

0.3 

0.2 

0.3 

0.3 

0.3 

0.1 

ET 

JG (50) 

0.3 

0.3 

0.3 

0.3 

0.4 

0.4 

0.4 

0.3 

0.4 

0.3 

0.0 


Table 9. A summary of the 0, and <f> values for the computation of 10/(0 , A)/ 0.96 Ed(0 + . A)] . 


Comparison 

0(412) 

CO 

0(490) 

0(510) 

0(555) 

0(665) 

0(683) 

*1-5 

*6-7 

*1-7 

*3/5 

SO 

JG (40) 

4.4 

4.8 

4.9 

5.1 

5.6 

8.7 

18.2 

5.0 

13.4 

7.4 

1.7 

LN 

JG (40) 

6.5 

8.5 

5.8 

6.2 

8.3 

20.5 

24.2 

7.1 

22.3 

11.4 

4.5 

WP 

JG (40) 

1.3 

1.1 

1.1 

1.3 

1.2 

1.7 

1.8 

1.2 

1.8 

1.4 

0.7 

MN 

JG (4°) 

1.7 

1.5 

2.2 

2.3 

1.8 

2.8 

1.9 

1.9 

2.3 

2.0 

1.1 

A verage 

3.5 

4.0 

3.5 

3.7 

4.2 

8.4 

11.5 

3.8 

10.0 

5.5 

2.0 

SO 

JS (40) 

10.5 

12.0 

12.5 

12.4 

12.7 

41.7 

43.9 

12.0 

47.7 

19.5 

2.3 

LN 

JS (40) 

22.5 

20.3 

21.1 

20.6 

22.0 

74.9 

65.6 

21.3 

68.0 

34.4 1 

6.9 

WP 

JS (40) 

4.2 

4.1 

4.1 

3.4 

5.4 

5.3 

4.8 

4.2 

5.0 

4.5 

3.5 

MN 

JS (40) 

4.8 

4.4 

3.2 

3.8 

3.1 

9.9 

10.4 

3.9 

10.1 

5.6 

1.5 

A verage 

10.5 

10.2 

10.2 

10.0 

10.8 

32.9 

31.2 

10.3 

32.7 

16.0 

3.6 

SO 

GS (40) 

7.6 

8.9 

9.7 

9.6 

9.9 

35.5 

34.5 

9.1 

40.0 

15.6 

1.9 

LN 

GS (40) 

26.9 

25.3 

23.5 

22.9 

28.6 

81.5 

77.5 

25.5 

78.2 

40.3 

10.6 

WP 

GS (40) 

3.8 

4.3 

4.1 

3.4 

5.6 

5.6 

5.7 

4.3 

5.6 

4.6 

3.2 

MN 

GS (40) 

5.4 

5.1 

4.2 

4.4 

4.7 

11.2 

11.0 

4.7 

11.1 

6.6 

2.3 

Average 

1 10.9 

10.9 

10.4 

10.1 

12.2 

33.4 

32.2 

10.9 

33.7 

16.8 

4.5 

DO 

JG (50) 

7.2 

8.4 

7.8 

7.9 

9.0 

15.6 

22.1 

8.1 

18.8 

11.1 

2.6 

SC 

JG (50) 

1.5 

1.3 

1.7 

1.8 

1.5 

2.2 

1.9 

1.6 

2.1 

1.7 

0.9 

WC 

JG (50) 

j 3.7 

4.0 

4.2 

4.3 

4.5 

5.7 

13.3 

4.2 

9.5 

5.7 

1.1 

FF 

JG (50) 

5.9 

6.9 

6.3 

6.4 

7.3 

14.1 

14.6 

6.6 

14.3 

8.8 

2.6 

OT 

JG (50) 

10.0 

11.2 

11.5 

10.8 

11.6 

16.4 

22.3 

11.0 

19.4 

13.4 

1.6 

MT 

JG (50) 

! 3.2 

3.1 

3.5 

3.9 

4.2 

8.5 

8.6 

3.6 

8.6 

5.0 

2.1 

ET 

JG (50) 

3.0 

4.4 

2.8 

3.1 

3.8 

7.4 

14.6 

3.4 

11.0 

5.6 

1.9 
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variable for estimating chlorophyll a concentration from in 
situ optical measurements. The OC2v2 algorithm, for ex- 
ample, uses log 10 [F rs (490)/i? rs (555)] as the input vari- 
able. In addition, the remote sensing reflectance is one of 
the primary variables relating the AOPs of seawater to the 
inherent optical properties (IOPs). 

The excellent agreement between the processors in the 
calculation of F^(0 + , A) (Sect. 1.5.4) suggests much of the 
variability in the derivation of F rs (A) will be due to the dif- 
ferences in calculating Ljy(A). Because the water-leaving 
radiance is computed as a simple constant times the up- 
welled radiance at null depth, Lw( A) = 0.54L U (0~A), the 
differences in the L u (0‘A) comparisons will determine most 
of the jRrs(A) differences. The exceptions will be when the 
differences in Ed( 0 + ,A) were maximal, which will increase 
the differences above the levels found for L u (0"A). 

The uncertainties for the R rs (\) intercomparisons are 
given in Table 10. As anticipated, the levels of agreement 
closely follow the values established for L u ( 0", A) (Table 4): 

1. The smallest differences are associated with the JG 
intercomparisons (the JS and GS intercomparisons 
are, on average, a little more than twice as large); 

2. The maximum differences are in the red part of the 
spectrum (665 nm), with minimal differences in the 
blue-green (490-510 nm) which increase towards the 
blue (412-443 nm); 

3. The band-ratio results are always the smallest dif- 
ferences (except for some of the Satlantic intercom- 
parisons); 

4. The shallow coastal data (particulary the WiSPER 
data) have smaller differences than the deep ocean 
data (most notably in the red wavelengths); and 

5. The results as a function of increasing chlorophyll a 
concentration, are not as simply defined as before, 
because of the elevated differences in Ed( 0 + ,A) for 
the OT regime. 

The 2.5% intercomparison objective is satisfied for nu- 
merous comparisons noted earlier. For comparisons with 
the Satlantic processor, only the blue-green WiSPER and 
some of the band-ratio uncertainties are below 2.5% — most 
of the results are, on average, significantly above the thresh- 
old. For the JG comparisons, all of the WiSPER, and some 
of the miniNESS differences, are below 2.5%, as are all of 
the band-ratio results. Many of the shallow coastal results 
satisfy the objective, but the winch and crane average re- 
sults do not (except the band-ratio results), because of the 
elevated differences in the estimation of Ed( 0 + ,A). Aside 
from the band-ratio results, the only concentration regime 
with uncertainties below 2.5% is the mesotrophic regime 
(but only for the blue-green wavelengths). 

1.5.6 [£if] N 

The so-called normalized water-leaving radiance (Gor- 
don and Clark 1981), [L w (A)] N , is defined as the hypo- 
thetical water-leaving radiance that would be measured in 


tbe absence of any atmospheric loss with a zenith sun at 
the mean Earth-sun distance. The latter is accomplished 
by adjusting R r S (A) with the (time dependent) mean ex- 
traterrestrial solar irradiance, Fo(A,gQ, which is usually 
formulated to depend on the sequential day of the year 
(SDY), d. 

The intercomparison of normalized water-leaving radi- 
ances for the different processors and data categories is 
presented in Table 11. Because [Lvv(A)] n is calculated 
directly from F rs (A), the uncertainties for the former are 
expected to follow those of the latter (Table 10). Although 
there is general agreement in the patterns between the 
two, in terms of the spectral shapes and the relationships 
between the analysis categories (e.g., the concentration 
regimes), the magnitude of the uncertainties differ signifi- 
cantly, with the exception of the JG band-ratio results. 

The most likely explanation for the increase in un- 
certainties is in the chosen values for the Fo(A) correc- 
tion term: the JRC processor uses the Neckel and Labs 
(1984) values interpolated to specific center wavelengths 
(Sect. 3. 3. 2. 8), while the GSFC processor uses the Neckel 
and Labs (1984) values weighted by the spectral response 
of the SeaWiFS bands (Sect. 2.4.4). Figure 12 shows the 
relationship between the GSFC and JRC F 0 (A) values (the 
GSFC processor does not include a value for 683 nm, be- 
cause this is not a SeaWiFS wavelength). There is a strong 
bias between the two processors with maximum differences 
in the blue (6.5% at 412 nm), intermediate differences in 
the green (3.2% at 510 nm), and minimum differences in 
the red (1.9% at 665 nm). The band-ratio results for the 
JRC and GSFC intercomparison of F 0 (A) are very good — 
the maximum percent difference is 0.4%. The JG R r S (A) 
band ratios are all less than 2.5%, so the excellent F 0 (A) 
band ratios ensure very good [Lw(A)] n band ratios. 



JRCFo(A) [mW crn-2 

Fig. 12. A comparison of the JRC and GSFC de- 
terminations of Fq(A) (for the 50-cast data set). 
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Table 10. A summary of the 0, 'I', and $ values for the computation of R rs 


Comparison 

0(412) 

0(443) 

0(490) 

0(510) 

0(555) 

^(665) 

0(683) 

'I' 1-5 

$6-7 

01-7 

$3/5 

so 

JG (40) 

3.8 

4.1 

3.4 

3.2 

2.8 

26.8 

12.4 

3.5 

19.6 

8.1 

1.0 

LN 

JG (40) 

3.7 

3.6 

2.6 

2.5 

4.0 

25.1 

23.8 

3.2 

24.8 

9.1 

2.1 

WP 

JG (40) 

1.7 

1.7 

1.5 

1.5 

1.5 

2.3 

3.5 

1.6 

2.9 

1.9 

0.2 

MN 

JG (40) 

3.0 

2.8 

2.3 

2.3 

2.4 

6.0 

8.1 

2.6 

7.0 

3.8 

0.4 

Average 

3.0 

3.0 

2.5 

2.4 

2.7 

15.0 

11.9 

2.7 

13.6 

5.7 

0.9 

so 

JS (40) 

13.4 

13.2 

10.0 

9.6 

8.6 

46.3 

32.8 

10.9 

40.9 

18.0 

2.0 

LN 

JS (40) 

7.0 

7.0 

6.5 

6.7 

8.4 

44.7 

43.5 

7.1 

42.5 

17.1 

7.2 

WP 

JS (40) 

2.0 

1.8 

1.6 

1.6 

1.7 

3.3 

4.5 

1.8 

3.9 

2.4 

0.2 

MN 

JS (40) 

7.5 

6.2 

4.6 

4.5 

4.3 

11.3 

10.7 

5.4 

11.0 

7.0 

0.3 

A verage 

7.5 

7.1 

5.7 

5.6 

5.8 

26.4 

22.9 

6.3 

24.6 

11.1 

2.4 

SO 

GS (40) 

9.9 

10.4 

8.0 

7.4 

7.2 

28.8 

25.8 

8.6 

27.1 

13.0 

1.6 

LN 

GS (40) 

7.8 

8.3 

7.3 

7.3 

10.5 

47.3 

48.3 

8.2 

46.4 

19.0 

8.2 

WP 

GS (40) 

2.3 

2.1 

1.8 

1.8 

1.8 

4.2 

4.7 

2.0 

4.4 

2.7 

0.3 

MN 

GS (40) 

7.4 

5.8 

3.9 

3.7 

3.5 

10.9 

10.2 

4.9 

10.6 

6.5 

0.6 

Average 

6.9 

6.7 

5.2 

5.1 

5.7 

22.8 

22.2 

5.9 

22.1 

10.3 

2.7 

DO 

JG (50) 

3.5 

3.6 

3.1 

3.2 

3.9 

25.8 

21.6 

3.5 

23.8 

9.2 

1.7 

SC 

JG (50) 

2.3 

2.2 

1.9 

1.9 

2.0 

4.1 

5.8 

2.1 

4.9 

2.9 

0.3 

WC 

JG (50) 

3.3 

3.5 

3.1 

3.1 

3.0 

17.2 

12.1 

3.2 

14.6 

6.5 

0.7 

FF 

JG (50) 

2.9 

2.7 

2.2 

2.4 

3.3 

16.7 

18.0 

2.7 

17.7 

6.9 

1.5 

OT 

JG (50) 

2.9 

2.9 

3.1 

3.7 

4.6 

24.0 

26.8 

3.4 

25.4 

9.7 

1.8 

MT 

JG (50) 

2.4 

2.5 

2.1 

2.0 

2.7 

13.3 

13.8 

2.3 

13.9 

5.5 

0.9 

ET 

JG (50) 

4.1 

4.1 

3.1 

3.0 

2.6 

15.9 

7.6 

3.4 

11.8 

5.8 

0.7 


Table 11. A summary of the ip-. 
Comparison •0(412) 0(443) 

SO JG (40) 8.0 5.2 

LN JG (40) 6.1 4.8 

WP JG (40) 7.2 4.5 

MN JG (40) 8.2 5.6 

Average 7.4 5.0 

SO JS (40) 15.0 13.8 

LN JS (40) 10.8 8.8 

WP JS (40) 4.7 4.0 

MN JS (40) 10.9 10.3 

Average 10.3 9.2 

SO GS (40) 8.6 9.6 

LN GS (40) 7.6 8.3 

WP GS (40) 3.6 2.7 

MN GS (40) 6.3 6.4 

Average 6.5 6.7 

DO JG (50) 7.0 4.9 

SC JG (50) 7.7 5.0 

WC JG (50) 7.6 5.1 

FF JG (50) 7.0 4.8 

OT JG (50) 6.3 4.2 

MT JG (50) 7.0 4.9 

ET JG (50) 8.5 5.7 


0, and $ values for the computation of \L\ 
0(490) 0(510) 0(555) 0(665) 0(683) 


8.0 

5.2 

3.9 

5.6 

3.2 

29.0 

6.1 

4.8 

3.5 

4.9 

4.3 

27.1 

7.2 

4.5 

3.5 

4.0 

3.6 

2.9 

8.2 

5.6 

4.8 

4.5 

4.8 

8.1 

7.4 

5.0 

3.9 

4.7 

4.0 

16.8 

15.0 

13.8 

10.8 

9.6 

8.5 

45.4 

10.8 

8.8 

8.4 

7.3 

9.4 

44.9 

4.7 

4.0 

4.0 

2.1 

2.6 

4.2 

10.9 

10.3 

8.3 

3.6 

5.3 

9.8 

10.3 

9.2 

7.9 

5.7 

6.5 

26.1 

8.6 

9.6 

7.7 

6.4 

6.3 

27.8 

7.6 

8.3 

7.5 

6.7 

10.2 

47.9 

3.6 

2.7 

2.3 

4.1 

2.9 

5.9 

6.3 

6.4 

4.7 

3.8 

3.7 

10.2 

6.5 

6.7 

5.6 

5.3 

5.8 

23.0 

7.0 

4.9 

3.9 

5.7 

4.5 

27.9 

7.7 

5.0 

4.2 

4.2 

4.2 

5.5 

7.6 

5.1 

4.1 

5.2 

4.0 

19.1 

7.0 

4.8 

3.9 

5.0 

4.7 

18.5 

6.3 

4.2 

3.9 

6.0 

5.1 

26.2 

7.0 

4.9 

3.9 

4.5 

4.3 

15.3 

8.5 

5.7 

4.2 

5.2 

3.8 

17.3 


w(A)] N - 

) frl-5 
5.2 
4.7 

4.6 

5.6 
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All of the individual spectral differences in Fig. 12 lie 
o^ straight lines, which suggests there is a deterministic, 
at significant, difference in how the Fq( A) terms are cho- 
sen (this is also inferred by the band-ratio results). A 
comparison between the Satlantic and JRC processor is, 
however, somewhat different. The band-ratio results are 
all elevated, the overall average is 2.2%, and the maxi- 
mum individual wavelength differences range from 7.7% at 
412 nm, 3.2% at 510 nm, and 2.8% at 665 nm. When these 
are combined with the R r S (A) differences, the [Lh/(A)] n 
uncertainties are all larger. Very similar results are seen 
for the Satlantic and GSFC comparisons. 

Very few of the [Lvt/(A)] N uncertainties are less than 
2.5%. The only notable exceptions are the JG band-ratio 
results, and a few of the WP results. The range of differ- 
ences is frequently small for the JG results, e.g., the DO, 
SC, WC, and FF results at 412 nm range from 7. 0-7. 7%, so 
if the Fo (A) calculations can be reconciled, there is every 
reason to believe these uncertainties can be significantly re- 
duced. The latter is only partly true for the Satlantic com- 
parisons, because the R rs ( A) uncertainties were frequently 
large, but there is certainly a chance to improve the WP 
intercomparisons. 

1.5.7 -Ejtr) and Q n ( O') 

L u { 0”,A) and £ U (0",A) are related by the Q-factor. 
This permits an alternative approach for the calculation of 
Lw{ A), based on applying the Q-factor to measurements of 
the irradiance reflectance, denoted R(z, A) and defined as 
E u {z, X)/E d (z, A). The uncertainties associated with this 
approach are only explored to the extent that intercom- 
parisons of E u ( 0”,A), R(0~ i A), and Q n (0~iA) are made 
for the JRC and GSFC processors (the Satlantic processor 
does not compute the Q-factor). 

The E u ( 0“,A) uncertainties are presented in Table 12. 
The best agreement occurs for the shallow coastal data 
(WP and MN), which have overall uncertainties below 

2.5% in the blue-green part of the spectrum, slightly above 
2.5% in the red domain, and band ratios well below 2.5%. 
The deep ocean data are a mixture of low (LN) and high 
(SO) uncertainties, with universally high values in the red 
wavelengths. The poorer results for SO data are hard to 
explain, because this instrument produces more data per 
unit depth interval than LN. It is interesting to note a 
similar result for the L u (0",A) results (Table 4), although 
it is not as significant. In comparison, the JG results 
for E d ( 0“,A), show the SO uncertainties are substantially 
lower than the LN uncertainties (Table 5). 

The (unexpectedly) larger uncertainty with the SO data 
is the reason the FF category has a smaller average uncer- 
tainty than the WC category, and it produces a clear trend 
of increasing uncertainty with increasing chlorophyll a con- 
centration (the SO data represent a significant portion of 
the ET data, 60% of all casts, as shown in Table 3). This 
is a clear example of how biases in sampling (and then 


processing results) can influence the relationships between 
the different analysis categories. 

The intercomparison results for determining R( A) are 
presented in Table 13. These are necessarily a combination 
of the E u ( CT, A) and E d (0 ~ , A) results, and there is a gen- 
eral elevation in uncertainties, even in the band-ratio re- 
sults. Although the increase in uncertainties is sometimes 
small, there are some notable changes in the relationships 
among the categories: a) the SO results are now lower than 
the LN results, so the WC uncertainties are less than the 
FF uncertainties; and b) the uncertainties as a function of 
concentration regime are more uniform. 

The Q-factor is the ratio between the upwelled irradi- 
ance and radiance, so intercomparisons of Q n (CT, A), which 
are presented in Table 14, will be a combination of Ta- 
bles 12 and 4, respectively. The patterns of agreement 
for Q n ( 0~,A) closely follow the patterns established for 
E u ( 0~,A) except there is greater uniformity between the 
instrument types (the WC and FF results are very simi- 
lar), there are much lower uncertainties in the red part of 
the spectrum, and there is less distinction as a function of 
concentration regime (the still anomalously high SO un- 
certainties elevate the ET results). 

Out of all the in-water variables, the overall uncertain- 
ties for Q n (CT, A) are the lowest, and many of the various 
categories of uncertainty are to within 2.5%. This is a 
direct consequence of the two input variables, L u ( CT,A) 
and E u (0~ , A), being upwelled measurements, so they are 
largely free of any surface effects (which raise the uncer- 
tainties in the measurements). 

1.5.8 Ancillary Variables 

Variables not immediately associated with the optical 
fields, but important to the derivation of optical parame- 
ters, plus those variables by which optical parameters are 
evaluated are important to the understanding of a total un- 
certainty budget. Examples of such variables are the solar 

zenith angle, 9 , and the total chlorophyll a concentration 
derived from the OC2v2 algorithm. 

The Satlantic data submission did not include the so- 
lar zenith angle, so only JRC and GSFC 9 values were 
intercompared. For all of the standard (40-cast) data, the 
uncertainty in 0 had an overall average value of 0.7%; the 
minimum-to-maximum differences for these data were 0.0- 
2.3%. For the 10 casts used to extend the standard data 
set, all of which were deep casts in oligotrophic conditions, 
the average uncertainty in 9 was 1.1%, and the minimum- 
to-maximum differences were 0. 3-5.0% (most of the vari- 
ance was attributed to slightly different reference times 
used in the computations). 

The intercomparison of the OC2v2 determination of to- 
tal chlorophyll a concentration between the three data pro- 
cessors and with the high performance liquid chromatog- 
raphy (HPLC) determination of total chlorophyll a is pre- 
sented in Table 15. The most obvious result is that the 
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Table 12. A summary of the 0 . T. and 4> values for the computation of £0(0". A). 


Comparison 

■0(412) 

0(443) 

0(490) 

0(510) 

0(555) 

0(665) 

0(683) 

*1-5 

*6-7 

*1-7 

*3/5 

so 

JG (40) 

4.4 

4.7 

3.9 

4.2 

4.6 

30.2 

12.3 

4.4 

21.3 

9.2 

1.4 

LN 

JG (40) 

1.7 

1.8 

1.9 

1.9 

2.0 

29.0 

27.0 

1.8 

28.0 

9.3 

0.8 

WP 

JG (40) 

1.1 

1.1 

0.9 

0.9 

1.0 

2.1 

3.1 

1.0 

2.6 

1.5 

0.2 

MN 

JG (40) 

2.5 

2.6 

2.0 

2.1 

2.2 

3.8 

3.6 

2.3 

3.7 

2.7 

0.3 

Average 

2.4 

2.5 

2.2 

2.3 

2.4 

16.3 

11.5 

2.4 

13.9 

5.7 

0.7 

DO 

JG (50) 

2.3 

2.3 

2.1 

2.4 

2.9 

32.0 

25.6 

2.4 

28.8 

9.9 

1.4 

SC 

JG (50) 

1.8 

1.8 

1.5 

1.5 

1.6 

2.9 

3.4 

1.6 

3.2 

2.1 

0.2 

wc 

JG (50) 

2.5 

2.6 

2.2 

2.3 

2.8 

22.4 

15.4 

2.5 

18.9 

7.2 

1.1 

FF 

JG (50) 

1.7 

1.7 

1.6 

1.7 

2.0 

18.6 

17.8 

1.8 

18.2 

6.5 

0.8 

OT 

JG (50) 

0.6 

0.5 

0.6 

0.9 

1.9 

34.0 

34.8 

0.9 

34.4 

10.5 

1.7 

MT 

JG (50) 

2.1 

2.1 

1.8 

1.9 

2.1 

14.7 

13.0 

2.0 

13.8 

5.4 

i 0.4 

ET 

JG (50) 

3.3 

3.7 

3.0 

3.1 

3.3 

16.8 

6.4 

3.3 

11.6 

5.7 

1.0 


Table 13. A summary of the 'll), T. and 4> values for the computation of R( 0 .A). 


Comparison 

0(412) 

0(443) 

0(490) 

0(510) 

0(555) 

'0(665) 

OO 

CO 

*1-5 

4/6-7 

*1-7 

*3/5 

so 

JG (40) 

3.2 

3.4 

3.3 

4.1 

4.6 

21.8 

21.4 

3.7 

21.6 

8.8 

2.1 

LN 

JG (40) 

7.6 

9.1 

7.0 

7.6 

9.0 

15.2 

26.1 

8.0 

20.7 

11.6 

3.7 

WP 

JG (40) 

1.3 

1.3 

1.1 

1.4 

1.1 

2.4 

3.3 

1.2 

2.8 

1.7 

0.6 

MN 

JG (40) 

2.0 

2.3 

3.2 

3.4 

3.0 

5.0 

3.9 

2.8 

4.4 

3.2 

1.2 

Average 

3.5 

4.0 

3.7 

4.1 

4.4 

in 

13.7 

3.9 

12.4 

6.4 

1.9 

DO 

JG (50) 

3.8 

4.3 

3.6 

4.2 

5.2 

24.5 

28.3 

4.2 

26.4 

10.6 

2.6 

sc 

JG (50) 

1.7 

1.8 

2.2 

2.4 

2.0 

3.7 

3.6 

2.0 

3.6 

2.5 

0.9 

wc 

JG (50) 

2.1 

2.1 

2.0 

2.5 

2.8 

18.8 

19.4 

2.3 

19.1 

7.1 

1.6 

FF 

JG (50) 

j 3.7 

4.4 

3.9 

4.3 

4.9 

14.0 

17.5 

4.2 

15.8 

7.5 

2.2 

OT 

JG (50) 

2.3 

2.5 

2.2 

2.2 

3.2 

32.4 

36.1 

2.5 

34.2 

11.6 

1.9 

MT 

JG (50) 

3.9 

3.6 

4.1 

4.7 

4.6 

10.5 

11.3 

4.2 

10.9 

6.1 

1.9 

ET 

JG (50) 

2T 

3.6 

2.2 

2.8 

3.6 

10.5 

13.6 

2.9 

12.0 

5.5 

2.1 


Table 14. A summary of the i b, T, and 4> values for the computation of Q n (0 .A). 


Comparison 

0(412) 

0(443) 

0(490) 

0(510) 

0(555) 

0(665) 

0(683) 

*1-5 

4> 6 -7 

*1-7 

*3/5 

so 

JG (40) 

3.6 

3.6 

3.5 

4.1 

4.1 

6.2 

9.1 

3.8 

7.6 

4.9 

0.7 

LN 

JG (40) 

2.6 

2.2 

2.0 

2.4 

2.9 

6.2 

6.4 

2.4 

6.3 

3.5 

1.4 

WP 

JG (40) 

1.0 

0.9 

0.7 

0.7 

0.7 

1.2 

1.3 

0.8 

1.2 

0.9 

0.1 

MN 

JG (40) 

1.6 

1.2 

1.2 

1.4 

1.6 

3.7 

4.8 

1.4 

4.3 

2.2 

0.4 

Average 

2.2 

2.0 

1.9 

2.1 

2.3 

4.3 

5.4 

2.1 

4.9 

2.9 

0.7 

DO 

JG (50) 

2.5 

2.5 

2.5 

2.9 

3.3 

7.5 

9.5 

2.7 

8.5 

4.4 

1.1 

sc 

JG (50) 

1.3 

1.1 

1.0 

1.1 

1.2 

2.5 

3.0 

1.1 

2.7 

1.6 

0.3 

wc 

JG (50) 

2.1 

2.1 

2.0 

2.3 

2.4 

5.7 

7.3 

2.2 

6.5 

3.4 

0.6 

FF 

JG (50) 

2.0 

1.8 

1.8 

2.0 

2.5 

5.3 

6.6 

2.0 

6.0 

3.1 

0.9 

OT 

JG (50) 

1.6 

1.5 

1.7 

2.0 

2.7 

9.3 

11.5 

1.9 

10.4 

4.3 

1.3 

MT 

JG (50) 

1.6 

1.6 

1.6 

1.8 

2.1 

4.6 

5.2 

1.7 

4.9 

2.6 

0.7 

ET 

JG (50) 

3.0 

2.7 

2.5 

2.8 

2.7 

3.5 

5.6 

2.7 

4.6 

3.3 

0.5 





S. Hooker, G. Zibordi, J-F. Berthon, D. D’Alimonte, S. Maritorena, S. McLean, and J. Sildam 


OC2v2 algorithm does a poor job and this is irrespective 
of the data processor being considered. Nonetheless, there 
are substantial differences between the data types, with 
the winch and crane results proving superior to the cor- 
responding free-fall results, i.e., the SO results are better 
than the LN results, and the WP results are superior to 
the MN results. 

Table 15. The intercomparison of the OC2v2 algo- 
rithm results between processors and with respect 
to the HPLC determination of Cr a f° r the standard 
(40-cast) data set. 


Data 

Type 

OC2v2 

HPLC 

JG 

JS 

GS 

j 

G 

s 

SO 

4.3 

4.3 

5.7 

68.2 

68.3 

67.8 

- LN 

6.9 

26.4 

29.2 

86.1 

80.4 

95.2 

WP 

5.0 

9.2 

12.0 

95.7 

97.1 

106.3 

MN 

6.9 

2.6 

9.2 

110.3 

109.7 

110.3 

Average 

5.8 

10.5 

14.1 

90.1 

88.9 

94.9 


It is also clear that the deep ocean results (SO and LN) 
are better than the shallow coastal results (WP and MN). 
This is expected, because the OC2v2 algorithm was formu- 
lated from Case-1 (deep ocean) data (O’Reilly et al. 1998), 
and, in general, is not suitable for coastal waters. The un- 
certainty in the OC2v2 algorithm, based on the original 
data set used to create it, is approximately 30.8%. If the 
AMT data set is used as a separate control to evaluate this 
uncertainty, the uncertainty in the OC2v2 algorithm based 
on the full AMT matrix (not just the data subset used for 
DARR-00) is approximately 34.4% when no filtering for 
possibly spurious observations is applied to the data; if 24 
stations, which appear as outliers are removed from the 
analysis (these usually correspond to unusual phytoplank- 
ton assemblages, highly heterogeneous water masses, etc.), 
the uncertainty is about 30.2%. 

The good agreement between the original uncertainty 
estimated of the OC2v2 algorithm and that achieved with 
the AMT data matrix suggests the large uncertainties seen 
with the SO and LN data in Table 15 are anomalously high 
and a consequence of dealing with a small subsample from 
a much larger database. 

1.5.9 VI Summary 

The standard results are hereafter referred to a s the 
VI results, and processing options not in keeping with the 
VI options are referred to as nonstandard. A summary of 
the VI processing uncertainties is presented in Table 16. 
The results are organized according to the primary vari- 
ables and the three concentration regimes. The latter are 
included to investigate whether or not there is any depen- 
dence with respect to chlorophyll a concentration as a func- 
tion of the processing options and, of course, the individ- 
ual variables. It is important to remember the results for 
the concentration regimes are derived from the extended 


(50-cast) data set and only involve the JRC and GSFC 
processors, whereas the individual parameter results are 
from the standard (40-cast) data set and all available pro- 
cessor combinations (the data sets involved are noted in 
the table alongside the two-letter processor codes). The 
principal VI results are as follows (Tables 4-14 and 16): 

1. The smallest uncertainties are associated with the 
JG intercomparisons (which, on average, are signif- 
icantly less than the .3 S or GS results). 

2. Many of the inter comparisons have maximum un- 
certainties in the red part of the spectrum (665 nm), 
with minimal uncertainties in the blue-green wave- 
lengths (490-510 nm) which increase towards the 
blue domain (412-443 nm). 

3. The band-ratio uncertainties are almost always the 
smallest, and are frequently less than 2.5%. 

4. Most of the uncertainties exceed the 2.5% inter- 
comparison objective, but the JG uncertainties for 
L u ( 0“), E d ( 0 + ), R rs , E u ( 0"), and Q n ( 0") are less 
than or close to 2.5% for many of the blue-green 
wavelengths. 

5. The shallow coastal data (in particular the WiS- 
PER data) have smaller uncertainties than the deep 
ocean data (most notably for the red wavelengths). 

6. The winch and crane data (especially the WiSPER 
data) have smaller uncertainties than the free-fall 
data (most notably in the red part of the spectrum). 

7. The uncertainties in the diffuse attenuation coeffi- 
cients, Kl and K d) exceed the uncertainties in the 
corresponding extrapolated (null depth) radiomet- 
ric measurements, L u ( 0") and E d ( 0"), respectively. 

8. As a function of increasing chlorophyll a concentra- 
tion (OT to MT to ET), there is a general increase 
in uncertainties in the blue-green part of the spec- 
trum, but a decrease in the red wavelengths, and 
a decrease in the band-ratio uncertainties (note the 
very large uncertainties for the K d OT results). 

In addition to following the evolution in the aforemen- 
tioned principal results, the subsequent analyses are con- 
cerned with the following objectives (documented below in 
Sects. 1.6. 1-1. 6. 4): 

a) Determine if results converge (i.e., the uncertainties 
approach zero) as the JRC and GSFC processing 
options are made more similar; 

b) Quantify the level of uncertainty associated with 
the processing options; and 

c) Establish which of the final uncertainties for the 
data products are within the 1.0% or 2.5% objec- 
tives set at the beginning of the DARR-00 activity. 

As an easily referenced indicator of the convergence pro- 
cess, the final row in Table 16 gives the averages of all the 
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Table 16 . A summary of the 0. T. and $ values for the computation of the VI primary variables. 


Comparison 

t/>(412) 

-0(443) 

0(490) 

-0(510) 

*0(555) 

*0(665) 

■0(683) 

*1-5 

*6-7 

*1-7 

*3/5 

MO') 
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1.6 
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4.6 
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0.7 

ET 
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3.0 

2.7 
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2.7 
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5.6 
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4.6 
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Table 17 . The GSFC and JRC processing options as a function of the processing version number. The selection 
of the extrapolation intervals is indicated with a one-letter code: G for GSFC and J for JRC, i.e., a G indicates 
the processing interval was set by the configuration for the GSFC data processor, and a J means it was set by 
the configuration for the JRC processor. 


Ver- 

sion 

No. 

Extrapolation Intervals 

Tilt 

Filtering 

E d ( 0 + ) 
Normalized 

Outlier 

Filtering 

Fitting 

Function 

so 

GSFC 
LN WP 

MN 

SO 

JRC 

LN WP 

MN 

GSFC 

JRC 

GSFC 

JRC 

GSFC 

JRC 

GSFC JRC 

1 

G 

G 

G 

G 

J 

J 

J 

j 

No 

Yes 

No 

Yes 

No 

Yes 

Linear Linear 

2 

G 

G 

J 

J 

G 

G 

J 

j 

No 

Yes 

No 

Yes 

No 

Yes 

Linear Linear 

3 

G 

G 

J 

J 

G 

G 

J 

j 

Yes 

Yes 

No 

No 

No 

Yes 

Linear Linear 

4 

G 

G 

J 

J 

G 

G 

J 

j 

Yes 

Yes 

No 

No 

Yesj 

Yes 

Robust Linear 

5 

G 

G 

J 

J 

G 

G 

J 

j 

No 

Yes 

No 

No 

No 

No 

Linear Linear 


t Partial outlier rejection from the robust fitting function which reduces the weight of data points with a large 


standard deviation. 


primary variables presented for the standard JG compar- 
isons, i.e., only the 40-cast data set is used. The conver- 
gence in the extended (50-cast) data set is discerned from 
the concentration regime results. 

1.6 RESULTS-NONSTANDARD 

Nonstandard options to the JRC and GSFC processors 
were executed to investigate differences in the VI inter- 
comparison by processing profile data with nonstandard 
options. The results from the subsequent changes to the 
standard processing options are identified as follows: 

V2 For the JRC processor, the only change with re- 
spect to the standard processing options was to 
use the extrapolation intervals for the SeaOPS and 
LoCNESS data as defined by the GSFC processing. 
Similarly, for the GSFC processor, the only change 
was to use the extrapolation intervals for the WiS- 
PER and miniNESS data as defined by the JRC 
processing. No changes were made to the Satlantic 
processing configuration. 

V3 In addition to the V2 changes, no E d { 0 + ) normal- 
ization was used with the JRC processor (the JRC 
processor still filtered the data for outliers and the 
GSFC processor did not), and the GSFC proces- 
sor used a 5° tilt filter (as did the JRC processor). 
No changes were made to the Satlantic processing 
configuration. 

The JRC and GSFC processors were very similar for the 
V3 configurations— the only difference between the two 
was the JRC processor used outlier filtering and the GSFC 
processor did not (which potentially affects a small amount 
of data in each cast). These options were then used as a 
baseline to explore two different aspects of the processing: 

V 4 In addition to the V3 changes, the GSFC processor 
was changed to use a robust linear fit, which weights 
the data in the extrapolation interval to minimize 
the significance of data points with higher variance. 


No changes were made to the Satlantic processing 
configuration. 

V 5 In addition to the V3 changes, the GSFC processor 
did not use a tilt filter, and the JRC processor did 
not filter the data for outliers. No changes were 
made to the Satlantic processing configuration. 

A summary of the options used with the GSFC and JRC 
processors for each processing version is given in Table 17. 

1.6.1 V2 Summary 

The only change made for the V2 processing was to use 
the same extrapolation intervals for the JRC and GSFC 
processors. A summary of the effects of this change is 
shown in Table 18. In comparison to the VI results (Ta- 
ble 16), there are substantial improvements in most of the 
JG variables. No changes were made to the Satlantic pro- 
cessor, and the JS and GS uncertainties are largely un- 
changed. The most significant aspects of the improvements 
are the uncertainties are lower at all wavelengths and they 
are spectrally flatter. Given the largest uncertainties for 
the VI results are in the red domain, the largest improve- 
ments with the V 2 processing are in the red domain. 

Because the JG values for the V2 results are spec- 
trally more uniform, * 1 - 5 , * 6 - 7 , and $ 3/5 are appropriate 
variables for summarizing the improvements derived from 
identical extrapolation intervals (the [ symbol represents 
a reduction in uncertainty with respect to the VI values): 


Lu(O-) 

* 1-5 

O 

II 

1.9)%, 

*6-7 " 

= 1.7(111.2)%, 

and 


$ 3/5 

= 0.1(1 

0.8)%; 





k l 

¥1-5 

= 0.8(| 

5.8)%, 

#6-7 = 

= 2.0(116.2)%, 

and 


$ 3/5 

= 0.4(| 

2.3)%; 





E d ( 0") 

* 1-5 

= 1.2(1 

2.5)%, 

*6-7 = 

= 2.5(1 

7.4)%, 

and 


$ 3/5 

= 0.3(1 

1.6)%; 





K d 

*1-5 

= 3.2(110.3)%, 

*6-7 = 

= 0.8(1 

4.0)%, 

and 


$ 3/5 

= 1-6(1 

6.1)%; 





E d ( 0 + ) 

* 1-5 

= 0.3() 

0.0)%, 

*6-7 “ 

= 0.3(1 

0.0)%, 

and 


$ 3/5 

- 0.1(1 

0.0)%; 
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Table 18 . A summary of the 0. T, and T values for the computation of the V2 primary variables. 
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11.9 

23.6 

13.7 

20.8 

16.6 


GS (40) 

18.9 

19.4 

23.8 

25.3 

20.7 

13.9 
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OT 
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0.8 
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JG (50) 

2.3 

3.0 

3.4 

3.2 

2.9 

1.0 

0.8 

3.0 

0.9 

2.3 

1.1 

ET 

JG (50) 

0.4 

0.8 

2.4 

1.5 

2.2 

0.4 

0.9 

1.5 

0.6 

1.2 

0,5 

E d { 0 + ) 

JG (40) 

0.2 

0.3 

0.2 

0.3 

0.3 

0.3 

0.2 

0.3 

0.3 

0.3 

0.1 


JS (40) 

0.6 

0.7 

0.8 

0.8 

0.9 

0.8 

0.7 

0.8 

0.8 

0.8 

0.2 


GS (40) 

0.6 

0.7 

0.6 

0.7 

0.8 

0.7 

0.6 

0.7 

0.6 

0.7 

0.2 

OT 

JG (50) 

1.9 

2.0 

2.2 

2.2 

2.3 

2.4 

2.4 

2.1 

2.4 

2.2 

0.1 

MT 

JG (50) 

0.2 

0.4 

0.2 

0.4 

0.3 

0.3 

0.2 

0.3 

0,3 

0.3 

0.1 

ET 

JG (50) 

0.3 

0.3 

0.3 

0.3 

0.3 

0.4 

0.3 

0.3 

0.4 

0.3 

0.0 
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JG (40) 
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1.3 

1.2 

1.3 

1.1 

3.2 

2.7 

1.3 

2.9 

1.7 

0.2 


JS (40) 

7.3 

7.0 

5.8 

5.7 

6.2 

23.0 

21.9 

6.4 

22.2 

10.6 

2.6 


GS (40) 

6.5 

6.4 

5.1 

4.9 

5.6 

22.5 

22.2 

5.7 

22.0 

10.1 

2.6 

OT 

JG (50) 

2.5 

2,6 

2.8 

2.8 

3.0 
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7.9 
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MT 
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1.0 

3.3 
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2.0 

1.4 

1.6 

1.7 

1.6 

0.1 
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4.9 
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3.8 
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0.4 
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9.2 
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7.0 

22.6 
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8.0 

21.9 
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4.2 
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6.5 
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5.1 

5.6 

22.4 
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8,3 

4.2 

OT 

JG (50) 
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4.3 

6.4 

4.6 

5.3 
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5.3 
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MT 
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4.0 
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4.0 
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4.0 

4.4 

0.4 

ET 
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3.9 
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0.3 
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0.7 
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0.6 
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0.1 

OT 
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0.4 

0.5 

0.6 

0.7 

4.7 
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0.5 
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MT 
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0.7 

0.7 

0.7 

0.7 

0.8 

3.1 
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0.7 

3.1 

1.4 

0.2 

ET 
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0.6 
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0.9 

0.6 

0.5 
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0.6 

0.1 
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MT 
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0.7 
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0.4 

0.3 

0.3 

0.3 

0.3 

0.4 

0.4 

0.3 

0.4 

0.3 

0.0 

All 

JG (40) 

1.7 

1.4 

1.4 

1.5 

1.3 

1.9 

1.7 

1.5 

1.9 

1.6 

0.4 
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R IS = 1.3(1 1.4)% : 4* 6 - 7 = 2.0(111.6)%, and 
$3/5 = 0.2(1 0.7)%; 

[L w ] n 4*!- 5 = 4.7(1 0.3)%. 4-6-7 = 3.8(113.0)%, and 

$3/5 = 0.4(4 o.6)% ; 

E u (0~) 4-i-s = 0.6(1 1.8)% s 4-6-7 = 2.3(111.6)%, and 
4-3/5 = 0.1(1 0.6)%; 

Q n { O') 4-i-s = 0.4(1 1.7)%, 4-6-7 = 1.0(1 3.9)%, and 
4-3/5 = 0.2(1 0.5)%; and 

All JG 4-i-s = 1.5(1 2.8)%. 4-6-7 = 1.9(1 8.7)%, and 
4-3/5 = 0.4(1 1.4)%. 

Note the Ed( 0 + ) uncertainties are the only ones that are 
unchanged. 

In terms of the individual parameters, the largest re- 
ductions are in the red wavelengths, except for the Kd 
values for which the largest change is in the blue-green do- 
main. The majority of the V2 uncertainties are less than 
2.5%, and almost all of the band-ratio uncertainties are 
less than 1.0%. The exceptions are Ed( 0“) at 683 nm, Kd 
for 412-555 nm, and [Lvv] n for 412-683 nm. This excel- 
lent level of agreement is well summarized by the overall 
JG averages (along the very bottom of Table 18) which 
are all less than 2.5% (note the overall band-ratio average 
is less than 1.0%). The excellent JG uncertainties are a 
result of using the same extrapolation intervals with the 
JRC and GSFC different processors. 

There is still a dependence in the uncertainties as a 
function of the concentration regimes. In general, the JG 
uncertainties decrease with increasing chlorophyll a con- 
centration. Exceptions to this are the R rs , [Lw/] n , blue- 
green Ed{0 + ), and blue-green L u ( 0") results for which the 
largest uncertainties are associated with the concentration 
regime end points (OT and ET). The E u { 0“) and the blue- 
green Q n (0') results show almost no dependence as a func- 
tion of concentration regime. 

1.6.2 V3 Summary 

The only changes made for the V3 processing with re- 
spect to the V2 processing was to remove the Ed( 0 + ) nor- 
malization usually used by the JRC processor, and to in- 
clude 5° tilt filtering for the GSFC processor. A summary 
of the effects of these changes is shown in Table 19. In 
comparison to the V2 results (Table 18), there are im- 
provements in most of the JG variables, but they are not as 
substantial a s the changes between the VI and V2 results — 
the Satlantic processor was not changed, so the JS and GS 
uncertainties are, again, largely unchanged. 

The most significant aspect of the improvements is the 
the majority of the JG uncertainties are at the 1.0% level 
or less (the | symbol represents an increase in uncertainty 
with respect to the V2 values): 

M0") * 1-5 = 0.3(| 0.2)%, *6-7 - 0.5(1 1.2)%, and 

* 3/5 = 0 . 0(1 0 . 1 )%; 

K l *!-5 - 0.2(1 0.6)%, * 6-7 = 0.3(1 1.7)%, and 
*3/5 = 0.0(1 0.4)%; 


E d { O') 

*1-5 

d 

II 

0.7)%, 

* 6-7 

= 1.4(4 

i.i)%, 

and 


4*3/5 

= 0.2(4 

0.1)%; 





Kd 

*1-5 

II 

O 

2.5)%, 

* 6-7 

d 

II 

0.4)%, 

and 


4*3/5 

d 

II 

1.2)%; 





E d { 0 + ) 

4*1-5 

cT 

d 

II 

0.1)%, 

* 6-7 

cT 

d 

II 

0.1)%, 

and 


4*3/5 

= 0.1(4 

0.0)%; 





Rrs 

4*1-5 

= 1-0(4 

0.3)%, 

* 6-7 

= 2.1(4 

0.1)%, and 


4*3/5 

cT 

d 

II 

0.0)%; 





M n 

4*1-5 

II 

0.1)%, 

* 6-7 

O 

CO 

II 

o 

oo 

and 

4*3/5 

d 

II 

0.0)%; 





E u (0-) 

*1-5 

II 

o 

to 

0.4)%, 

* 6-7 

II 

o 

cn 

oo 

and 


4*3/5 

II 

o 

d 

0.1)%; 





Qn((T) 

4*1-5 

= 0.2(4 

0.2)%, 

* 6-7 

II 

O 

CO 

0.7)%, 

and 


4*3/5 

= 0.0(4 

0.2)%; 

and 




All JG 

4*1-5 

II 

o 

0.6)%, 

* 6-7 

= 1.0(4 

0.9)%, 

and 


4*3/5 

= 0.1(4 

0.3)%. 






All of the JG band-ratio uncertainties are below 1.0% — 
in fact, the overall band- ratio average is 0.1%. Almost all 
of the spectral uncertainties are at the 1.0% level or less, 
except a) the red Ed(0~) and R rs values, b) some of the 
blue-green R rs values are at 1.1%, and c) all of the [Z>vv] N 
values are greater than 2.5%. All of the spectral overall 
averages are at the 1.0% level or less, except at 412 nm for 
which the overall average is 1.2%. 

There are still spectral dependencies in the results, but 
most of these are weak. For example, the difference be- 
tween the maximum and minimum L u ( 0“) uncertainty is 
0.3%. Exceptions to this are seen for Ed(0~) and R rs in 
the red domain, and [T\v] n in the blue-green. Overall, 
the difference between the maximum and minimum uncer- 
tainties for 412-683 nm is 0.5%, with most of the difference 
coming from the larger overall uncertainty at 412 nm. 

The convergence of the processors is well quantified by 
considering the average difference in the number of points 
in the extrapolation ranges, AiV, for the JRC and GSFC 
V3 processing of the SO, LN, WP, and MN data (the av- 
erage number of points in the extrapolation range is given 
in parentheses): 

SO AN = b (34), 

LN AN = 3 (61), 

WP AN = 2 (209), and 

MN AN = 0 (20). 

Although the AN values are small, the use of identical 
extrapolation ranges and very similar processing options 
does not result in exactly the same number of points in 
the extrapolation interval. The worst case is with Sea- 
OPS, but this is an instrument that is going up and down 
as a result of wave action (acting on the ship) as it is de- 
ployed. The surging motion can occur at an extrapolation 
boundary, and how the processors deal with the addition 
of new data after the boundary is first encountered can 
be different. The best case is with miniNESS data, where 
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Table 19. A summary of the t/j, 'F, and <f> values for the computation of the V3 primary variables. 


Comparison 

0(412) 

0(4 43) 

•0(490) 

0(510) 

•0(555) 

■0(665) 

0(683) 

Ms 

4 > 6-7 

$ 1-7 

$ 3/5 

MO') 

JG (40) 

0.4 

0.3 

0.2 

0.2 

0.2 

0.5 

0.4 

0.3 

0.5 

0.3 

0.0 


JS (40) 

6.8 

9.1 

5.0 

4.9 

5.3 

23.2 

21.9 

6.3 

22.3 

10.8 

2.4 


GS (40) 

6.6 

9.0 

4.9 

4.8 

5.2 

23.0 

21.8 

6.1 

22.2 

10.7 

2.4 

OT 

JG (50) 

0.1 

0.1 

0.1 

0.2 

0.3 

1.1 

1.3 

0.2 

1.2 

0.5 

; 0.1 

MT 

JG (50) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.2 

0.2 

0.1 

0.2 

0.1 

0.0 

ET 

JG (50) 

0.8 

0.7 

0.5 

0.5 

0.5 

1.1 

0.6 

0.6 

0.9 

0.7 

I o.i 

K b 

JG (40) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.3 

0.4 

0.2 

0.3 

0.2 

; 0.0 

OT 

JG (50) 

0.4 

0.5 

0.4 

0.4 

0.3 

0.8 

1.2 

0.4 

1.0 

0.6 

0.3 

MT 

JG (50) 

0.1 

0.1 

0.2 

0.2 

0.1 

0.3 

0.4 

0.1 

0.3 

0.2 

0.1 

ET 

JG (50) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.4 

0.3 

0.2 

0.4 

0.2 

0.0 

M O') 

JG (40) 

0.6 

0.6 

0.5 

0.5 

0.5 

1.2 

1.7 

0.5 

1.4 

0.8 

0.2 


JS (40) 

11.2 

11.2 

11.0 

10.9 

12.7 

32.4 

31.3 

11.4 

32.8 

16.9 

4.3 


GS (40) 

11.2 

11.2 

10.8 

10.8 

12.5 

32.9 

32.0 

11.3 

33.4 

17.0 

4.3 

OT 

JG (50) 

0.2 

0.2 

0.3 

0.3 

0.3 

0.6 

0.4 

0.3 

0.5 

0.3 

0.1 

MT 

JG (50) 

0.3 

0.3 

0.3 

0.3 

0.4 

0.7 

0.7 

0.3 

0.7 

0.4 

0.1 

ET 

JG (50) 

1.0 

1.0 

0.9 

0.7 

0.8 

2.0 

3.5 

0.9 

2.8 

1.4 

0.2 

K d 

JG (40) 

0.5 

0.6 

0.9 

0.8 

0.7 

0.4 

0.4 

0.7 

0.4 

0.6 

0.4 


JS (40) 

19.1 

20.3 

25.6 

26.5 

22.0 

14.1 

11.9 

23.8 

13.7 

20.9 

16.6 


GS (40) 

19.0 

20.3 

25.4 

26.3 

21.8 

14.0 

11.9 

23.6 

13.7 

20.7 

16.9 

OT 

JG (50) 

0.9 

1.1 

1.3 

0.9 

0.6 

0.2 

0.2 

1 1.0 

0.2 

0.7 

0.7 

MT 

JG (50) 

0.6 

0.7 

1.0 

1.2 

1.0 

0.4 

0.3 

0.9 

0.3 

0.7 

0.3 

ET 

JG (50) 

0.3 

0.4 

0.8 

0.4 

0.4 

0.4 

0.6 

0.5 

0.5 

0.5 

0.6 

M 0 + ) 

JG (40) 

0.1 

0.3 

0.2 

0.3 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.1 


JS (40) 

0.5 

0.6 

0.7 

0.7 

0.8 

0.7 

0.6 

0.7 

0.6 

0.6 

0.2 


GS (40) 

0.6 

0.7 

0.6 

0.7 

0.8 

0.7 

0.6 

0.7 

0.6 

0.7 

0.2 

OT 

JG (50) 

1.5 

1.6 

1.8 

1.8 

1.8 

1.9 

1.9 

1.7 

1.9 

1.8 

0.1 

MT 

JG (50) 

0.1 

0.3 

0.1 

0.3 

0.2 

0.3 

0.1 

0.2 

0.2 

0.2 

0.1 

ET 

JG (50) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.0 

•Rrs 

JG (40) 

1.1 

1.1 

1.0 

1.1 

1.0 

2.0 

2.2 

1.0 

2.1 

1.3 

0.2 


JS (40) 

7.2 

6.9 

5.6 

5.6 

6.1 

22.9 

21.9 

6.3 

22.1 

10.6 

2.6 


GS (40) 

6.6 

6.5 

5.1 

5.0 

5.6 

22.2 

21.8 

5.8 

21.6 

10.1 

2.6 

OT 

JG (50) 

1.8 

2.0 

2.1 

2.2 

2.2 

6.5 

5.0 

2.1 

5.7 

3.1 

0.4 

MT 

JG (50) 

| 0.8 

0.9 

0.8 

1.0 

0.8 

1.8 

2.1 

0.9 

2.0 

1.2 

0.2 

ET 

JG (50) 

1.5 

1.4 

1.2 

1.2 

1.2 

1.9 

1.2 

1.3 

1.5 

1.4 

0.1 

[M-] n 

JG (40) 

7.1 

4.2 

3.4 

4.8 

3.5 

3.0 


4.6 

3.0 

4.3 

0.4 

JS (40) 

1 10.3 

9.1 

7.9 

5.6 

6.8 

22.6 

21.5 

8.0 

21.8 

11.7 

4.1 


GS (40) 

6.2 

6.6 

5.6 

5.0 

5.6 

22.4 


5.8 

22.4 

8.3 

4.2 

OT 

JG (50) 

6.4 

4.0 

3.6 

5.6 

3.9 

3.6 


4.7 

3.6 

4.5 

0.3 

MT 

JG (50) 

6.8 

3.9 

3.4 

4.4 

3.3 

2.6 


4.4 

2.6 

4.1 

0.4 

ET 

JG (50) 

7.5 

4.6 

3.6 

5.3 

3.8 

3.6 


5.0 

3.6 

4.7 

0.3 

E u ( O') 

JG (40) 

0.3 

0.2 

0.2 

0.2 

0.2 

0.5 

0.4 

0.2 

0.5 

0.3 

0.0 

OT 

JG (50) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.2 

0.5 

0.0 

0.4 

0.1 

0.0 

MT 

JG (50) 

0.2 

0.2 

0.1 

0.1 

0.2 

0.4 

0.4 

0.2 

0.4 

0.2 

0.0 

ET 

JG (50) 

0.5 

0.4 

0.3 

0.3 

0.3 

0.8 

0.5 

0.4 

0.6 

0.4 

0.1 

Qn(O') 

JG (40) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.3 

0.3 

0.2 

0.3 

0.2 

0.0 


JG (50) 

0.1 

0.1 

0.1 

0.2 

0.2 

1.0 

1.4 

0.2 

1.2 

0.5 

0.1 


JG (50) 

0.2 

0.1 

0.1 

0.1 

0.2 

0.2 

0.3 

0.1 

0.2 

0.2 

0.0 


JG (50) 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.0 

All 

JG (40) 

1.2 

0.8 

0.8 

0.9 

0.7 

0.9 

0.7 

0.9 

1.0 

0.9 

0.1 
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the vertical resolution and extrapolation range are min- 
imized. Between these two, the processors are basically 
differing by one or two points at each extrapolation bound- 
ary. Differences at this level are to be expected, because 
the processors can be using different numerical precision 
in the calculation of depth, absolute or nearest neighbor 
limts can be employed for the boundaries, etc. 

Although there are variables for which the uncertainties 
either increase or decrease with increasing concentration, 
the majority do not show a significant dependence, i.e., 
the difference between the maximum and minimum uncer- 
tainty for the OT, MT, and ET results is usually less than 
1.0%. Notable exceptions are E d (0 + ), R rs , and [£w] N - 
The Ed{ 0 + ) result is associated with maximum uncertain- 
ties in the OT data which are composed solely of deep 
ocean data using unstabilized references. The R rs result is 
due to maximal uncertainties in the ET L u ( 0”) data and 
the OT Ed( 0 + ) data. The [Lw] N result is a consequence 
of the aforementioned differences in the Fq(\) coefficients. 

The DO, SC, WC, and FF partitioning of the V3 JG re- 
sults (using the 50-cast data set) are presented in Fig. 13. 
The dashed lines set the 2.5% intercomparison objective 
and the dotted lines the 1.0% objective (plotted within 
the 4% uncertainty range achieved with DARR-94, except 
for the normalized water-leaving radiance which has an ex- 
panded scale). In terms of the uncertainty levels, the vari- 
ables can be separated into three groups: the first group is 
composed of variables with uncertainties always less than 
1.0%, i.e., K l , E u ( 0"), and Q n (0 _ ); the second group has 
uncertainties frequently below 1%, with a few wavelengths 
slightly above 1%, i.e., L u ( 0“), Ed{ 0”), Kd, Ed{ 0 + ), and 
i?(0“); and the third group has uncertainties persistently 
above 1%, and even 2.5%, i.e., R rs and [Lvv] n - 

The uncertainties within the first group are so small 
and so similar, the differences within the various categories 
are probably not significant. There are, however, two no- 
table results associated with this group: a) the largest un- 
certainties are in the red wavelengths, and b) the Q n -factor 
is the only multivariate (i.e., normalized) parameter with 
uncertainties persistently below the 1.0% level. 

Within the second group, the differences between the 
various categories are more substantial: a) the shallow 
coastal data usually have lower uncertainties than the deep 
ocean data; b) the free-fall data usually have lower uncer- 
tainties than the winch and crane data; and c) the SC and 
FF L u ( 0"), E d ( 0"), E d { 0 + ), and R(0~) uncertainties are 
always less than 1.0% (even though the corresponding DO 
and WC values are not). 

The third group is important, because it consists of 
parameters central to calibration and validation activities, 
[Lw] n and £ rs , respectively. Unfortunately, both vari- 
ables are associated with the largest uncertainties. The 
normalized water-leaving radiances have large uncertain- 
ties, because of the differences in Fo(A), as discussed in 
Sect. 1.5.6. The R rs uncertainties are only large for the 
DO and WC categories— the SC uncertainties are always 


less than 1.0%, and the FF uncertainties are less than 1.0% 
in the blue-green domain (including the band-ratio value). 

There is also a general result supported by all three 
groups: the SC and FF uncertainties are less than the DO 
and WC values, respectively. The lower uncertainties with 
the SC data are a consequence of comparatively higher 
uncertainties for the SO data. The processing uncertainties 
for the JG (40-cast) data sets organized according to the 
four instrument types, are presented in Table 20 along with 
samples of the SO and LN uncertainties for the 50-cast 
data set. These samples show the differences between the 
two data sets for the SO and LN results are small and the 
general properties of the data (e.g., spectral dependence, 
overall averages, band ratios, etc.) are similar. 

For the L u ( 0"), £^(0"), R rs , E u ( 0 ), and Q n (0 ) var F 
ables, SO uncertainties are significantly larger than the 
LN, WP, and MN uncertainties. If the SO uncertainties 
are ignored, the latter three usually agree to within 0.5% 
(except for the red R rs wavelengths). 

The K l , E d { 0 + ), [Lw] n , and #( 0+ ) uncertainties in 
Table 20 do not show any significant dependence as a func- 
tion of the instrument type. Interestingly, the SO K L un- 
certainties are not distinguishable from the other instru- 
ment types even though the SO L u ( 0 ) uncertainties are. 
The low Ed( 0 + ) uncertainties suggest the difference be- 
tween stabilized and unstabilized references does not effect 
the processing uncertainties (when the processing options 
are made similar). The [Lw] n uncertainties do not show 
a strong sensitivity to instrument type, but the smallest 
values are usually associated with the LN data. In com- 
parison, the R( 0 + ) uncertainties also do not show a strong 
sensitivity to instrument type, but the largest values are 
usually associated with the LN data. 

As was seen with the Kl and L u ( 0 ) data, the Kd 
uncertainties do not follow the instrument dependence for 
E d ( 0“): the largest Kd uncertainties are associated with 
the WP data, whereas the largest £<*(0') uncertainties are 
for the SO data. Given the direct importance of the atten- 
uation coefficients in determining the extrapolated (null 
depth) values, it is interesting that the uncertainties for 
the two are not always commensurate. As with the VI 
and V2 results, the K d uncertainties in the red domain are 
significantly smaller than the blue- green uncertainties. 

The effect of processor convergence on the intercom- 
parison of the OC2v2 determination of total chlorophyll a 
between the JRC and GSFC data processors is as follows: 

SO 2.7%, 

LN 7.0%, 

WP 5.2%, and 
MN 7.2%, 

with an overall average of 5.5%. These are very similar to 
the VI results (Table 15), so even though the inter compar- 
ison of many of the primary optical variables has improved 
significantly, the net effect on the derived (OC2v2) chloro- 
phyll concentration is minimal. For comparison purposes, 
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□ H □ ■ 

DO SC WC FF 


0(412) 0(443) 0(490) 0(510) 0(555) 0(665) 0(683) 4>i- 5 \J> 6 - 7 0i- 7 $3/5 

Fig. 13. Average V3 processing uncertainties in the primary JG intercomparison variables for the DO (light 
gray bars). SC (dark gray bars). WC (white bars), and FF (black bars) data categories: a) L u (0"). b) Ki . 
c) Ed{ 0"), d) K d , and e) E d { 0 + ). 
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□ nan 

DO SC WC FF 


^(412) ^(443) -0(490) 0(510) ^(555) ^(665) 0(683) 4>i_ 5 4/ 6 - 7 * 1-7 $3/5 

Fig, 13 . (cont.) Average V3 processing uncertainties in the primary JG intercomparison variables for the 
DO (light gray bars), SC (dark gray bars), WC (white bars), and FF (black bars) data categories: f) R rs , g) 
[Lw( A)] n , h) £ u (0-), i) Q n { O'), and j) R( O'). 
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Table 20. A summary of the ?/;, 
40-cast data set, with samples from the 


and values for 
•50-cast data set, 


the V3 primary variables. The comparisons are for 
organized according to the four instrument types. 


the JG 


Comparison 

■0(412) 

■0(443) 

0(490) 

0(510) 

0(555) 

0(665) 

0(683) 

*1-5 

^6-7 

T 1-7 

$3/5 

L u { O') SO (40) 

1.3 

1.2 

0.9 

0.8 

0.8 

1.6 

1.0 

1.0 

1.3 

1.1 

0.1 

LN (40) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.4 

0.3 

0.1 

0.4 

0.2 

0.0 

WP (40) 

0.1 

0.0 

0.0 

0.0 

0.0 

0.1 

0.1 

0.0 

0.1 

0.1 

0.0 

MN (40) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

0.1 

0.0 

0.1 

0.0 

0.0 

SO (50) 

1.0 

0.9 

0.7 

0.6 

0.6 

1.3 

0.9 

0.8 

1.1 

0.9 

0.1 

LN (50) 

0.1 

0.1 

0.1 

0.2 

0.2 

1.0 

1.1 

0.2 

1.0 

0.4 

0.1 

K l SO (40) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.5 

0.4 

0.1 

0.4 

0.2 

0.1 

LN (40) 

0.1 

0.2 

0.2 

0.2 

0.2 

0.3 

0.5 

0.2 

0.4 

0.2 

0.1 

WP (40) 

0.3 

0.4 

0.3 

0.4 

0.4 

0.4 

0.4 

0.3 

0.4 

0.4 

0.0 

MN (40) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.0 

E d { 0") SO (40) 

1.4 

1.3 

1.1 

1.2 

1.2 

3.4 

5.3 

1.2 

4.4 

2.1 

0.1 

LN (40) 

0.4 

0.4 

0.5 

0.4 

0.4 

0.5 

0.7 

0.4 

0.6 

0.5 

0.3 

WP (40) 

0.4 

0.4 

0.4 

0.4 

0.5 

0.8 

0.7 

0.4 

0.7 

0.5 

0.2 

MN (40) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

0.1 

0.0 

0.1 

0.0 

0.0 

K d SO (40) 

0.2 

0.2 

0.4 

0.4 

0.4 

0.3 

0.7 

0.3 

0.5 

0.4 

0.1 

LN (40) 

0.7 

0.8 

1.6 

0.8 

0.5 

0.1 

0.2 

0.8 

0.2 

0.6 

1.0 

WP (40) 

1.1 

1.3 

1.6 

2.0 

2.0 

0.9 

0.8 

1.6 

0.8 

1.4 

0.4 

MN (40) 

0.1 

0.0 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.0 

SO (50) 

0.2 

0.3 

0.4 

0.4 

0.4 

0.4 

0.6 

0.3 

0.5 

0.4 

0.1 

LN (50) 

0.9 

1.1 

1.7 

1.0 

0.7 

0.2 

0.2 

1.0 

0.2 

0.8 

1.1 

E d { 0 + ) SO (40) 

0.1 

0.1 

0.2 

0.2 

0.2 

0.2 

0.2 

0.1 

0.2 

0.2 

0.0 

LN (40) 

0.2 

0.7 

0.2 

0.6 

0.4 

0.6 

0.2 

0.4 

0.4 

0.4 

0.2 

WP (40) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.0 

MN (40) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

R ts SO (40) 

2.1 

1.9 

1.6 

1.5 

1.4 

2.8 

1.9 

1.7 

2.3 

1.9 

0.2 

LN (40) 

0.7 

1.1 

0.9 

1.2 

0.9 

3.9 

5.2 

1.0 

4.6 

2.0 

0.4 

WP (40) 

0.8 

0.7 

0.8 

0.7 

0.8 

0.7 

0.8 

0.8 

0.7 

0.7 

0.1 

MN (40) 

0.7 

0.7 

0.7 

0.7 

0.8 

0.8 

0.8 

0.7 

0.8 

0.8 

0.1 

SO (50) 

2.8 

2.8 

2.7 

2.7 

2.7 

5.3 

2.7 

2.8 

4.0 

3.1 

0.3 

LN (50) 

0.7 

0.9 

0.8 

1.1 

0.8 

4.2 

4.6 

0.9 

4.4 

1.9 

0.3 

£ 

CO 

O 

8.0 

4.8 

3.8 

6.5 

4.1 

4.0 

0.0 

5.5 

4.0 

5.2 

0.3 

LN (40) 

6.5 

3.4 

3.0 

5.8 

3.2 

2.5 

0.0 

4.4 

2.5 

4.1 

0.6 

WP (40) 

7.0 

4.3 

3.3 

3.8 

3.4 

2.6 

0.0 

4.4 

2.6 

4.1 

0.3 

MN (40) 

6.7 

4.2 

3.6 

3.3 

3.4 

2.7 

0.0 

4.2 

2.7 

4.0 

0.4 

SO (50) 

7.4 

4.9 

4.3 

6.3 

4.6 

4.6 

0.0 

5.5 

4.6 

5.3 

0.3 

LN (50) 

6.6 

3.5 

3.0 

5.8 

3.2 

2.6 

0.0 

4.4 

2.6 

4.1 

0.4 

E u { O') SO (40) 

0.8 

0.7 

0.5 

0.4 

0.4 

1.3 

0.7 

0.5 

1.0 

0.7 

0.1 

LN (40) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

0.4 

0.0 

0.4 

0.1 

0.0 

WP (40) 

0.1 

0.1 

0.0 

0.0 

0.0 

0.1 

0.1 

0.1 

0.1 

0.1 

0.0 

MN (40) 

0.3 

0.3 

0.2 

0.3 

0.3 

0.5 

0.5 

0.3 

0.5 

0.3 

0.0 

Qn( O') so (40) 

0.5 

0.5 

0.4 

0.4 

0.4 

0.4 

0.4 

0.5 

0.4 

0.4 

0.0 

LN (40) 

0.0 

0.0 

0.0 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.0 

WP (40) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

0.1 

0.0 

0.1 

0.1 

0.0 

MN (40) 

0.2 

0.2 

0.2 

0.2 

0.3 

0.5 

0.5 

0.2 

0.5 

0.3 

0.0 

i?(0‘) SO (40) 

0.1 

0.1 

0.2 

0.2 

0.2 

0.2 

0.2 

0.1 

0.2 

0.2 

0.0 

LN (40) 

0.2 

0.7 

0.2 

0.6 

0.4 

0.6 

0.2 

0.4 

0.4 

0.4 

0.2 

WP (40) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.0 

MN (40) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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the JS and GS average uncertainties for the OC2v2 in- 
tercomparison are 12.2 and 14.2%, respectively, which are 
very similar to the VI results. 

1.6.3 V4 Summary 

The only change made for the V4 processing with re- 
spect to the V3 processing was to use a robust linear fit for 
the GSFC processor. The robust linear fit weights the data 
in the extrapolation interval to minimize the significance of 
data points with higher variance, which acts as a mild form 
of outlier rejection. The effects of this change are presented 
in Table 21. In comparison to the very low uncertainties 
achieved with the V3 results (Table 19), there are com- 
paratively large increases in the uncertainties of the JG 
variables — the Satlantic processor was not changed, so the 
JS and GS uncertainties are, again, largely unchanged. 

The most significant aspect of the change in uncertain- 
ties for the V4 processing is the majority of the JG uncer- 
tainties are no longer at the 1.0% level or less, although 
many are to within 2.5%: 


L u ( O') 

& 

Cn 

II 

O 

bo 

0.5)%, 

$6-7 

II 

to 

to 

1.7)%, 

and 


$3/5 — 0.5(T 

0.5)%; 





K l 

$1-5 = 1.5(T 

1.3)%, 

$6-7 

= 2.2(1 

1.9)%, 

and 


$3/5 = 1 • 1 (T 

1.1)%; 





E d { O') 

$1-5 = 2.8(T 

2.3)%, 

$6-7 

= 3.8(1 

2.4)%, 

and 


$3/5 = 2.2(T 

2.0)%; 





K d 

$1-5 = 6.4(1 

5.7)%, 

$6-7 

= 2.0(1 

1.6)%, 

and 


$3/5 = 8.9(1 

8.5)%; 





E d ( 0 + ) 

$1-5 = 0.2(T 

0.0)%, 

$6-7 

= 0.2(1 

0.0)%, 

and 


$3/5 = 0.1(1 

0.0)%; 





-Rrs 

$1-5 = 1.2(1 

0.2)%, 

$6-7 

= 3.0(1 

0.9)%, 

and 


$3/5 = 0.5(1 

0.3)%; 





[ l w]n 

$1-5 = 4.6(1 

0.0)%, 

$6-7 

= 4.1(1 

1.1)%, 

and 

$3/5 — 0.6(1 

0.2)%; 





£ 

o 

1 

$1-5 = 0.7(1 

0.5)%, 

$6-7 

= 2.5(1 

2.0)%, 

and 


$3/5 — 0.4(1 

0.4)%; 





Qn(O') 

$1-5 = 0.9(T 

0.7)%, 

$6-7 

= 2.3(1 

2.0)%, 

and 


$3/5 = 0.6(1 

0.6)%; 

and 




All JG 

$1-5 = 2.1(1 

1.2)%, 

$6-7 

= 2.5(1 

1.5)%, 

and 


$3/5 - 1.7(1 

1.6)%. 






The R rs and [Lw] n uncertainties are not very differ- 
ent in terms of the percent increase of the V4 results with 
respect to the V3 results. The E d ( 0 + ) uncertainties are 
unchanged, because there is nothing in the robust linear 
fitting that can influence the calculation of E d ( 0 + ). The 
R rs and [L w ] n are not strongly altered (from a percent 
change point of view), because the normalizing variable, 
E d (0 + ), is unchanged, and the absolute increases in the 
L u { 0") uncertainties are proportionately smaller, particu- 
larly for [Lw] n - 

Although the JG L u { 0") and K L uncertainties show 
similar increases, the increases in the E d (0 ) and K d un- 
certainties are considerably different. This shows the type 


of fit is sensitive to the amount of noise in the data — the 
E d (0~) data are noisier than the L u (0") data, because of 
wave focusing effects. This can also be seen in the con- 
centration regime data. The L u ( 0 ) uncertainties do not 
show a clear trend as a function of concentration, but the 
E d { 0") uncertainties decrease as concentration increases. 

In most cases, the band-ratio uncertainties are much 
less than the spectral uncertainties, and most are less than 
1.0%. The notable exception is the K d results for which 
the band-ratio uncertainty is very large and responsible for 
increasing the overall band-ratio uncertainty above 1.0% 
(but still within 2.5%). The latter is caused by the signif- 
icantly larger K d uncertainties in the blue-green domain. 

The sensitivity of the primary variables to chlorophyll a 
concentration is a mixture of increasing and decreasing un- 
certainties as a function of increasing concentration. Per- 
haps the most significant results in this regard are the very 
large K d uncertainties in the OT and even the MT regimes. 
This points to the difficulty of estimating the slope of the 
log-transformed light data in clear waters where surface 
effects add a considerable amount of noise to the data. 


1.6.4 V5 Summary 


The primary purpose of the V5 processing configura- 
tion was to investigate the importance of filtering on the 
uncertainties, so the GSFC processor was set to its usual 
configuration (no tilt filtering and no outlier filtering), and 
the JRC processor was changed slightly: tilt filtering was 
enabled (as usual), but outlier filtering was disabled (usu- 
ally outlier filtering was enabled). A summary of the pro- 
cessing options is given in Table 17, and the effects of the 
V5 changes is shown in Table 22. In comparison to the V3 
results (Table 19), there are increases in uncertainties for 
most of the JG variables, but they are not as substantial 
as the changes seen with the V4 results (Table 21) — the 
Satlantic processor was not modified, so the JS and GS 
uncertainties remain largely unchanged. 

The most significant aspect of the change in uncertain- 
ties for the V5 processing with respect to the V3 processing 
is there is an overall increase in uncertainty values, but the 
majority of the JG uncertainties for the blue-green wave- 
lengths are within, or close to, the 1.0% level, and most of 
the uncertainties in the red domain are less than 2.5%: 


L u { 0") $i- 5 = 0.5(T 0.2)%, 4/ 6- 7 - 1.7(T 1-2)%, and 
$3/5 = 0.1(T 0.1)%; 

Kl $ 1-5 = 0.7(T 0.5)%, $ 6-7 = 1.9(1 1-6)%, and 
$3/5 = 0.3(T 0.3)%; 

E d { 0 ") $ 1-5 - 1-2(T 0.7)%, $ 6 -7 - 2.4(T 1-0)%, and 
$3/5 - 0.3(1 0.1)%; 

K d $ 1-5 = 3.0(T 2.3)%, $ 6-7 = 0.8(T 0.4)%, and 
$ 3/5 = i-5(T 1.1)%; 

E d { 0 + ) $ 1-5 = 0.2(T 0.0)%, $ 6-7 = 0.2(T 0.0)%, and 
$3/5 - 0.1(T 0.0)%; 
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Table 21. A summary of the ip, and $ values for the computation of the V4 primary variables. 


Comparison 

1^(412) 

ip( 443) 

0(490) 

ip( 510) 

-0(555) 

0(665) 

-0(683) 

Ms 

^ 6-7 

T1-7 

$ 3/5 

MO') 

JG (40) 

0.9 

0.8 

0.8 

0.8 

0.8 

2.2 

2.2 

0.8 

2.2 

1.2 

0.5 


JS (40) 

6.8 

9.1 

5.0 

4.9 

5.3 

23.2 

21.9 

6.3 

22.3 

10.8 

2.4 


GS (40) 

6.7 

9.0 

4.9 

4.8 

5.3 

23.4 

22.1 

6.2 

22.5 

10.8 

2.3 

OT 

JG (50) 

0.4 

0.6 

0.6 

0.7 

1.2 

6.2 

7.1 

0.7 

6.7 

2.4 

1.2 

MT 

JG (50) 

0.5 

0.5 

0.6 

0.6 

0.7 

2.0 

1.9 

0.6 

1.9 

1.0 

0.4 

ET 

JG (50) 

1.5 

1.3 

1.0 

1.0 

0.9 

2.3 

2.0 

1.1 

2.1 

1.4 

0.5 

K l 

JG (40) 

0.9 

1.2 

1.7 

1.6 

1.9 

1.7 

2.6 

1.5 

2.2 

1.7 

1.1 

OT 

JG (50) 

1.1 

2.0 

2.4 

2.0 

3.0 

6.8 

9.0 

2.1 

7.9 

3.8 

3.3 

MT 

JG (50) 

1.1 

1.2 

2.0 

2.1 

2.2 

1.8 

2.5 

1.7 

2.1 

1.8 

0.8 

ET 

JG (50) 

0.7 

0.9 

0.9 

0.9 

1.0 

1.2 

1.4 

0.9 

1.3 

1.0 

0.6 

MO') 

JG (40) 

2.6 

2.8 

2.7 

3.0 

2.7 

3.3 

4.2 

2.8 

3.8 

3.0 

2.2 


JS (40) 

11.2 

11.2 

11.0 

10.9 

12.7 

32.4 

31.3 

11.4 

32.8 

16.9 

4.3 


GS (40) 

11.1 

11.1 

10.7 

11.0 

12.4 

33.8 

32.2 

11.3 

33.8 

17.0 

5.2 

OT 

JG (50) 

4.0 

4.2 

4.7 

5.0 

5.1 

4.0 

5.3 

4.6 

4.6 

4.6 

2.7 

MT 

JG (50) 

2.8 

2.6 

2.8 

3.5 

3.4 

3.8 

4.1 

3.0 

3.9 

3.3 

2.3 

ET 

JG (50) 

2.0 

2.7 

2.0 

1.7 

1.5 

2.7 

4.7 

2.0 

3.7 

2.5 

1.6 

Kd 

JG (40) 

4.2 

5.4 

10.1 

6.7 

5.7 

1.7 

2.3 

6.4 

2.0 

5.1 

8.9 


JS (40) 

19.1 

20.3 

25.6 

26.5 

22.0 

14.1 

11.9 

23.8 

13.7 

20.9 

16.6 


GS (40) 

19.5 

18.9 

28.4 

26.3 

23.5 

14.6 

12.3 

24.2 

14.2 

21.3 

21.2 

OT 

JG (50) 

19.5 

20.1 

19.4 

10.6 

7.4 

1.3 

1.1 

15.4 

1.2 

11.3 

13.2 

MT 

JG (50) 

5.8 

7.2 

12.8 

9.0 

8.0 

2.4 

3.2 

8.6 

2.8 

6.8 

11.8 

ET 

JG (50) 

1.1 

2.0 

3.2 

2.6 

2.4 

0.9 

1.3 

2.3 

1.1 

1.9 

2.4 

E d (0 + ) 

JG (40) 

0.1 

0.3 

0.2 

0.3 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.1 


JS (40) 

0.5 

0.6 

0.7 

0.7 

0.8 

0.7 

0.6 

0.7 

0.6 

0.6 

0.2 


GS (40) 

0.6 

0.7 

0.6 

0.7 

0.8 

0.7 

0.6 

0.7 

0.6 

0.7 

0.2 

OT 

JG (50) 

1.5 

1.6 

1.8 

1.8 

1.8 

1.9 

1.9 

1.7 

1.9 

1.8 

0.1 

MT 

JG (50) 

0.1 

0.3 

0.1 

0.3 

0.2 

0.3 

0.1 

0.2 

0.2 

0.2 

0.1 

ET 

JG (50) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.0 

Rrs 

JG (40) 

1.2 

1.2 

1.0 

1.3 

1.2 

3.0 

3.1 

1.2 

3.0 

1.7 

0.5 


JS (40) 

7.2 

6.9 

5.6 

5.6 

6.1 

22.9 

21.9 

6.3 

22.1 

10.6 

2.6 


GS (40) 

6.7 

6.5 

5.2 

5.0 

5.6 

23.1 

22.4 

5.8 

22.4 

10.3 

2.6 

OT 

JG (50) 

1.7 

1.9 

2.2 

2.4 

3.0 

6.5 

7.7 

2.2 

7.1 

3.6 

1.2 

MT 

JG (50) 

1.0 

1.1 

0.9 

1.2 

1.1 

2.9 

3.3 

1.1 

3.1 

1.6 

0.3 

ET 

JG (50) 

■ 1.7 

1.5 

1.3 

1.5 

1.3 

2.9 

1.9 

1.5 

2.4 

1.7 

0.5 

[ m] n 

JG (40) 

7.0 

4.0 

3.4 

4.9 

3.7 

4.1 


4.6 

4.1 

4.5 

0.6 

JS (40) 

10.3 

9.1 

7.9 

5.6 

6.8 

22.6 

21.5 

8.0 

21.8 

11.7 

4.1 


GS (40) 

6.3 

6.7 

5.6 

5.0 

5.6 

22.7 


5.8 

22.7 

8.4 

4.1 

OT 

JG (50) 

6.2 

3.7 

3.3 

5.4 

4.6 

8.0 


4.6 

8.0 

5.2 

1.5 

MT 

JG (50) 

6.9 

4.0 

3.4 

4.5 

3.6 

3.7 


4.5 

3.7 

4.4 

0.5 

ET 

JG (50) 

7.2 

4.2 

3.7 

5.5 

3.6 

4.2 


4.8 

4.2 

4.7 

0.5 

E u(O') 

JG (40) 

0.7 

0.7 

0.7 

0.6 

0.7 

3.1 

2.0 

0.7 

2.5 

1.2 

0.4 

OT 

JG (50) 

0.4 

0.4 

0.3 

0.4 

0.4 

8.2 

7.0 

0.4 

7.6 

2.4 

0.4 

MT 

JG (50) 

0.7 

0.7 

0.8 

0.7 

0.8 

2.9 

2.1 

0.8 

2.5 

1.3 

0.4 

ET 

JG (50) 

0.6 

0.7 

0.6 

0.5 

0.6 

2.5 

1.2 

0.6 

1.9 

1.0 

0.4 

Qn(O') 

JG (40) 

0.9 

0.9 

0.9 

0.9 

0.9 

2.1 

2.4 

0.9 

2.3 

1.3 

0.6 

OT 

JG (50) 

0.6 

0.5 

0.7 

0.7 

1.2 

3.4 

4.2 

0.7 

3.8 

1.6 

1.3 

MT 

JG (50) 

0.6 

0.8 

1.0 

0.9 

1.0 

2.4 

2.9 

0.8 

2.6 

1.3 

0.5 

ET 

JG (50) 

1.2 

1.2 

0.9 

1.0 

0.8 

1.5 

1.6 

1.0 

1.5 

1.2 

0.6 

All 

JG (40) 

2.1 

1.9 

2.4 

2.2 

2.0 

2.4 

2.4 

2.1 

2.5 

2.2 

1.7 
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Table 22. A summary of the 1 /;, 'L, and values for the computation of the V5 primary variables. 


Comparison 

*0(412) 

-0(443) 

0(490) 

0(510) 

0(555) 

665) 

V>(683) 

*1-5 

$6-7 

*1-7 

$3/5 

L u { O') 

JG (40) 

0.6 

0.6 

0.4 

0.4 

0.5 

1.8 

1.7 

0.5 

1.7 

0.8 

0.1 


JS (40) 

6.8 

9.1 

5.0 

4.9 

5.3 

23.1 

21.8 

6.2 

22.3 

10.8 

2.4 


GS (40) 

6.5 

8.9 

4.8 

4.7 

5.2 

23.0 

21.7 

6.0 

21.9 

10.5 

2.4 

OT 

JG (50) 

0.4 

0.4 

0.4 

0.5 

0.6 

4.0 

5.0 

0.4 

4.5 

1.6 

0.4 

MT 

JG (50) 

0.5 

0.5 

0.4 

0.4 

0.4 

2.1 

1.9 

0.4 

2.0 

0.9 

0.1 

ET 

JG (50) 

0.9 

0.8 

0.6 

0.6 

0.5 

1.0 

0.7 

0.7 

0.8 

0.7 

0.1 

K l 

JG (40) 

0.7 

0.7 

0.7 

0.6 

0.6 

1.6 

2.1 

0.7 

1.9 

1.0 

0.3 

OT 

JG (50) 

1.2 

1.5 

1.5 

1.7 

1.3 

4.9 

6.3 

1.4 

5.6 

2.6 

0.8 

MT 

JG (50) 

1.0 

1.0 

1.0 

0.9 

0.8 

2.1 

2.7 

0.9 

2.4 

1.3 

0.5 

ET 

JG (50) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.4 

0.3 

0.2 

0.3 

0.2 

0.0 

E d ( O') 

JG (40) 

1.2 

1.2 

1.2 

1.1 

1.4 

2.1 

2.8 

1.2 

2.4 

1.6 

0.3 


JS (40) 

11.3 

11.3 

11.0 

11.0 

12.7 

32.3 

31.3 

11.4 

32.8 

16.9 

4.3 


GS (40) 

10.7 

10.6 

10.1 

10.2 

11.8 

32.9 

31.6 

10.7 

33.2 

16.5 

4.3 

OT 

JG (50) 

2.4 

2.5 

2.5 

2.6 

2.7 

3.4 

3.0 

2.5 

3.2 

2.7 

0.3 

MT 

JG (50) 

1.1 

1.1 

1.0 

1.0 

1.3 

2.1 

2.0 

1.1 

2.0 

1.4 

0.3 

ET 

JG (50) 

0.8 

0.8 

1.0 

0.8 

1.0 

1.6 

3.8 

0.9 

2.7 

1.4 

0.1 

K d 

JG (40) 

2.8 

3.3 

4.0 

3.0 

2.6 

0.8 

0.8 

3.0 

0.8 

2.4 

1.5 


JS (40) 

19.1 

20.4 

25.7 

26.6 

22.1 

14.1 

11.9 

23.8 

13.7 

20.9 

16.6 


GS (40) 

18.9 

19.4 

23.8 

25.3 

20.7 

13.9 

11.6 

22.6 

13.5 

20.1 

16.8 

OT 

JG (50) 

8.3 

9.5 

9.3 

5.8 

3.5 

0.9 

0.6 

7.3 

0.7 

5.4 

5.9 

MT 

JG (50) 

2.3 

3.0 

3.3 

3.0 

2.7 

1.0 

0.8 

2.9 

0.9 

2.3 

1.1 

ET 

JG (50) 

0.2 

0.5 

2.1 

1.3 

1.8 

0.3 

0.8 

1.2 

0.5 

1.0 

0.3 

E d ( 0 + ) 

JG (40) 

0.1 

0.2 

0.1 

0.2 

0.2 

0.2 

0.1 

0.2 

0.2 

0.2 

0.1 


JS (40) 

0.5 

0.6 

0.7 

0.7 

0.8 

0.7 

0.6 

0.7 

0.6 

0.7 

0.2 


GS (40) 

0.6 

0.7 

0.6 

0.7 

0.8 

0.7 

0.6 

0.7 

0.6 

0.7 

0.2 

OT 

JG (50) 

1.5 

1.6 

1.7 

1.8 

1.8 

1.9 

1.9 

1.7 

1.9 

1.8 

0.1 

MT 

JG (50) 

0.1 

0.3 

0.1 

0.3 

0.2 

0.3 

0.1 

0.2 

0.2 

0.2 

0.1 

ET 

JG (50) 

0.1 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.0 

R rs 

JG (40) 

1.3 

1.2 

1.1 

1.2 

1.0 

3.1 

2.6 

1.2 

2.8 

1.6 

0.2 


JS (40) 

7.2 

6.8 

5.6 

5.6 

6.1 

22.9 

21.8 

6.3 

22.1 

10.5 

2.6 


GS (40) 

6.5 

6.4 

5.1 

4.9 

5.6 

22.5 

22.2 

5.7 

22.0 

10.1 

2.6 

OT 

JG (50) 

2.1 

2.3 

2.4 

2.4 

2.5 

8.1 

6.7 

2.3 

7.4 

3.8 

0.7 

MT 

JG (50) 

1.2 

1.1 

1.0 

1.2 

1.0 

3.2 

2.9 

1.1 

3.1 

1.6 

0.3 

ET 

JG (50) 

1.6 

1.5 

1.3 

1.3 

1.2 

1.9 

1.2 

1.4 

1.5 

1.4 

0.1 

[Lw} n 

JG (40) 

7.2 

4.2 

3.5 

4.9 

3.6 

3.8 


4.7 

3.8 

4.5 

0.4 

JS (40) 

10.3 

9.1 

7.9 

5.6 

6.8 

22.6 

21.5 

8.0 

21.8 

11.7 

4.1 


GS (40) 

6.1 

6.5 

5.4 

5.1 

5.6 

22.4 


5.7 

22.4 

8.3 

4.2 

OT 

JG (50) 

6.7 

4.3 

3.9 

5.9 

4.2 

4.8 


5.0 

4.8 

5.0 

0.5 

MT 

JG (50) 

7.0 

4.0 

3.5 

4.5 

3.4 

3.9 


4.5 

3.9 

4.4 

0.4 

ET 

JG (50) 

7.6 

4.7 

3.7 

5.4 

3.9 

3.6 


5.0 

3.6 

4.8 

0.3 

E u ( O') 

JG (40) 

0.7 

0.6 

0.5 

0.5 

0.6 

2.4 

2.1 

0.6 

2.3 

1.1 

0.1 

OT 

JG (50) 

0.2 

0.2 

0.2 

0.3 

0.4 

4.5 

4.8 

0.3 

4.6 

1.5 

0.2 

MT 

JG (50) 

0.7 

0.7 

0.6 

0.7 

0.8 

3.1 

3.0 

0.7 

3.1 

1.4 

0.2 

ET 

JG (50) 

0.6 

0.6 

0.4 

0.4 

0.3 

0.8 

0.5 

0.5 

0.7 

0.5 

0.1 

Qn( on 

JG (40) 

0.4 

0.4 

0.4 

0.5 

0.5 

1.0 

1.0 

0.4 

1.0 

0.6 

0.2 

OT 

JG (50) 

0.3 

0.3 

0.3 

0.4 

0.4 

2.4 

2.8 

0.3 

2.6 

1.0 

0.3 

MT 

JG (50) 

0.4 

0.5 

0.5 

0.6 

0.7 

1.1 

1.2 

0.5 

1.2 

0.7 

0.2 

ET 

JG (50) 

0.4 

0.3 

0.3 

0.3 

0.3 

0.4 

0.3 

0.3 

0.4 

0.3 

0.0 

All 

JG (40) 

1.7 

1.4 

1.3 

1.4 

1.2 

1.9 

1.7 

1.4 

1.9 

1.5 

0.4 
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Rrs 

*1-5 

= 1.2(t 

0.2)%, 

*6-7 

= 2.8(T 

0.7)%, 

and 


$3/5 

= 0.2(T 

0.0)%; 





[£w] n 

*1-5 

- 4- 7( T 

0.1)%, 

*6-7 

= 3.8(T 

0.8)%, 

and 

$3/5 

= 0.4(1 

0.0)%; 





E u (0~) 

*1-5 

= 0.6(T 

0.4)%, 

*6-7 

= 2.3(1 

1.8)%, 

and 


*1*3/5 

= 0.1(1 

0.1)%; 





Q n(0-) 

*1-5 

= 0.4(T 

0.2)%, 

*6-7 

= 1.0(1 

0.7)%, 

and 


$3/5 

= 0.2(T 

0.2)%; 

and 




All JG 

$1-5 

= l.4(T 

0.5)%, 

*6-7 

= 1.9(T 

0.9)%, 

and 


$3/5 

= 0.4(T 

0.3)%. 





Note that in 

many cases, the 

increases in uncertainties 


due to filtering differences is approximately one half (or 
more) of the uncertainty levels achieved when the proces- 
sors were the most similar (V3). 

As was seen with the V4 and V3 comparisons, the V5 
Ed( 0 + ), R rs , and [Lw] n uncertainties are very similar to 
the V3 results, particularly in terms of the percent in- 
creases of the V5 results with respect to the V3 results. 
Again, the Ed{ 0 + ) uncertainties are actually unchanged, 
and the R rs and [Lw] n are not strongly altered (from 
a percent change point of view), because the normalizing 
variable, Ed( 0 + ) 5 is unchanged, and the absolute increases 
in the L u ( 0”) uncertainties are proportionately smaller, 
particularly for [Lw'] n - 

The Ki uncertainties are similar to the L u ( 0") un- 
certainties, both in magnitude and spectral dependence. 
The Kd uncertainties are larger than the Ed( 0") uncer- 
tainty values, and have an inverted spectral dependence, 

i.e., when one is maximal, the other is minimal. With the 
exception of the Kd results, the band-ratio uncertainties 
are much less than the spectral uncertainties, and always 
less than 1.0%. 

1.7 DISCUSSION 

The most significant changes in the application of the 
various JRC and GSFC processing options occur from the 
progression from the VI, V2, and then V3 options. A 
summary of the evolution in uncertainties for the primary 
intercomparison variables is shown in Fig. 14 (the dashed 
and dotted lines set the 2.5% and 1.0% intercomparison 
objectives, respectively). The figure shows a significant 
reduction in processing uncertainties as a function of the 
processing versions (VI to V3) except for the Ed(0 + ) and 
[Lw(A)] n variables. The former are small to begin with 
and stay small, and the latter are large to begin with and 
stay large. Another exception is some of the band-ratio 
results. The L u ( O’), R rs , [Lw(A)] n , and E u ( 0") band- 
ratio uncertainties are less than 1.0% for the VI processing 
and remain so up through the V3 processing. 

Other aspects to note from Fig. 14 are as follows: 

1. The most significant improvement is seen in the red 
wavelengths, except for Kd which has the largest 
improvement in the blue-green domain; 


2. All of the uncertainties are reduced to within 1.0% 
for the V3 processing except for the Ed{ 0”) red 
wavelengths, all the [L\y(A)] N values, and the R vs 
values in the red domain; 

3. In terms of overall averages, the 1.0% objective is 
achieved at all wavelengths except at 412 nm which 
is slightly above 1.0%; and 

4. The largest improvement occurs for the VI to V2 
progression, i.e., making the extrapolation ranges 
the same has the strongest influence in reducing 
processing uncertainties. 

The uncertainties of the three processors based on the 
V3 processing of four of the most important primary vari- 
ables are presented in Fig. 15. The L u (0”) results (Fig. 15a) 
show the lowest uncertainties with the Satlantic proces- 
sor are for the shallow coastal data, and in all cases, the 
largest uncertainties are for the red wavelengths — none of 
the Satlantic uncertainties are below 1.0%, although, the 
blue-green WiSPER wavelengths are within 2.5%. The av- 
erage band-ratio uncertainty, $3/5, for both Satlantic in- 
tercomparisons are, however, below 2.5% (Table 19). The 
JG uncertainties are always below 2.5%, and mostly below 
1.0% except for the blue and red SeaOPS data. The Kd un- 
certainties (Fig. 15b) show much higher overall uncertain- 
ties with respect to the Satlantic processor, particularly 
for the LoCNESS data set where every wavelength has an 
uncertainty above 25%. None of the Satlantic uncertain- 
ties are below 2.5%, whereas all of the JG uncertainties 
are, and most of the JG uncertainties are to within 1.0% 
(the primary exception are the WiSPER data). 

The R rs uncertainties (Fig. 15c) are similar to the 
L u ( 0~) results, especially for the Satlantic intercompar- 
isons. The lowest uncertainties with the Satlantic proces- 
sor are for the shallow coastal data, and in all cases, the 
largest uncertainties are for the red wavelengths — none of 
the Satlantic uncertainties are below 1.0%, although, the 
blue-green WiSPER wavelengths are within 2.5%. Again, 
the overall band-ratio results are much improved and are 
a little larger than 2.5% (Table 19). Almost all of the JG 
uncertainties are below 2.5% (except for red deep ocean 
wavelengths), and the majority of the uncertainties are be- 
low 1.0% (all of the shallow coastal wavelengths are within 
1.0%). 

In terms of the identified processing configurations, the 
JRC and GSFC processors were the most similar for the V3 
options, and the minimum uncertainties, both individually 
and from an overall perspective, were seen with the V3 
results. 

The V4 processing results show the type of fit used to 
estimate the slope of the log-transformed data in the ex- 
trapolation interval is correlated with the amount of noise 
in the data, thus, the noisier Ed(0~) data (from wave fo- 
cusing effects) are more sensitive to the extrapolation pro- 
cess. This can also be seen in the concentration regime 
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Fig. 14. (cont.) Average uncertainties in the primary intercomparison variables for the VI (light gray bars). 
V2 (dark gray bars), and V3 (black bars) .JG processing: f) R rs: g) [Lw(A)] n , h) £J U (0"). i) Q„(0"). and j) 
all the variables. 
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Instrument and Intercomparison 



Fig. 15 . The average V3 processing uncertainties for the JG, JS, and GS intercomparisons of a) L u { 0”) 
and b) Kd . The dashed line sets the 2.5% intercomparison objective and the dotted line the 1.0% objective. 
Uncertainties above 25% are shown clipped at the 25% level, and the maximum uncertainty achieved is given 
at the top of the panel to the side of the clipped bars. The individual wavelengths are shown as the sequential 
bars with varying intensities of gray (going from left to right, blue is light gray and red is dark gray). The 
instrument (and deployment location) codes are given along the top of the subpanels. The overall average of 
all the data (40 casts) is given in the right-most three intercomparisons. 


43 






Results of the Second SeaWiFS Data Anal}- sis Round Robin. March 2000 (DARR-00) 
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JG JS GS JG JS GS JG JS GS JG JS GS JG JS GS 

Instrument and Intercomparison 



Instrument and Intercomparison 

Fig. 15 (cont.). The average V3 processing uncertainties for the JG. JS, and GS intercomparisons of c) 
R rs and d) [Lw(A)] n - The dashed line sets the 2.5% intercomparison objective and the dotted line the 
1.0% objective. Uncertainties above 25% are shown clipped at the 25% level, and the maximum uncertainty 
achieved is given at the top of the panel to the side of the clipped bars. The individual wavelengths are shown 
as the sequential bars with varying intensities of gray (going from left to right, blue is light gray and red is 
dark gray). The instrument (and deployment location) codes are given along the top of the subpanels. The 
overall average of all the data (40 casts) is given in the right-most three intercomparisons. 
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data, wherein the Ed(0~) uncertainties normally increase 
as the chlorophyll concentration decreases. 

The V5 results show that differences in filtering tech- 
niques can account for a significant percentage of the un- 
certainty budgets. In many cases, the increases in uncer- 
tainties due to filtering differences was approximately one 
half (or more) of the uncertainty levels achieved when the 
processors were the most similar (V3). 

The overall sensitivity of the primary variables to chlo- 
rophyll a concentration depends primarily on the process- 
ing options. As the JG processors converged, any depen- 
dence on concentration weakened except for the aforemen- 
tioned larger Kd uncertainties in clear waters (the OT and 
even the MT regimes). 

In most instances, the band-ratio uncertainties are 
much less than the spectral uncertainties. Furthermore, 
once the processors were made more similar (i.e., con- 
verged), almost all of the band-ratio uncertainties were 
less than 1.0%. The largest band-ratio uncertainties are 
associated with the K d results, which are maximal in the 
blue-green wavelength domain. 

1.8 CONCLUSIONS 

From the strict point of view of data processing, the 
most obvious conclusion from this study is that when sep- 
arate groups process data in ways they consider appropri- 
ate, this can lead to unacceptable uncertainties (greater 
than 2.5%) in the final products used for calibration and 
validation activities. Although in situ data are probably 
the best approach for vicarious calibrations, it is important 
to note the uncertainties associated with the products de- 
rived from in situ data can be relatively high if different 
data processors are used. 


There are several other aspects of this study that are 
relevant to calibration and validation activities: 

1. The persistently lower L u (0~) uncertainties with re- 
spect to the [Lw(A)] n uncertainties suggest the 
validation of satellite radiometric data and the vi- 
carious calibration of space sensors might be best 
accomplished with Lvv(A) data collected at a time 
very close to the satellite overpass (which may be 
used without any normalization). 

2. Although a database constructed with processed 
data from a wide source of contributors will have 
substantially higher uncertainties than a database 
constructed with raw data processed with a single 
processor, the OC2v2 results are not very depen- 
dent on the processing options, because band ratios 
regularly provide reduced uncertainties with respect 
to individual spectral uncertainties. 

3. Data sets for mean F 0 (A) used in the computation of 
[L w { A)] N should be chosen in agreement with those 
used for satellite data processing to remove system- 
atic uncertainties in vicarious calibrations and in 
product validations. 

4. Because of the large uncertainties associated with 
the selection of extrapolation intervals, techniques 
should be developed to support the objective selec- 
tion of the limits for extrapolation intervals, as an 
attempt to standardize the selection of parameters 
used for the determination of subsurface values. 

5. Data collection from stabilized (i.e., fixed or gim- 
baled) platforms should be emphasized, and the ver- 
tical resolution of the in-water instruments should 
be maximized (either with a low descent speed or a 
high data rate). 
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Abstract 


The GSFC data processing system was designed to calculate apparent optical properties from underwater profiles 
of upward and downward irradiance, as well as upwelling radiance, from a wide variety of open ocean regimes. 
The full processing involves three steps: a) the calibration of the data, b) the establishment of the extrapolation 
interval, and c) the derivation of the final data products. In the first step, special files containing dark voltages 
are used to remove detector offsets and apply a pressure tare to the depth sensor. The second step is done 
interactively while viewing all the optical and ancillary (sensor tilts, fluorescence, etc.) fields. The final step is 
usually executed in batch mode using a file of extrapolation limits created in the second step. All of the code 
is written in the IDL programming environment. 


2.1 INTRODUCTION 

The GSFC processing software was developed to sup- 
port the activities of the SeaWiFS Field Team whose main 
objectives are a) the vicarious calibration of the space- 
borne sensor, b) the validation of the algorithms used to 
generate the SeaWiFS data products, and c) the evalu- 
ation of the instruments and protocols used for in-water 
radiometric measurements. Between September 1995 and 
April 1999, most of the field campaigns conducted by the 
SeaWiFS Field Team were performed during AMT cruises 
(Aiken et al. 2000). The codes described here were initially 
developed to process the AMT radiometric data for Sea- 
WiFS calibration and validation purposes. The processing 
system was designed to produce, in near-real time, all the 
quantities (direct or indirect) that can be derived from ra- 
diance or irradiance optical profiles conducted in the water 
column. 

During the AMT cruises, the SeaWiFS Field Team de- 
ployed a variety of commercial and prototype optical sam- 
pling systems built by Satlantic, Inc. These instruments 
included winch and crane, as well as free-fall deployment 
platforms, with 2 or 3 optical sensors on each system, and 
7 or 13 channels in each sensor. The former are capable 
of sampling the closest to the surface (especially during 
the up cast when the instruments are emerging from the 
mostly undisturbed deep part of the cast), and the sensors 
are located at very nearly the same depth: the latter cannot 


return near-surface measurements for all sensors because 
of the tilting of the profiler during its initial righting pe- 
riod immediately after release, and the sensors at each end 
of the profiler are usually separated from one another by a 
vertical distance of approximately 1-2 m. The primary wa- 
ter type sampled during the AMT cruises was clear, deep 
ocean, Case-1 conditions, although a large range (approx- 
imately three decades) in chlorophyll a concentrations was 
sampled after six cruises. 

The diversity in sampling equipment, a preponderance 
of deep casts (as much as 200 m) in clear water (Case-1 
conditions), the availability of coincident solar reference 
data, and the high quality of the in-water data were the 
governing forces behind the design of the GSFC processor. 
There was no need for a self-shading correction capability, 
a bottom correction scheme, or any other data handling 
requirements usually associated with the coastal (predom- 
inantly Case-2) environment, so these were not included as 
processing options. 

For the sake of efficiency, the same core of code was 
adapted for all instruments and configurations, which also 
share a variety of utility programs. The processing code is 
written in IDL and has three main components: 

■ A calibration module for the conversion of the raw 
data (counts) to physical units, which also generates 
an ASCII file containing the calibrated data for each 
cast plus a vertically binned file; 
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■ A visualization module that plots the various data 
fields as a function of depth; and 

■ A module that generates the surface values, that is, 
it extrapolates the underwater data to the surface 
and computes quantities relevant to ocean color and 
marine optics studies. 

For the DARR-00 exercise, some minor modifications were 
made to the codes, so the processing system could handle 
the formats of the JRC data. 

2.2 DATA COLLECTION 

The emphasis during SeaWiFS Field Team deployments 
is to use instruments with absolute calibrations, traceable 
to the National Institute of Standards and Technology 
(NIST), and state-of-the art sampling and optical mea- 
surement protocols. For each deployment, every effort is 
made to limit, a s much as possible, experimental (ship 
shadow, superstructure reflections, instrument instability, 
etc.) and environmental (sea and sky state) perturbations 
in the measurements. Under most circumstances, the data 
were collected in Case-1 waters (Morel and Prieur 1977) 
during stable sky conditions. To ensure a full optical char- 
acterization of the euphotic zone, the profiles were usually 
deployed to a depth that included the 1.0% light level and 
the deep chlorophyll maximum. 

For each cast, two instruments are used: one collects 
the in-water data in the euphotic zone as a profiling sys- 
tem, and the second one simultaneously records the above- 
water incident solar irradiance. The free-fall profiling sys- 
tems are deployed far enough from the ship to avoid ship 
perturbations, the winch system uses a crane with suffi- 
cient reach over the side of the ship to significantly re- 
duce contamination of the in-water data from ship per- 
turbations, and for both systems, the surface reference is 
installed in a location chosen to minimize the possible in- 
fluence of the ship’s superstructure. 

In addition to the in situ measurements, dark signals 
are recorded at least once a day for all light sensors. Dark 
levels for sensors are generally recorded for 3 min before or 
after a series of in-water profiles. While on deck, the pro- 
filer is kept under a white sheet that is regularly rinsed off 
to ensure the temperature is more similar to that encoun- 
tered during sampling. The dark data for the pressure 
sensor is also used as a pressure tare to ensure accurate 
depth readings during the in- water cast. 

The radiometric measurements always include the a) 
upwelled radiance, L u (z,\)\ b) downward irradiance, 
Ed(z,X)\ and c) incident total (direct plus diffuse) solar 
irradiance, Ed{ 0 + ,A). For most cruises, a third sensor is 
used to measure the upward irradiance, E u (z,X). Ancil- 
lary in-water measurements are used to ensure the quality 
of the data or to understand the light measurements in 
terms of the basic water column properties. These extra 
measurements include the two- axis tilts of the sensor, flu- 
orescence, plus water temperature and conductivity. 


2.3 DATA FORMAT 

The SeaWiFS Field Team has developed its own data 
acquisition software, so the format of the raw data are 
different from that of the optical equipment manufacturer. 
The raw data (profile, surface reference, and dark measure- 
ments) are stored in ASCII, tab-delimited files composed of 
a header section followed by a line-by-line data block con- 
taining the digitized measurements (Fig. 16). Keywords 
are used at the start of each header entry which are termi- 
nated with the at-sign (@). This architecture follows the 
design used with SeaBASS, which is the primary in situ 
database for the SeaWiFS Project (Hooker et al. 1994). 

The header section includes a number of values associ- 
ated with the configuration of the data acquisition module, 
the so-called DATA- 100, plus a variety of entries which 
identify the optical sensors used, experimental informa- 
tion, calibration documentation, etc. Note the use of an 
explicitly decoded Gregorian date — the month is spelled 
out — to prevent any time confusion. The data block be- 
gins with a header line and is followed by the sequentially 
numbered data lines. Explicit keywords for the beginning 
and ending of the data block are not needed, because it 
starts after the /end_header(9 keyword, and ends with the 
(platform-dependent) end-of-file (EOF) character. 

A data line begins with a sequential record number 
and a Greenwich Mean Time (GMT) time stamp based on 
the SDY. The format for the latter is yyyy dddhhmm ss . ss, 
where yyyy is the year, ddd is the SDY, hh is the (24-hour) 
hour of the day, mm is the minutes of the hour, and ss . ss 
is the seconds of the minute (to within 0.01s). The digi- 
tized sensor counts appear after the time stamp, which are 
followed by three quality control variables: Frame count 
is a sequential (circular) counter provided by the analog- 
to-digital converter unit, Checksum, and BadChecksum pa- 
rameters. 

2.4 PROCESSING STEPS 

For most applications, the raw data are considered as 
the lowest level, the so-called level-0 data. The basic archi- 
tecture of the GSFC processor is to convert these data to 
calibrated values, the so-called level-1 data. The level -1 
data are then used to calculate the final data products. 
Typically, four raw data files are used to process a cast: the 
in-water data, the solar reference data, the dark measure- 
ments of the in- water profiler, and the dark measurements 
of the surface reference for that cast or station. 

2.4.1 Data Conversion 

The first step of the GSFC data processing system is 
the transformation of the raw data into physical units for 
all measured parameters. For each parameter, the stability 
of the dark measurements is controlled (coefficient of vari- 
ation, i.e., standard deviation divided by the mean, less 
than 0.01% for the radiometric data) and the mean dark 
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/begin Jieader {cr} 

/sat_pref ix=SAT(cr) 

/ instrument_id=OCP(cr) 

/ serial_no=0011 (cr) 

/ samples_avg=l(cr) 

/ start_date=1999MAY10(cr) 

/ af f iliations=GSFC(cr) 

/ invesigators=STANF0RD_B_H00KER(cr) 

/ experiment=AMT_8 (cr) 

/ data_type=CAST(cr) 

/ column_headers=Sample_No , GMT.Ti me .Stamp , Edl , ... , Lul , ... , Eul , ... ,Temp_irr , Pressure , Tilt _X , Tilt _Y , 
Temp.w , Cond , Fluor , FrameCount , Checksum , BadChecksum (cr) 

/sequence_number=10 (cr) 

/ad_port_l=0CI_ED_050 (cr) 

/ad.port _2=0CR_LU_036 ( cr) 

/ad_port_3=0CI_EU_048(cr) 

/end _header@( cr) 

Sample_No GMT_Time_Stamp Edl . . . Lul . . . Eul Temp.irr Pressure Tilt JC Tilt.Y 
Temp.w Cond Fluor FrameCount Checksum BadChecksum (cr) 


1 

1999 130 1414 36.72. 

. .44295. 

. . 50021 

50441 

40064 

32233 

32696 

32459 

37013 

56515 

33350 

54 

62 0 (cr) 

2 

1999 130 1414 36.89 

. .41733. 

. . 50340 

50421 

40065 

32258 

32927 

32500 

37017 

56519 

33356 

55 

68 0(cr) 

3 

1999 130 1414 37.05. 

. .42239. 

. . 50774 

50401 

40064 

32285 

33133 

32585 

37023 

56525 

33354 

56 

78 0 (cr) 

4 

1999 130 1414 37.20. 

. .40577. 

. .51217 

50429 

40065 

32312 

33254 

32737 

37026 

56525 

33351 

57 

220 0 (cr) 

5 

1999 130 1414 37.36 

. .43097. 

..51209 

50429 

40066 

32341 

33390 

32923 

37028 

56527 

33355 

58 

97 0 (cr) 

6 

1999 130 1414 37.56. 

. .46355. 

. . 50924 

50403 

40066 

32369 

33547 

33105 

37029 

56530 

33364 

59 

100 0 (cr) 

7 

1999 130 1414 37.72. 

. .48256 . 

. . 50629 

50382 

40066 

32395 

33643 

33230 

37032 

56531 

33356 

60 

2 0 (cr) 

8 

1999 1301414 37.87. 

. .48696. 

..50102 

50347 

40067 

32422 

33698 

33349 

37034 

56533 

33361 

61 

30 0 (cr) 

9 

1999 130 1414 38.02. 

. .50337. 

. .49611 

50220 

40067 

32449 

33762 

33428 

37034 

56534 

33364 

62 

233 0 (cr) 


Fig. 16. An example SeaWiFS Field Team LoCNESS data file. The /beginJtieader keyword denotes the 
start of the header section which ends with the /end_header@ keyword. The first line after the latter is a 
header line for the data block, and all subsequent lines are data lines. Explicit carriage returns, (cr), are 
shown, because some lines are displayed as wrapped entries (they continue onto a following line), and others 
as abbreviated entries (they contain the . . ” sequences, which indicate omitted values). 


value is computed. The calibration files for each of the 
sensors involved in a particular instrument are read, and 
the relevant calibration coefficients are extracted. 

The raw data are converted to physical units based on 
a formulation given by the manufacturer of the sensor. For 
example, the fth sample for a radiometric sensor at center 
wavelength A, are converted to physical units according to 
the following equation: 

©(A ,U) = C c (\) //(A) [V (A, ti) - D(X)], (4) 

where 0(A ,t z ) is the calibrated value from the GSFC pro- 
cessor, for which 0 is replaced by Ed(z), E u (z). L u (z). 
or Ed{ 0 + ); C c (A) is the calibration coefficient; //(A) is the 
immersion factorj; V(\.t l ) is the raw voltage (in digital 
counts) measured by the instrument at time ti (which also 


f The coefficient accounting for the change in sensor respon- 
sivity when the in-air calibration is applied to in-water mea- 
surements {If — 1 for above- water sensors). 


sets the depth, z); and 5(A) is the average dark value (in 
digital counts) measured during a special dark cast with 
the caps on the radiometer. 

Of the ancillary data available with the SeaOPS and 
LoCNESS profiling systems, the tilt of the profiler with 
respect to the vertical, <p, is important in ensuring the 
quality of the optical data: 


ip = tan 1 ^/tan 2 (p x T tan 2 (p y (5) 

where <p x and ip y are the x- and y - axis tilts, respectively, 
with respect to the vertical coordinate. 

The dark calculations and the transformation into phys- 
ical units are done for the profiler and the reference suc- 
cessively. The two calibrated data sets are then combined 
into a unique ASCII file. The matchup between the in- 
water and reference data are based on the time stamp of 
each sensor record which must agree within a time interval, 
St, defined by the user. The default value for St is set for 
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an agreement of less than 0.1 s. This value is based on the 
sampling frequency of the light sensors which is 6 Hz, thus, 
each sample is generated and acquired within 167 ms, so 
a 6t value of 0.1 s is approximately one-half the sampling 
interval. 

An ASCII file containing the entire (unbinned) profile 
is then generated. It includes depth, radiometric measure- 
ments from the profiler and the reference, and all other 
available ancillary data. With the exception of the tempo- 
ral match-up filter between the in-water and surface ref- 
erence data, no other filtering is applied to the unbinned 
level -1 data. After the combined (unbinned) file is com- 
pleted, a vertically binned file is also generated. The ver- 
tical bin size, Sz, is defined by the user (the default value 
for 5z is 1 m). For this binned file, filtering of the data can 
be applied based on the two-axis tilts of the instruments 
during the cast. Threshold values for these angles can be 
set by the user resulting in the rejection of the data line if 
the threshold value is exceeded (the default tilt threshold 
is 10°). 

2.4.2 Data Visualization 

Although the profiles can be optionally displayed dur- 
ing the calibration step, the GSFC processing system has 
a specific data visualization module that allows the cali- 
brated data to be plotted for whatever depth range the 
user specifies. The default graphic window displays six 
panels: 1) the L u (z, A) profile, 2) the E d (z, A) profile, 3) 
the Ed{ 0 + , A) data (plotted as a function of time, but cor- 
related with the corresponding in- water measurements), 
4) the two-axis tilt profiles, 5) the temperature and con- 
ductivity profiles, and, 6) the fluorescence profile. Other 
combinations of plots are possible and user definable. 

The visualization module permits a visual inspection 
of all the measured parameters. Because the data can be 
displayed for any desired depth interval, the visualization 
module is very useful for defining what depth range should 
be used to perform subsurface extrapolations to the sur- 
face, which is required to determine the values of some 
parameters at null depth, z — 0". In the GSFC processing 
system, the selection of these depth ranges by the data an- 
alyst is a necessary step before any extrapolations can be 
achieved — there is no automated method for such a selec- 
tion. 

The selection of the data is based on the examination of 
the upper part of the profile (0 to 5-30 m depending on the 
clarity of the water). The data analyst selects the depth 
range to be used in the extrapolations (defined by an upper 
and a lower depth threshold) based on the characteristics 
of the radiometric and ancillary data. The stability of the 
E d (0 + , A) data (no clouds or other perturbations), the lin- 
earity of the L u (z, A) data (in semi-log space), two-axis 
tilts, the characteristics of the mixed-layei depth, and the 
features of the fluorescence profile are all taken into ac- 
count in that selection. There is some obvious subjectivity 


in this part of the process, because two different analysts 
might interpret the profiles in different ways, which would 
probably result in different depth threshold selections. 

A completely objective or automated procedure in the 
depth selection process is almost impossible to achieve and, 
in any case, is not considered desireable for two reasons: 
a) in the vicarious calibration process, each in situ ob- 
servation is potentially very precious, because only a few 
matchups are usually derived from each week of field work 
(Hooker and McClain 2000), so each cast must be pro- 
cessed to the highest standard possible; and b) automated 
procedures tend to produce biases or nonstochastic offsets 
in the data, which would add an unnatural variance to the 
algorithm validation database. 

Generally, the data in the first 1-2 m below the surface 
cannot be safely used, because of excessive sensor tilts and 
surface contamination effects (bubbles, foam, wave focus- 
ing, etc.). When correctly trimmed, free-fall profilers gen- 
erally reach their vertical stability after 1-2 m of descent. 
Aside from the characteristics of the E d ( 0 + ,A) data, the 
lower depth to be used in the extrapolation range depends 
mostly on the clarity of the water. In clear water, the lower 
limit can be as deep as 20 m, while in more turbid waters, 
this limit can be very shallow (less than 4 m). 

One limitation of the GSFC data processing system 
is that it was not designed to generate accurate results for 
wavelengths which are strongly influenced by water absorp- 
tion (i.e., A > 570 nm). When water absorption is strong, 
the red bands require very specific, and generally very shal- 
low, depth ranges for the extrapolations to be meaningful. 
Because the SeaWiFS validation effort was initially focused 
on the first five SeaWiFS bands, no specific processing was 
developed for the red channels. Furthermore, as good data 
are generally collected once the L u sensor is 2-3 m deep in 
the water column, this may frequently be too deep to al- 
low a relevant extrapolation of the red channels. The null 
depth values of the red channels are, therefore, acknowl- 
edged to be less accurate than the blue and green bands. 

2.4.3 In- Water Data Products 

Once the depths for the extrapolation intervals have 
been selected, the values at null depth are calculated for 
the L u (z, A), E u (z, A), and E d (z,X) data. This step also 
involves the computation of various other derived param- 
eters above or beneath the sea surface interface. These 
parameters, and how they are calculated, are described 
below. 

The first step of the in--water computations consists of 
a series of regression calculations performed, at each wave- 
length, between the depth (corrected for any difference 
between the position of the pressure sensor and that of 
the individual light sensors) and the log-transformed (nat- 
ural logs) radiometric data. These regressions generate 
the L u (0“,A), E d ( 0",A), and E u ( 0‘, A) values, as well as 
the diffuse attenuation coefficients for the subsurface re- 
gion (the extrapolation range): a) downward irradiance, 
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K Ed (zo, A); b) upwelled radiance. Kl u { z o?A); and c) up- 
ward irradiance, K £ v (z 0: A). 

This approach assumes a good linear relationship be- 
tween the log-transformed (natural log) radiometric data 
and the depth, such that 


R rs (A), that can be considered independent of the varia- 
tions of the incident solar radiation: 


^rs(A) 


Lw(X) 

E d ( 0 + ,A)‘ 


( 10 ) 


In [ 0 ( 2 , A)] = In [©(so, A)] - K e (zo,X)dz, (6) 


where, again, 0 is either L u , E d , or E u and dz — z ~ zq. 

The slope of the regression provides the diffuse atten- 
uation coefficients, while the exponential of the intercept 
gives the L u ( 0". A), ^(O - , A), and £^(0”, A) values. By de- 
fault, the processing program uses a classic least-squares 
method for the regressions, but a robust least absolute de- 
viation method (Press et al. 1992), can also be used. For 
the latter, possible outliers (noisy points with a larger vari- 
ance) have a reduced weight in the analysis. 

The subsurface Q-factor at null depth derived from the 
in- water (nadir-viewing) instruments, Q n (0",A), is com- 
puted as 


Qn{ OTA) - 


Eu( OTA) 
L u ( OTA)' 


(7) 


Another commonly used quantity is the normalized 
water-leaving radiance defined for in situ measurements 
as: 

[Lw( A)] n = F 0 (A)/? rs (A), (11) 

where Fo is the mean extraterrestrial solar irradiance cor- 
rected for the Earth-sun distance which depends on the 
SDY, d, i.e., F 0 = Fo{X)Sq(cI). The Earth-sun distance 
correction is computed as 


So(d) = 


1.0 + 0.0167 cos 


27 r(d - 3) 
365 


( 12 ) 


The Fo(A) values used with the GSFC processor are in- 
terpolated from the Neckel and Labs (1984) data set, but 
weighted by the spectral response of the SeaWiFS bands. 


2.4.4 Above- Water Data Products 

The L u { 0", A) value is transmitted through the sea sur- 
face to derive the upwelled radiance leaving the sea surface, 
the so-called water- leaving radiance, L\\;( A), according to 
Austin (1974): 

Lw{ A) - Lu ( 0-,A), (8) 

7 Rl(A) 

where p{ A, $) is the Fresnel surface reflectance at the x) 
nadir angle and n w (X) is the refractive index of seawater 
(all the in-water instruments use nadir-viewing radiance 
sensors for which *0 = 0°). 

Except for a very rough sea surface or a high sun zenith 
angle, the transmission term in (8), (l — p(X.d))n~ 2 (X), 
varies only slightly in magnitude and wavelength (Austin 
1980). The coefficient 0.54 has been shown (Mobley 1999) 
to be appropriate under most situations for transmitting 
the normal radiance from below, to above, the sea surface 
(Mobley 1999), so (8) is rewritten as: 


2.4.5 Ancillary Data Products 

The extrapolation module also calculates several other 
parameters: 

1. The solar zenith angle (computed from the average 
time measured during the extrapolation interval): 

2. The thickness of the mixed-layer depth (as deter- 
mined by a drop of 0.1°C in four consecutive lm 
depth bins) and the mean mixed-layer depth tem- 
perature: 

3. The surface PAR values; 

4. The depths of the 10%, 1%, and 0.1% surface PAR 
values; 

5. The depth of the chlorophyll maximum (determined 
from the fluorescence profile); and 

6. The surface chlorophyll a concentration as estimated 
by the OC2 and OC4 algorithms. 

For any particular cruise or experiment, the surface data 
and other derived quantities are saved in ASCII files. 

2.5 SUMMARY 


Lw{ A) = 0.54 L u (0~ , A). (9) 

The only directly measured above- water parameter is 
the incident solar irradiance. The stability of the light 
field for the time period matching the chosen extrapo- 
lation range is characterized by the coefficient of varia- 
tion in the corresponding E d ( 0 + ,A) values. The spectral 
means of the E d ( 0 + ,A) values for the extrapolation time 
period, E d { 0 + , A), are used for the normalization of Lw( A) 
which results in a quantity, the remote-sensing reflectance, 


The GSFC data processing system has been developed 
to analyze the radiometric and ancillary data collected for 
calibration and validation purposes by the SeaWiFS Field 
Team. The processor includes 3 major steps: a) the conver- 
sion of the raw data to (calibrated) physical units, b) the 
visualization of the calibrated data and, c) the calculation 
of the surface values for ocean color quantities and other 
ancillary data. The codes can process data collected by 
S at 1 antic OCR-200, OCI-200, OCR-1000, and OCM000 
sensors. Several options can be selected throughout the 
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processing procedure: tilt filtering threshold, data visu- 
alization scaling, vertical bin size, extrapolation interval 
depth range, and the type of regression used for the ex- 
trapolation. The final step of the processing generates 
the following products: L u (0~,A), Ed{ 0~, A), £<u(0~,A), 
E d ( 0 + ,A), K d ( A), K l ( A), K u ( A), R rs (\), [L w ( A)] N , and 
Q n (0 , A). Several ancillary data (temperature, conduc- 
tivity, fluorescence, etc.) are used to estimate PAR, MLD, 
DCM, sun zenith angle, etc. The processor is written in 
the IDL programming environment. 
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Abstract 

The JRC Data Processor was developed for the analysis of field data collected in coastal and marine regions 
in support of ocean color calibration and validation activities. The code specifically supports the calibration, 
formatting, visualization, and processing of in-water optical profiles of L u (z , A), Ed{z. A), and E u (z. A) collected 
with the Satlantic seven-channel OCR-200 and OCT200 sensor series. The user can define a variety of options 
which permit the selection of the extrapolation interval for the computation of the near-surface values, and 
the removal of outliers through statistical filtering. Relevant features include the possibility to minimize mea- 
surement perturbations such as instrument self shading, bottom effects, and superstructure shading (although 
the latter is restricted to measurements performed at the Acqua Alta Oceanographic Tower site in the northern 
Adriatic Sea). Each step of the processing is supported by the graphic presentation of intermediate results. Out- 
put products are the subsurface values L^(0“,A), Ed( 0",A), and E u ( 0",A); the diffuse attenuation coefficients 
K d ( A), Kl{ A), and K u { A): plus the derived products R( 0", A), R rs (A), [Lw(A)] n , and Q n (0“,A). 


3.1 INTRODUCTION 

The JRC Data Processor is a comprehensive software 
package for field data analysis. It was developed within 
the Coastal Atmosphere and Sea Time-Series (CoASTS) 
project (Zibordi et al. 2001) to support ocean color cal- 
ibration and validation activities. The processor ensures 
handling of marine and atmospheric measurements and, 
specifically, the analysis of in-water optical profile data. 

An important feature of the processor is the capabil- 
ity for removing measurement perturbations, like instru- 
ment self shading, bottom effects, and tower shading [the 
latter correction is specific to data collected at the Acqua 
Alta Oceanographic Tower site (Zibordi et al. 2001)] . This 
makes the code well suited for the processing of oceanic, as 
well as coastal, data. The processor is composed of a data 
calibration and formatting program plus a data filtering 
and processing program, and uses a graphical user inter- 
face (GUI) to interactively define the processing settings. 

The data calibration and formatting program has five 
main functions: 

1. Ingestion of ASCII raw (i.e., level -0) data in differ- 
ent formats; 

2. Calibration of the data, including dark correction 
for light sensors and pressure tare for depth sensors; 


3. Screening of data as a function of a user-selectable 
tilt threshold; 

4. Visualization of data as a function of depth; and 

5. Creation of the level -1 (output) data files. 

The DARR-00 exercise also required the implementation of 
specific routines for ingesting the GSFC ASCII raw data. 

The data filtering and processing program has six main 
functions: 

i) Definition of extrapolation intervals for the compu- 
tation of subsurface values and bottom reflectances; 

ii) Selection of normalization methods for in-water 
data with respect to above- water reference measure- 
ments; 

iii) Filtering of outlier data, in the subsurface extrapo- 
lation interval, and removal of data below the noise 
threshold (defined as the radiance or irradiance 
given by 5 digital counts); 

iv) Correction for perturbation effects; 

v) Visualization of data as a function of depth together 
with extrapolated (surface and bottom) values; and 

vi) Creation of the level -2 (output) binned data and 
products. 


52 


S. Hooker, G. Zibordi, J-F. Berthon, D. D’Alimonte, S. Maritorena, S. McLean, and J. Sildam 


All of the output files are in ASCII format, and can be 
quickly and easily used by the user. 


3.2 DATA REQUIREMENTS 

The JRC optical processor was primarily designed to 
support the analysis of data collected by Satlantic optical 
systems composed of multiple seven-channel OCR-200 and 
OCT200 sensors to measure L u (z, A), Ed(z, A), E u (z, A), 
and Ed{ 0 + ,A). The analysis of in- water profile data re- 
quires the following ancillary data: 

a) Dark data for the calibration process (including set- 
ting the pressure tare); 

b) Seawater absorption, a(A); 

c) Seawater single scattering albedo, u;(A); and 

d) The diffuse- to-direct solar irradiance ratio, r(A). 
The latter is obtained from 


r{ A) = 


Ej{ 0 + ,A) 

E d (0 + ,A)-E t (0+,A)’ 


where 3(A,^) is the calibrated value from the JRC pro- 
cessor, for which 3 is replaced by Ed(z ), E u (z ), L u (z), or 
£d(0 + )r; C c (A) is the calibration coefficient; //(A) is the 
immersion factor; V(X ,U) is the raw voltage (in digital 
counts) measured by the instrument at time U (which also 
sets the depth, z); and D( A) is the average dark value 
(in digital counts) measured at least once a day during a 
special dark cast with the system in air and caps on the 
radiometer. 

3. 3. 1.3 Data Screening 

An optional capability for screening the data as a func- 
tion of the instrument tilt (i.e., removing data above a tilt 
threshold value which can be set by the user) is computed 

from ip = tan -1 y^tan 2 <p x + tan 2 (p y , where p x and <p y are 
the two-axis tilts with respect to the vertical coordinate. 
The threshold for data rejection, </? 0 , is user-selectable and 
has a default value of </? Q — 5° . 


where F z (0 + , A) is the diffuse solar irradiance (for the cor- 
rection of perturbation effects). 

3.3 DATA PROCESSING STEPS 

Data processing includes two steps: i) data calibra- 
tion and formatting to convert the raw data from counts 
(level -0 data) to physical units (level -1 data); and ii) data 
filtering and processing to remove outliers and create the 
level -2 data products. During each step, for each optical 
profile, the user can define processing parameters, which 
are stored in an information file and can be reloaded for 
subsequent automated reprocessing. 

3.3.1 Calibration and Formatting 

Data calibration and formatting involves five steps: in- 
gestion, calibration, screening, display, and output. 

3. 3. 1.1 Data Ingestion 

Input data for the calibration and formatting program 
include several ASCII files: i) the in-water radiance and ir- 
radiance files; ii) the reference above-water irradiance file; 
iii) the calibration files for the in-water and above-water 
devices (i.e., radiometers, pressure sensors, and tilt sen- 
sors); and iv) the dark data for the in- water and above- 
water radiometers. Specific ingestion routines ensure in- 
gestion of profile and reference data in different formats. 

3. 3. 1.2 Data Calibration 

The absolute calibration of the primary radiometric 
quantities, L u (z, A), Ed(z,X), E u (z, A), and F^(0 + ,A), are 
obtained from: 

3(A ,U) = C c (A)If(A)[V(A,t i ) - 5(A)], (14) 


3. 3. 1.4 Data Display 

Calibrated profiles, i.e., L u (z, A), Ed{z, A), and 
E u (z, A), can be displayed, together with tilts as a func- 
tion of depth in logarithmic scale. Ed{ 0 + , A) is displayed 
as a function of the depth of the coincident in-water pro- 
file data. This easily ensures the identification of artifacts, 
which may affect the optical profile data. 

3. 3. 1.5 Output Data Format 

The output data are stored using the Coastal Region 
Long-Term Measurements for Colour Remote Sensing De- 
velopment and Validation (COLORS) data formatf which 
is composed of four blocks of information: 

a) The header block contains geographic and tempo- 
ral information, in addition to codes identifying the 
methods used during data collection; 

b) The calibration block contains information on the 
measurement systems, the sensors heads, and the 
calibration coefficients for each sensor; 

c) The data block description provides the list of quan- 
tities in the data block and the units assigned to 
each of them; and 

d) The data block contains the in- water data as a func- 
tion of depth together with the related above-water 
reference data, tilts, and time stamps. 

An example of the COLORS data format, for optical pro- 
file data, is given in Fig. 17. 

f The COLORS data can be reviewed on the World Wide Web 
at http : //marine . ie/datacentre/pro jects/ colors/ . 
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/begin_header( cr) 

/ site_code=P(cr) 

/method_code=04A005 , 04A004 , 04A003 , 04A003(cr) 

/parameters=subsurf ace_upwelling_spectral_radiance , subsurf ace_downwelling_spectral_ir radiance , 
surf ace_downwelling_spectral_irradiance , subsurf ace_upwelling_irradiance (cr) 

/latitude [deg] =45 . 3139(cr) 

/longitude [deg] =12 . 5083(er) 

/mean.sounding [m] =17 . 00(cr) 

/sampling_date=121 1 1999 (cr) 

/ start_time [GMT] =09 : 11: 15{cr) 

/end_time [GMT] =09 : 11 : 36 (cr) 

/station_ref=P630302P( cr) 

/calibration_block(cr) 

instrument _id=SATLANTIC 0CP#016 , SATLANTIC MVD#054(cr) 
calibrat ion_date=0107 1999 (cr) 

Lu_s ensor _id=067 ( cr) 

Ed-sensor_id=097 (cr) 

Es_sensor_id=129 (cr) 

Eu_sensor_id=109(cr) 

Lu_depth_of f =-0 . 098 (cr) 

Ed-depth_of f =-0 . 258(cr) 

Es-depth_of f =0 . 0(cr) 

Eu_depth_of f =0 . 763(cr) 


! Variable 

Depth 

LU412 

Lambda 

-999.9 

412.5 

CoeO 

32588.7 

32775.6 

Coel 
6 . 4830e-003 
1.7770e-006 

Coe2 

0.00 

1.75 

Coe3 

0 

0 

Coe4 (cr) 
0(cr) 
0 (cr) 

ED412 

412.3 

32783.0 

6. 1497e-005 

1.52 

0 

0(cr) 

ES412 

412.1 

32779.1 

9 . 6073e-005 

1.00 

0 

0 (cr) 

EU412 

412.5 

32770.5 

4 . 1442e-006 

1.52 

0 

0(cr) 

EU685 

683.2 

32774.0 

2 . 2731e-006 

1.38 

0 

0{cr) 


/datablock.description(cr) 

/fields=seq_number, sampling-depth [m] , sampling-time [hh : mm : ss : ssss] , SUR-ATX [deg] , . . . (cr) 

/n_f ields=35 (cr) 

/n_records=94 (cr) 

/ def ault=-9 ,-9.99, -99 : -99 : -99 : -9999 ,-9.99,-9.99,-9.99, -9 .99,-9. 99999999 , . . . , -9 . 99999999(cr) 
/delimiter=comma(cr) 

! Related CoASTS f iles=p63d3a . ocp,p63s3b . ocp, p63d3a .mvd,p63s3b .mvd( cr) 

! Calibration f iles=ocpd. cal ,mvdd. cal (cr) 

/end_header (cr) 

/begin_data( cr) 

1.0. 21.09:11:19: 4060 ,0.00,0.00, 4.06,-1.30,0. 00685499 , . . . , 0 . 00735924(cr) 

2.0. 86.09:11:20:5680,0.00,0.00, 4.64, 0.58,0.00557378, . . . , 0 . 00507873(cr) 

3.0. 99.09:11:20:7380,0.00,0.00, 3.65, 2.18,0.00542140, ... , 0 . 00458938(cr) 
4,1.08,09:11:20:8980,0.00,0.00, 2.28, 2.65,0.00511975, ... , 0 . 00405297(cr) 
5,1.23,09:11:21:0680,0.00,0.00, 1.02, 2.08,0.00472793, . . . ,0 . 00356675(cr) 

6,1.37,09:11:21: 2390 ,0.00,0.00,-0.32, 1.11,0. 00436408 , . . . , 0 . 00312759(cr) 


/ end_data( cr) 
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Fig. 17. An example COLORS data file. The header section begins with the /begin_header keyword 
and ends with /end-header. Similarly, the data block begins with the /begin.data keyword and ends with 
/ end-data. 
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3.3.2 Filtering and Processing 

The primary purpose of the filtering options with the 
processor are to deal with data outliers and to correct for 
in situ perturbations to the light field. 

3. 3. 2.1 Extrapolation Intervals 

The appropriate selection of the extrapolation inter- 
val, is the key element in extracting accurate subsurface 
data from in-water profiles. In coastal water, where op- 
tical stratifications may occur in the first few meters of 
the water column, an inappropriate selection of the ex- 
trapolation interval can induce large uncertainties in the 
final products. Although default values can be defined, the 
identification of an optimum extrapolation range satisfying 
the requirement for a linear decay of In [3fy, A)] is made 
by the operator on a cast-by-cast basis. 

In addition to the subsurface extrapolation interval, the 
selection of a bottom extrapolation interval is required to 
estimate bottom reflectance from profile data, when the 
correction for bottom effects is needed. 

3. 3. 2. 2 Normalization of In- Water Data 

In-water radiometric quantities can be normalized with 
respect to Ed(0 + , A,t), with t explicitly expressing the de- 
pendence with time, according to 


3(z, A, to) 


3(z, A ,t) 


Ed( 0 + , A, t Q ) 
£ d (0 + ,A,t) ’ 


(15) 


where 3 ( 2 , A, to) identifies the radiometric quantity as it 
would have been taken at all depths z at the same time 
to; Fd(0 + ,A,t) is the solar irradiance taken at the same 
time, t, as the in- water data; and Ed( 0 + , A, to) is the above 
water irradiance at time to (where to is generally chosen 
to coincide with the start of the cast). 

Alternatively, F^(0 + ,A,t o ) can be defined as the me- 
dian for data records in the extrapolation interval. The 
normalization process can be skipped if data are not avail- 
able or are strongly affected by noise (e.g., ship motion). 
When Fd(0 + ,A,t) is affected by harmonic noise (induced 
by the regular motion of the ship), Ed( 0 + , A, t) can be lin- 
early fitted as a function of t, and the fitted values can be 
used in the normalization of the profile data. For simplic- 
ity, the variable t is omitted in the following discussion. 


3. 3. 2. 3 Filtering 

Data extrapolation in the subsurface layer is carried out 
using a least-squares linear fit. This technique may pro- 
duce biased results in the presence of outliers, e.g., pro- 
duced by wave focusing effects. This bias may become 
larger as the number of points per unit depth decreases in 
the optical profile. Because of this, two complementary fil- 
ters have been implemented. The first ensures removal of 


points distant by more than no from the linear regression 
line obtained with all profile data, where n is a user-defined 
multiplier, and 0 is the standard deviation in the distance 
of the data points from the regression line. 

The second filter computes successive regressions by 
adding depth increments to a defined basic extrapolation 
interval, and compares the slope of the resulting regres- 
sion line with the slope of the regression line obtained for 
the basic extrapolation interval. The latter filter may re- 
move all points below a certain depth within the chosen 
fitting interval, if for a given depth increment, the slope 
of the fitting line changes by more than a selected percent 
difference with respect to the initial value. 

Regression trials, to verify filtering in the extrapolation 
interval, are possible choosing a specific in-water quantity 
(i.e., L u or Ed) and a specific wavelength, e.g., £<*(665). 
The use of Ed also allows a direct comparison of the in- 
water extrapolated value £ c j(0~,A) with the solar refer- 
ence value, Ed( 0 + , A), to immediately identify difficult sit- 
uations (i.e., those showing very different Ed( 0",A) and 
£<*(0 + ,A) values). The use of a channel in the red part 
of the spectrum, where seawater is characterized by high 
absorption and light levels exhibit a rapid drop to noise 
levels as a function of depth, helps in excluding noisy data 
from the fitting interval chosen for processing. 

All records (i.e., all quantities and channels) taken at 
the same time as data flagged by the filter are removed. 
The assumption made here is that all data collected at the 
same time are affected by the same perturbations. The 
input parameters are user selectable, so the filtering tech- 
niques can be used while retaining complete control in the 
removal of the outliers. 

3. 3. 2. 4 Removing In Situ Perturbations 

Instrument self-shading, bottom reflection, and tower- 
shading effects can be removed through single or multiple 
selections. The correction factors are calculated assuming 
all the perturbations are independent from each other. 

Self-shading correction factors are computed using the 
scheme proposed by Gordon and Ding (1992) and the pa- 
rameterizations suggested by Zibordi and Ferrari (1995) 
and Mueller and Austin (1995). The computations are 
carried out using the solar zenith angle, 9 , as well as a( A), 
r( A), and the ratio between the instrument diameter and 
the sensor entrance optics diameter. 

The correction factors for bottom reflection effects are 
estimated with a simple analytical model (Zibordi et al. 
2001) adapted from Maritorena et al. (1994). Compu- 
tations are performed using: a) the bottom reflectance, 
p( A) = E u (z}>, \)/Ed(zi, A), obtained from each profile at 
a depth z\> near the bottom; and b) the mean diffuse at- 
tenuation coefficient £<*( A) for the depth interval z 0 — Zb, 
where z 0 is the subsurface depth identifying the upper- 
most value of the extrapolation interval. 

Tower-shading correction factors are obtained from 
look-up tables (Doyle and Zibordi 1998) derived from 
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Monte Carlo simulations of the radiance and irradiance 
fields at the AAOT site for the specific point where the 
WiSPER measurements are made (Zibordi et al. 1999). 
For each center wavelength, A, the correction factors are 
indexed by discrete values of 0, a( A), u;(A), and r( A). 


3. 3. 2. 5 Fitting 

Subsurface 3(0“ , A) optical values ( L u , E u , and E d ). are 
obtained from the exponent of the intercept given by the 
least-squares linear regression of ln[3(z,A)] (computed 
from data not removed by filters) versus 2 within the ex- 
trapolation interval identified by z$ < 2 < z\. Generally 
for the CoASTS profiles, 0.3 < zq < 1 m and 2.5 < z\ < 

4.5 m. The negative value of the slope of the regression fit 
is the diffuse attenuation coefficient A3 (A). 

An option, called “multicast processing” enables the 
processing of data from a set of profiles, collected within a 
user-defined time interval (generally spanning a few min- 
utes), as a single optical profile. This procedure increases 
the statistical significance of regression analysis applied to 
subsurface data characterized by a low depth resolution 
(which is a usual characteristic of free-fall instruments) 
and affected by large perturbations (e.g., caused by sur- 
face effects). 

3. 3. 2. 6 Data Visualization 

For each specific quantity and for all wavelengths, it 
is possible to display the optical data and accompanying 
tilts as a function of depth. Regression lines for the near- 
surface extrapolation layer and the bottom layer are also 
shown. In addition, arrows at 0 depth indicate the subsur- 
face extrapolated value, the subsurface extrapolated value 
corrected for perturbation effects, and in the specific case 
of E d ( 0",A), also the corresponding E d ( 0 + ,A). 


3. 3. 2. 7 In- Water Data Products 


The major in- water products obtained from optical pro- 
files are L u (0“,A), E d ( 0“,A), E u ( 0 , A), and K d ( A) com- 
puted from regressions in the extrapolation interval. Ad- 
ditional products are the subsurface Q-factor at nadir view, 
which is computed as Q n (Q~,^) — E u ( 0 , A)/L u (0 ,A), 
and the subsurface irradiance reflectance, R( CT,A), which 
is given by 


R{ O’, A) 


E u {0’, A) 
E d { 0“,A)' 


(16) 


3. 3. 2. 8 Above- Water Data Products 

The above- water products are a) the total solar irra- 
diance measured at the start of the extrapolation range, 
£ d (0 + ,A,t 0 ); b) the remote sensing reflectance, R rs ( A) = 
L w (\)/E d (tf , A, t 0 ); an d c) the normalized water-leaving 
radiance, [Lw(A)] n = F 0 (A) i? rs (A). For the latter, the 
mean extraterrestrial solar irradiance, Fq, is corrected for 


the Earth-sun distance which depends on the SDY, d , i.e., 
F 0 = Fo(\)$o(d) using 


60(d) 


1.0 + 0.0167 cos 


6.283(d - 3) 


365 


(17) 


The Fq(A) values used with the JRC processor are inter- 
polated from Neckel and Labs (1984) to the center wave- 
lengths of the sensors involved. 

Alternatively, if E d ( 0 + , A) is not directly measured (with 
a separate sensor), F rs (A) and [Lw(A)] n are derived using 
E d ( 0 + , A) values computed from E d (0~ , X) according to 
E d { 0",A) - 0.49 F n (0”, A) 


E d ( 0 + ,A) = 


1 


Ps (A) 


(18) 


where p s ( A) is the sea surface reflectance for downward 
irradiance. The latter is given by 


PW + 0.66 r(A) 
= 1 + r(A) 


(19) 


where p(6) is the Fresnel reflectance of the sea surface, and 
the constant 0.66 is the sea surface albedo under diffuse 
illumination. For very clear-sky conditions. r(A) 1 and 
Ps( A) -> p(0). 


3.4 SUMMARY 

The JRC Data Processing System has been designed 
to perform absolute calibration, formatting, visualization, 
and processing of in-water optical profile data. Data han- 
dling is specifically optimized for OCT200 and OCR-200 
S at 1 antic seven-channel radiometers. The processor, devel- 
oped to ensure the analysis of both coastal and open ocean 
data, has options for the following: 

1. The selection of the extrapolation interval chosen to 
satisfy linear decay of the logarithm of profile data 
at all wavelengths; 

2. The normalization of the in- water profile data with 
respect to the above-water downward irradiance: 
and 

3. The capability of combining multiple casts to in- 
crease the statistical significance of single profile 
data collected with low vertical resolution. 

Additional relevant processing options are: 

4. The filtering of outliers through statistical analysis; 
and 

5. The removal of measurement perturbations for self- 
shading, bottom effects, and superstructure shading 
(only for data collected at the AAOT site in the 
northern Adriatic Sea) making use of atmospheric 
parameters as the ratio between diffuse and direct 
atmospheric irradiance, and marine parameters as 
seawater absorption and single scattering albedo. 
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Chapter 4 


The ProSoft Optical Data Processor 

Juri Sildam 
Scott McLean 
Satlantic, Inc . 

Halifax , Canada 

Abstract 

The ProSoft data processor was developed for the analysis of optical data collected with instruments manu- 
factured by Satlantic, Inc., including profiling and fixed-depth (buoy) systems. The data processing steps are 
organized into different levels which follow, as much as possible, the procedures given in the Ocean Optics Proto- 
cols for SeaWiFS Validation, Revision 1 (Mueller and Austin 1995). The processing involves three steps wherein 
the level -0 raw data (usually binary) are converted to calibrated values (engineering units) using one of several 
options for dark current correction, the level -1 data are averaged and corrected for acquisition and performance 
problems (if any), and the level -2 data are used to compute the final data products, e.g., R rs , [Lw(A)] n , PAR, 
etc. Processing can be done manually (and interactively) by the user, guided by the ProSoft processing steps, or 
in a fully automated batch mode. A notable difference between ProSoft and other data processors is the diffuse 
attenuation coefficient can be estimated using a traditional estimation of the slope of the optical variables within 
a user-defined extrapolation interval or using a new methodology based on multiresolution wavelet analysis. 


4.1 INTRODUCTION 

ProSoft is a data analysis package for processing op- 
tical data collected from oceanographic profiling systems. 
The primary goal in the development of ProSoft was to 
create a package capable of processing optical profile data 
in an automated manner to remove the subjective aspects 
of data processing, thereby ensuring any two investigators 
would get the same derived values from the same data set. 
ProSoft allows users to read, process, visualize, and pro- 
duce numerous data products from raw data sets collected 
from their instrumentation. 

The ProSoft program was developed over the past ten 
years by Dalhousie University (by the authors) and dis- 
tributed to users as source codef. ProSoft currently con- 
sists of about 220 program modules and 60 support files 
linked by a GUI. The current version of ProSoft is 6.3d, 
and runs under MATLABf version 6. 

This chapter is intended to give an overview of the main 
data processing steps carried out by ProSoft. It is as- 
sumed that the radiometric data have been collected with 


f The ProSoft program is currently supported by Satlantic, 
Inc. and is freely available at the following Internet address: 
ftp . satlantic . com/pub/sensors/sof twar 2 /prosof t. 

t MATL AB is a registered trademark of MathWorks, Inc. (Nat- 
ick, Massachusetts). 


an optical instrument whose raw data format is in compli- 
ance with the Satlantic Instrument Files Standard (SIFS) 
format § or with GSFC level -0 ASCII files. It is also as- 
sumed that the radiometric data have been collected as 
a vertical profile of the water column with coincident ref- 
erence measurements from an instrument floating at the 
ocean surface. 

The ProSoft data processing steps are organized into 
different levels which follow, as much as possible, the pro- 
cedures given in the Ocean Optics Protocols for SeaWiFS 
Validation, Revision 1 (Mueller and Austin 1995). The 
processing blocks which exactly follow the Mueller and 
Austin (1995) recommendations are explicitly indicated. 
During processing, detailed notes, equations, and refer- 
ences for the various modules can be obtained by clicking 
on the Help button in any of the ProSoft (internal) menus. 

4.2 DATA PROCESSING 

ProSoft processing is segmented into four different lev- 
els: a) level -0 data are raw (usually binary) files from 
an instrument; b) level-1 data are calibrated (engineer- 
ing units), but unaveraged instrument data; c) level -2 data 


§ Satlantic Instrument Files Standard, Satlantic ’s Data For- 
mat Standard for Calibration and Telemetry Definition Files, 
2000, Version 6. 
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are averaged and corrected instrument data; and d) level -3 
data are the final products, i.e., R rs , [L\y (A)] n> PAR, etc. 
Processing can be done manually (and interactively) by 
the user, guided by the ProSoft processing steps, or in a 
fully automated batch mode. User preferences, settings, 
and instrument configuration parameters are set in a con- 
figuration file (prosof t . cfg). ProSoft will output data 
files in three formats; tab-delimited ASCII text in ProSoft 
format, SeaBASS text format, or MATLAB binary format. 
The latter can be reingested by ProSoft much faster than 
the ASCII formats. 

4.2.1 Level-1 Data Processing 

Level -1 data processing begins with reading the raw 
level -0 data produced by Sat View (the Satlantic data ac- 
quisition program) compliant with SIFS 2000 (or GSFC 
level -0 data). These files typically contain detailed header 
information on the cast, in addition to the binary telemetry 
from one or more instrument systems. An instrument sys- 
tem usually consists of a profiling optical system (profiler) 
and a surface reference at some fixed depth or location 
above the sea surface (reference). ProSoft scans the data 
files and determines the number and type of instruments, 
prompting the user for the appropriate calibration files (or 
automatically loading the most recent file for that system 
in the current directory if the program is being used in an 
automated mode). The calibration information for each 
channel is then applied to the raw data to convert it to 
engineering units. 

4. 2. 1.1 Optical Data Calibration 

Standard optical sensor data formats are processed dif- 
ferently, based on the capabilities of the various types of 
acquisition systems. These are referred to as 0PTIC1 (high 
resolution gain switching 24-bit systems), 0PTIC2 (stan- 
dard 12-, 16-, 24-, or 32-bit systems), and 0PTIC3 (hyper- 
spectral systems with adaptive integration). Application 
of the calibration data to all optical and ancillary sensors 
is carried out in accordance with the procedures detailed 
in SIFS. 

Optical data are converted into engineering units using; 

6(A ,U) = C c (\) I f (X) (V(A,i,) - D(A)], (20) 

where 6(A,U) is the calibrated value from the Satlantic 
processor, for which 0 is replaced by Ed(z) : E u (z ), L u {z ), 
or Ed{ 0 + ); C c (A) is the calibration coefficient; //(A) is the 
immersion factor; V(A,U) is the raw voltage (in digital 
counts) measured by the instrument at time t{ (which also 
sets the depth, 2); and D(A) is the dark value (in digi- 
tal counts). Within the ProSoft nomenclature, the dark 
value in (20) is referred to as the dark current , and can 
be obtained from a calibration file or from one of the dark 
current correction (DCC) methods described below. 


ProSoft currently implements DCC options other than 
calibration dark or shutter dark (which are explained be- 
low) only in the OPTIC 1 (high resolution 24-bit systems) 
fitting mode.f The dark current can change under varying 
thermal conditions on the high resolution systems and has 
to be adjusted accordingly. 

In the 0PTIC1 fitting mode, DCC other than that based 
on the calibration file (except HyperTSRB) can only be 
used for measurement frames obtained with a gain switch 
value of 1 or higher. In such a case, an analysis is first made 
of the measured frames according to the gain switches 
and applies the DCC according to the required correction 
scheme (e.g., PRO, BIN, or others). In the 0PTIC2 fitting 
mode, there are only two options available: CAL or NUL. 

4-2. 1.1.1 CAL Darks 

The CAL DCC method is given by (20) and is based 
on using the dark values recorded in the calibration files, 
D C (A), for the instrument systems: 

©(A, U) « C c (X)I f (X)[V(X : t l ) - D c ( A)]. (21) 

This is the default mode for both OPTIC 1 and 0PTIC2 data 
types. 

4 . 2. 1.1. 2 NUL Darks 

The NUL dark correction scheme is a special mode in 
which no darks are subtracted during data calibration (note 
that in the 0PTIC1 fitting mode, CAL darks are still sub- 
tracted for the data frames with a gain switch of zero): 

©(A, U) = C' c (A) //(A) V (A, ti). (22) 

The NUL DDC scheme for 0PTIC3 data is explained in a 
later section. 

4.2. 1.1.3 PRO Darks 

The PRO DCC option is a special mode for OPTIC 1 data 
where a separate dark cast is taken. Immediately before 
or after a regular cast (profiler measurements with no caps 
on the sensors), a dark cast (profiler measurements with 
caps on the sensors) can be measured. Both profiles are 
first processed with the NUL option (i.e., no darks are sub- 
tracted), and then a second dark-corrected light profile, 
0(2, A), is obtained as a difference between the regular 
light profile, £(2, A) and the binned dark profile, (©(2, A)): 

0(2, A) = £(2, A) - (33(2, A)), and 

1 JL (23) 

(®(*,A)> = mE®( 2 - a ) 

i— 1 

where the correspondence between the light profile and 
the binned dark profile occurs within plus or minus half a 
binning interval. Dark data binning (averaging) is carried 
out within 1 m depth intervals. 


| Explanations for the fitting modes are contained within SIFS. 
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4.2. 1 . 1.4 BIN Darks 

If the profiler attains a depth for which all the optical 
sensors reach their dark levels, the darks can be computed 
from the average of a number of samples at the bottom 
of the profile. For each wavelength, the value of DCC is 
obtained from a layer, defined by z m [ n and 2 max , where the 
average minimum light values, (£ m i n (z, A)), are observed: 

6 ( 2 , A) - £( 2 , A) - (£ mi n(*,A)>, (24) 

where z = 2 m j n , . . . , 2 max . 

4.2. 1.1.5 LUB Darks 

For each wavelength, the DCC factors for L u ( A) values 
are obtained from a layer wherein average minimum light 
values are observed (24). At the same time, Ed(\) values 
are dark corrected with a calibrated dark value, £> C (A): 

0 ( 2 , A) = £( 2 , A) — S c (^)i (25) 

where 2 — 2 m j n , . . . ,z max . 

4 . 2. 1.1 . 6 LOG Darks 

Prior to, or just after, a light profile (profiler measure- 
ments with no caps on), darks can be measured on the deck 
with the profiler caps in place. Both data files are first pro- 
cessed with the NUL option (i.e., no darks are subtracted), 
then the LOG dark-corrected light profile is obtained as a 
difference between the original light profile, £( 2 , A), and 
the dark values, which are obtained at each wavelength as 
the result of time averaging all the calibrated dark data, 
2 (A): 

6 ( 2 , A) = £( 2 , A) - ®(A). (26) 

Note that in this case, the PRESSURE TARE value is not 
subtracted at the time of data processing with NUL option. 
The value of PRESSURE TARE is estimated as an average 
pressure of the measurements carried out on deck. The 
PRESSURE TARE value is subtracted during LOG processing. 

4.2. 1.1.6 DCC of 0PTIC3 (Hyper spectral) Data 

Hyperspectral data are usually dark corrected with the 
values obtained from so-called shutter darks to obtain the 
most accurate correction. Shutter darks are continuously 
recorded during the measurements by occulting the input 
fiber with an optical shutter, typically after every five light 
samples. To calibrate hyperspectral data, the NUL DCC 
option should be used first, and the hyperspectral calibra- 
tion and subsequent DCC factors are determined in the 
following steps: 

1 . Correct the shutter dark counts obtained from a log 
file, D( A), by the dark offset, D 0 (A): 

D s ( A) = D( A) - D 0 ( A). (27) 

where D s ( A) is the corrected shutter dark values. 


2 . Convert the digital (voltage) dark counts into en- 
gineering units in accordance with the calibration 
equations. The calibration equation for hyperspec- 
tral data is: 

2(A) = C c (A)/ / (A)D 5 (A)L (28) 

1 2 

where C c (A) is the calibration coefficient, //(A) is 
the immersion coefficient, £2 is the integration time 
during calibration, and t\ is the integration time 
during the measurement. The former three are ob- 
tained from a calibration file, and the latter is ob- 
tained from the same log file as the optical data. 

3 . Check the sequence of frame numbers, and flag the 
frames (i.e., set to a value of -999) that are not in 
sequence. 

4. Deglitch the dark data using a first difference filter 
(optional). 

5. Smooth the shutter darks using a running boxcar 
filter. 

6 . Interpolate shutter darks as a function of measure- 
ment time to match the number of dark and light 
data measurements. 

7. Correct the light data using the calibrated and cor- 
rected shutter darks: 

6 ( 2 , A) - £( 2 , A) - D(A). (29) 

8 . Correct the light data using the thermal responsiv- 
ity correction: 

0 ( 2 , A) = 0 ( 2 , A) [l + 0.01(ciA 3 4- C 2 A 2 + 

1 ( oU ) 

c 3 A + c 4 )(T - 20)] , 

where ci, C 2 , c 3 , and C 4 are constants, and T is the 
temperature of the radiance or irradiance spectro- 
graph (d = 6.79131 x 10" 9 , c 2 = -1.09902 x 1CT 5 , 
c 3 = 6.51646 x 10" 3 , and c 4 = -1.31056). 

4. 2. 1.2 Range Check 

Data that is suspect, out of range, or invalid is au- 
tomatically flagged (i.e., set to a value of -999). Sensor 
values of -999 are not used for any computations and are 
ignored. On graphs, these values are not shown and, thus, 
are displayed as a missing section of data. The thresh- 
old limits for various sensors, which determines whether 
or not a particular value will be out of range, are set in the 
configuration file (prosof t . cf g). 

4. 2. 1.3 Frame Check 

ProSoft checks whether a frame counter was present. 
If it was, the sequence of frames, which should increase 
monotonically until the maximum value is reached where- 
upon the frame counter resets, is checked. Any frames out 
of sequence are flagged. 
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4. 2. 1.4 Rate Check 

If the reference and profiler measurements are com- 
bined into a single data file, the data rates of both instru- 
ments (i.e., the number of frames per second) are checked 
(where time stamps are not available in the data). If the 
rates of both instruments are not equal, the following pro- 
cessing steps are executed (time stamps are optional dur- 
ing data acquisition, so the rate check and frame sequence 
procedures do not rely on them): 

1. Information about the sequence of the reference and 
profiler frames is obtained from a raw binary file 
(the sequence can be obtained, because profiler and 
reference frames have unique frame counters): 

2. Blank frames (i.e., frames with the variables set to 
values of -999) are added to the frames with the 
lowest data rate. This ensures that the number of 
frames corresponds to the measurement rates of the 
instruments: and 

3. The blank frames that belong to the instrument 
with the lowest sampling frequency are removed. 

Note that some blank frames could remain if the sequence 
of the reference and profiler data frames was not in strict 
accordance with the respective sampling frequencies. 

4. 2. 1.5 Profiler Data 

Profiler data are subjected to a number of corrections 
and checks before the level -1 data files are created. These 
include a) pressure tare (offset) correction, b) pressure 
spike check, c) wavelength sequence check, and d) a tilt 
quality control edit. 

4. 2. 1.5.1 Pressure Tare Correction 

Pressure sensors have a small offset value that varies 
with atmospheric pressure. The offset is removed by a pres- 
sure tare, which is stored in the header records by SatView 
when the data are acquired. The pressure tare is defined 
as the pressure value that is measured when the Ed sensor 
on the instrument system is located just at the surface. 
This value has to be subtracted later from the measured 
pressure values of a profile. If the profiler measurements 
were logged with SatView', the pressure tare value is op- 
tionally inserted into the file header. If no information is 
obtained, the pressure tare value is set to zero. 

4-2. 1.5. 2 Wavelength Sequence Check 

The wavelength sequence check resequences the optical 
data so that the data columns are organized from left to 
right as a function of increasing wavelength. 

4-2. 1.5.3 Pressure Spike Check 

All frames with obvious pressure spikes are flagged. A 
deglitching procedure is applied using a threshold factor of 
100 (Sect. 4. 2. 1.6.2). 


4. 2. 1.5. 4 Tilt Edit 

Profiler measurements can be edited interactively or 
automatically to remove any data records with unaccept- 
able tilt values. In manual mode, a plot is shown with 
profiler tilts (if available) and profiler velocity as a func- 
tion of depth. With the interactive graphical editor, the 
user defines an upper (surface) z mm and lower z max depth 
(the z coordinate is directed downwards). Data above z m \ n 
or below z max axe removed. 

A tilt threshold, Tiltlimit (the default value is 
(y? max =5° in the prosoft.cfg file), is usually used for es- 
timating z m i n and z m ax- In the case of a missing (or in- 
operative) tilt sensor, the free-fall velocity of the profiler 
can be used (zero or negative velocity is assumed to be 
associated with high profiler tilts). An automated editor 
module determines z m in and z max in 10% of the upper and 
lower part of the profiler measurements. The upper depth 
and the lower depth are defined as the first depth level (z) 
where the following conditions are satisfied: 

1. ip(z) < V?max 5 

2. z m i n < 2 < 0.1 Zmax for the top layer of the profile: 
and 

3. (zmax — O.lZmax) < z < z max for the bottom layer 
of the profile. 

If a reference instrument is present, then the matching 
records in the reference data are also removed. 

4.2. 1.5.5 Output Files with the PRO Extension 

The results of processing level -0 profiler data are saved 
in tab-delimited MAT (an optional MATLAB internal for- 
mat) and ASCII formats. 

4. 2. 1.6 Reference Data 

Reference data are checked before the level -1 data files 
are created. These checks include the wavelength sequence 
check and spike deglitching. 

4. 2. 1.6.1 Wavelength Sequence Check 

The wavelength sequence check resequences the optical 
data so the columns are organized from left to right as a 
function of increasing wavelength. 

4. 2. 1.6. 2 Data Deglitching 

The (optional) deglitching module is applied to data to 
remove obvious glitches in sensor performance. For refer- 
ence data, the process begins with the application of a first 
difference operator to the time series of N data frames at 
each channel: 

^(A^i) “ ^d(b ) A £i+ 1 ) — Ed(0 , A .ti). (31) 

where t t sets the sample index (in time or by frame count 
number), and the calculations occur over the interval i — 
1,...,V- 1. 
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The standard deviation for each channel, a( A), is then 
calculated as: 


cr(A) - 


1 

TV -2 


N - 1 


0.5 




(32) 


i— 1 


where 5(A) is the average difference in the Ed{ 0 + ,A) time 
series (of TV frames). Any frames with differences more 
than 10<r from the average value are flagged. 

The general process described above is used with other 
deglitching requirements, except a new threshold value is 
used to match the type of deglitching under consideration 
(e.g., pressure spike removal uses a threshold of 100). 

4.2. 1.6.3 Output Files with the REF Extension 

The results of processing leveTO profiler data are saved 
in ASCII and in MAT formats. 

4.2.2 Level-2 Data Processing 

The level -2 ProSoft data processing produces averaged 
(binned) profiles from the calibrated level -1 data as level - 
2a output. The data from the profiler sensors are also 
adjusted to the same reference depth (to account for offsets 
in depth between the sensors). In level -2b, the reference 
correction of the profile for surface irradiance fluctuations 
can be applied. The smoothed and adjusted data are then 
ready for the production of the level -3 data products. 


is performed as a function of time. If a profiler, for some 
reason, was kept at a fixed depth level for a long time, how- 
ever, binning is still carried out as a function of time. A 
criterion for depth or time binning is based on the effective 
depth resolution. 

The effective depth resolution is the ratio of the number 
of data points measured between the maximum and mini- 
mum pressure indices to the maximum change in pressure. 
If the effective depth resolution is less than 1 mm, binning 
is carried out as a function of time; otherwise binning is 
carried out as a function of pressure. 


4.2.2. 1.4 Calculation of Time and Velocity 

In addition to the measured variables, level-2 data in- 
cludes the variables that are derived from the measured 
ones: time (t) and velocity (u). Note that time stamps, 
which can be added to a hie during data acquisition, are 
optional. For this reason, ProSoft cannot rely solely on 
time stamps, but rather has to generate the time for each 
data frame using the start time, to (in hours), and the 
sampling frequency, f s (in hertz): 


2 — 1 

U = to + 3600 f 3 ’ 


(33) 


where the calculations are performed over the interval 1 = 
l,...,iV. 

The velocity at each data frame is computed as 


v(ti) = /. [z(t l+1 ) - z(u)}, (34) 


4. 2. 2.1 Level-2a Processing 

The purpose of the level -2a processing is to provide or 
prepare for most of the binning and filtering options. 

4.2.2. 1.1 Calculation of Tilts 

From the level -1 data, the x and y components of tilt 
measurements are obtain ed in degrees. The tilts are calcu- 
lated using = tan -1 ^/tan 2 <p x + tan 2 c^ y , where <p x and 
{p y are the two- axis tilts with respect to the vertical. 

4 . 2. 2. 1.2 Constant Wavelength Intervals 

Interpolation to a constant wavelength interval can be 
carried out only with hyperspectral data. The interpo- 
lation interval is defined in the prosoft. cfg file by the 
DELTA _WL variable. If DELTA.WL = 0, no wavelength inter- 
polation is performed. 

4.2.2. 1.3 Depth and Temporal Binning 

Data processing from level- 1 to level -2a performs data 
averaging (referred to in ProSoft as “binning ) as a func- 
tion of pressure (depth) or time. If both profiler and ref- 
erence or only profiler data is available, then the binning 
is performed as a function of pressure. Otherwise, binning 


and calculated over the interval i = 1, . . . ,TV. After the 
computation is made, the velocity is smoothed using a box- 
car filter with a filter width of M. The value of M is the 
product of VEL_SM00TH_TIME times f s (VEL_SM00TH_TIME 
has units of seconds and is set within the prosoft. cfg 
file). 

4.2.2. 1.5 Binning 

Binning is performed in several steps. In the first step, 
each instrument (i.e., profiler or reference) form an index 
vector that includes the sequential indices of the frames. 
The frame indices, F t) depend on the measurement rate of 
an instrument in the following way: 

Fi(i) = l + (»-l)/r, * = ( 35 ) 

where N is the number of frames (measurements) carried 
out at a maximum rate, and f r is an integer part of the 
rates ratio given by MAXIMUM_RATE// S (j) for j — 1, . . . ,TV/ 
and Nj is the number of instruments (usually TV/ = 2). If 
the rates of all instruments are equal, Fi(i) — i over the 
interval i — 1, ... ,TV. 

In the second binning step, a coordinate vector must be 
chosen (i.e., time or pressure) corresponding to the binning 
that will be performed (Sect. 4. 2. 2. 1.3). If depth (pressure) 
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binning is selected, the data are interpolated to ensure the 
difference between sequential frames is a constant (the de- 
fault resampling is 0.1m). 

The third binning step adjusts the location between the 
Ed sensor and the other optical sensors based on the depth 
differences between the sensors. For example, the pressure 
measured at the L u sensor is the pressure measured at 
the Ed sensor minus the PRO_LU_DIST value given in the 
prosof t . cf g file. 

The binning along the chosen coordinate vector is car- 
ried out in the fourth step. Currently, it is assumed that 
profiler measurements are always carried out at a higher 
rate, and prior to binning, the optical data are log trans- 
formed (whether or not the data are log transformed de- 
pends on switches given in the prosof t . cfg file). 

In the fifth and final step, the profiler optical data, 
other than Ed, are adjusted to match the depth of the first 
Ed measurement (in the prosof t. cfg file, the respective 
sensor distance offsets are given by PR0_SENS0RNAME_DIST, 
where SENSORNAME is, for example, Lu or Eu). Extrapola- 
tion is based on a linear regression fit of the log-transformed 
data versus pressure. The number of points used in the lin- 
ear regression calculation is defined in the prosof t . cfg file 
by the variable NUM_K_BINS. 

4. 2. 2. 2 Level-2b Processing 

The purpose of the level -2b data processing is to com- 
ply with the SeaWiFS Ocean Optics Protocols (Mueller 
and Austin 1995) and as formulated in Smith and Baker 
(1984 and 1986), by weighting the in- water profiler data by 
the surface irradiance. The weighting is carried out using: 

= S( t ,A,t) Ed{0 . jX t) • (36) 

but note that this procedure is different from the alternate 
method also given in Mueller and Austin (1995) using the 
solar irradiance measured when the profiler was at the sur- 
face (z = 0”)f. The reason for this difference is mainly 
related to the problem- that £^(0", A) is often unavailable, 
so ProSoft uses (36) instead. 

Level -2a to level -2b processing automatically follows 
the level -1 to level -2a processing if the REF_C0RRECT_PR0 
variable is enabled (turned ON) in the prosoft. cfg file. 
Note also that the weighting given by (36) often introduces 
significant variations in the profiler data, particularly when 
clouds are present, so the default is for REF_C0RRECT_PR0 
to be disabled (set OFF). 

4.2.3 Level-3 Data Processing 

The level -3 data processing is based mostly on the 
level -2 data. An exception is the estimation of the dif- 
fuse attenuation coefficient and the related level -3 prod- 
ucts from the level -1 (unaveraged) radiometric data. In 


t See Mueller and Austin (1995), Sect. 6.1.3, page 49. 


general, level -3 data processing can be divided into two 
parts: a) calculation of the products at the surface, and 
b) calculation of the products based only on the profiler 
measurements. 

The calculation of the products at the surface involves 
three steps: 

1. Propagation of all the radiometric data to the null 
depth, 2 = 0~, level; 

2. Propagation of the radiometric data through the 
air-sea interface and, thus, above the surface, 2 -- 
0 + ; and 

3. Calculation of the various products, such as, the 
remote sensing reflectance, reflectance, normalized 
water-leaving radiance, diffuse attenuation coeffi- 
cient, PAR, and pigment concentration. 

In addition to the surface products, level -3 (vertical) 
profile products are calculated and involve the estimation 
of the diffuse attenuation coefficient, PAR, pigment con- 
centration, remote sensing reflectance, and reflectance. 

4.2.3. 1 Surface K 

Surface values of K( A) and the primary radiometric 
variables, i.e., E d ( A) and L u ( A), are obtained from a least- 
squares regression fit to the first N (looking down from 
the surface) log-transformed radiometric and depth (or 
pressure) data points, expressed here in terms of the data 
indicesf: 

ln[s(zj,A)] = In [0(29, A)] — K< 5 (zo,\)zi, (37) 

where 2* sets the sample index in terms of the depth (or 
the frame count), zq denotes the surface value (or initial 
frame count), and the calculations occur over the interval 
i = 1 , . . . , N — 1 . 

4. 2. 3. 2 K Profile 

A vertical profile of K( A) is obtained from a least- 
squares regression fit which is applied to the N data points 
at a time 

In [ &(z l , A)] = In [ 0(2 m , A)] - 

A)(2* — z m ) , 

where i = m + 1 , . . . ,m -f N and m = 1 , . . . ,M - N , and 
M is the number of data points at a profile. 

4. 2. 3. 3 K from Fixed Depth Measurements 

In some cases, K( A) has to be estimated from a time 
series of upwelling radiance measurements, collected at a 


J Formulated from (31) on page 50 of Mueller and Austin 
(1995), which is based on the Smith and Baker (1984 and 
1986) traditional method of K analysis. 


62 


S. Hooker, G. Zibordi, J-F. Berthon, D. D’Alimonte, S. Maritorena, S. McLean, and J. Sildam 


fixed depth (approximately 1 m) below the surface (this 
type of measurement is a common feature of nonprofiling 
instruments, like surface buoys). In this case, ProSoft esti- 
mates the spectrum of K( A) based on an empirical model 
(valid for Case-1 conditions): 

1. Compute K( 490) and K( 520) from L u ( 443) and 
L u (550) using the Austin and Petzold (1981) al- 
gorithms: 


AT(490) = 0.022 + 0.0883 
and 

K( 520) - 0.044 + 0.0663 


L u (443) 
L u { 550) 


- 1.491 


M 443) 

M550) 


- 1.398 


(39) 


2. The Morel (1988) Case-1 empirical model for esti- 
mating K is inverted to obtain the pigment concen- 
tration C (chlorophyll a plus phaeopigments): 


C490 = 

C520 — 
and 

C = 


£(490) — AT (490) 
Xc( 490) 

£(520) - K w { 520) 
Xc(520) 

^490 + C*520 


,e" 1 (490) 


1 e~ 1 (520) 


(40) 


as, Lw — L u { 0“,A)(l - p(A, $))n” 2 (A), where p( A,tf) is 
fixed at 0.021 for normal incidence, and n^(A) = 1.345 for 
seawater (Austin 1974). For the latter, if an above- water 
reference is unavailable, the incident solar irradiance is de- 
rived from the profiler Ed(z,\) measurements: 

E d ( 0 + ,A) = (1 -f a) E d (0~, A), (42) 

where a is the Fresnel reflection albedo for irradiance from 
sun and sky (Gordon et al. 1988) which is fixed at 0.043. 
For these formulations, p and a are assumed independent 
of wind speed (Priesendorfer and Mobley 1986). 

4. 2. 3. 5 Normalized Water-Leaving Radiance 

The normalized water-leaving radiance is estimated us- 
ing (63) of the SeaWiFS Ocean Optics Protocols (Mueller 
and Austin 1995): 

[Lh/(A)] n = L w J^- y , (43) 

where, in ProSoft, the Neckel and Labs (1984) mean ex- 
traterrestrial solar irradiance data are reinterpolated to a 
0.1 nm resolution and then smoothed to match the 10 nm 
bandwidth of the radiometer detectors. The Lw{ A) term 
in (43) is computed following the discussion in Sect. 4. 2. 3. 4. 

4. 2. 3. 6 PAR 

PAR, in units of quanta cm -2 s _1 , is estimated using: 


where AT (A) represents the diffuse attenuation co- 
efficient for pure water, and the coefficients Xc(^) 
and e(A) are regression coefficients determined by 
statistical analysis of £( A) in Case-1 water for the 
wavelength range 400-700 nm (Morel 1988). 

3. The algorithm for computing £( A) is as follows: 

I<{ A) = K w ( A) + X c(A )C e ( A ). (41) 

To calculate £( A) within the 700-800 nm range, Prosoft 
uses the following steps: 

i) The values of AT (A) within 700-800 nm are ob- 
tained from Smith and Baker (1981); 

ii) It is assumed that Xc(^) and e(A) decrease linearly 
from the values given at Xc(?00) and e(700) to 0 at 
800 nm; and 

iii) Using the above assumptions and AT (A) within 
700-800 nm, K(X) can be estimated for the 700- 
800 nm interval. 

4. 2. 3. 4 Remote Sensing Reflectance 

The remote sensing reflectance is calculated using the 
ratio of the water-leaving radiance and incident solar ir- 
radiance, A rs (A) = Lw(X)/ Ed(0 + , A). Lw( A) is defined 


7 700 A 

PAR - / — E d {\) d\, (44) 

J 400 

where h is Plank’s constant, and c is the speed of light. 
The PAR formulation (44) is used for estimating PAR at 
the surface as well as for vertical profiles of PAR. ProSoft 
uses the MATLAB trapezoidal numerical integration rou- 
tine trapz.m to perform the integration in (44). Prior to 
integration, E d ( A) is interpolated onto a constant 1 nm in- 
terval and a minimum-maximum wavelength check is per- 
formed: 

if Ai > 400 => £**(400) = E d { A x ), and 
if A/v < 700 => £^(700) = E d { A^v). 

4. 2. 3. 7 Energy Flux 

The energy flux, £#, is computed as: 

r700 

F e - 100 / E d (\)d\, (46) 

^400 

ProSoft uses the MATLAB trapezoidal numerical integra- 
tion routine trapz.m to perform the integration in (46). 
Prior to integration, E d ( A) is interpolated onto a constant 
1 nm interval, and a minimum-maximum wavelength check 
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is performed as in (45). Fe is given in units of W m~ 2 , and 
because E d is in units of \i Wcm" 2 nm _1 , a scaling factor 
of 100 is used in (46). The energy flux equation (46) is 
used for estimating Fe at the surface, as well as for verti- 
cal profiles of Fe- 

4. 2. 3. 8 Estimating K Using Wavelet Analysis 

This section describes a new methodology for comput- 
ing K using the K from multiresolution (KMR) wavelet 
analysis method (Sildam et al. 1998). This level -3 module 
uses high resolution level -1 (unaveraged) data instead of 
the standard (Smith and Baker 1984 and 1986) method 
which uses level -2 (averaged) data. 

The following modules in ProSoft are experimental and 
are intended for advanced users. These subroutines do not 
follow the SeaWiFS Ocean Optics Protocols (Mueller and 
Austin 1995) and should not be used for submitting data 
to SeaBASS until further validation has been done. 

Estimation of K from optical profiles involves ambi- 
guity in the choice of the vertical scale. Usually, a noisy 
profile of K is obtained if K is estimated directly from: 


4 .2. 3. 8.1 Mathematics of the KMR Method 

Multiresolution orthogonal Haar wavelet transforms de- 
compose the K(z , A) profile into M profiles at M scales, so 
that K(z m ,X) at different scales are related to each other 
as: 

K(z m . i,A) - K(z mt A) + K'(z m , A), (48) 

where m — 1, ... ,M and K'(z m . A) is the so-called detail. 
K(z m , A) is given by: 

oc 

K(Zm,\) — ^ ^ n 0 7 n,n(z') . ( 49, 1 

71= — OC 

where 

U*m,n — A"(^ ) < Pm.n (z^dz . 

1 J /z-n2 m 5z\ ] / N 

0m:n = — 2 ^ — )' and (50 ' 

f 1, 0 < 2 < 1; 

<j>(z) — s 

(0, otherwise. 


K(z, A) 


-1 A E d 

E d (z, A) Az 


(47) 


where A E d = Ed(z n +i . A) — E d (z n ,\) is calculated at 
the resolution of the measurements, Az. If A E d is esti- 
mated at an increased vertical scale mAz, i.e., A Ed — 
Ed{zn+m, A) - E d (z n . A) , smoother — but less detailed — 
profiles of K are obtained. Note that a similar problem 
holds also for estimating K using (38). The objective of 
the KMR methodology is to define an optimal depth in- 
terval for the determination of K that is objectively and 
automatically determined. 

Estimation of the K profile from discrete irradiance or 
radiance measurements makes the assumption that any ob- 
served profile can be divided into optically homogeneous 
layers. Usually the choice of the thickness of these layers 
is subjective. Moreover, it varies with wavelength, depth, 
time, and the geographical location of the measurements. 
Multiresolution wavelet decomposition of a K(z, A) profile 
(that is calculated with a fixed high vertical resolution) 
provides a hierarchy of profiles of K(z m: A) which are ob- 
tained by smoothing K(z, A) at different vertical scales. 
An optimal choice of the averaging scale at any depth level 
could then be based on a statistical approach (i.e., on the 
analysis of the K frequency distribution). 

In practice, having a few K estimates at a fixed depth, 
results in an inconclusive analysis of the K frequency dis- 
tribution. Frequency analyses, however, can be carried 
out over all observed depth levels if K is centralized and 
normalized by its standard deviation. The results of such 
an analysis show that in optically clear waters, recursive 
wavelet decomposition of the K profile combined with me- 
dian averaging can be used to obtain K profiles, which 
asymptotically reach the state of maximum probability 
(Sildam et al. 1998). 


In (49), the depth location and scale are given by the n and 
m indices, respectively, <p m ,n is the so-called smoothing 
function , and </>(z) is Haar’s scaling function. 

In the KMR method, at every depth level z, m — 
1,...,A/ estimates of K(z m , A) can be obtained. From 
these M profiles, a median profile (K(z. A)) can be eas- 
ily obtained. The median profile can be decomposed once 
again using the Haar scaling function, and a new median 
profile can be obtained. This iterative procedure can be 
repeated several times until an asymptotic K( A) is reached. 


4- 2. 3. 8. 2 Practical Steps of the KMR Method 

The practical steps of the KMR method are as follows 

1. Estimate the K(z, A) profile using 


K(z, A) 


1 ln E d (zj, A) 
Az _ Ed(zi + n , A) _ 


(51) 


where Az = z l+n - z 7 represents the maximum al- 
lowed resolution (if i — 1, the maximum allowed 
resolution is equal to the depth resolution of the 
measurements). 

2. Apply the Haar multiresolution wavelet decompo- 
sition (50) to K(z, A). As a result, obtain M pro- 
files of AT(z m , A), each calculated with a depth scale 
given by (zi - z l+ i)2 m , where m = 1, . . . ,M. Note 
that n and m represent the minimum and maximum 
number of points, respectively, used in estimation of 
K(Zm, A). 

3. From the K(z rn , A) profiles (m — 1 ,...,M), esti- 
mate a median profile, (K (z, A)). 

4. Repeat steps 2 and 3, each time using newly ob- 
tained (K(z, A)) until an asymptotic limit is reached. 
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4.3 SUMMARY 

The standard (40-cast) DARR-00 data were processed 
with ProSoft. Level -2 and level -3 products were submitted 
including [L w (A)] n . For all the data, the K analyses (us- 
ing the methods described in Sect. 4. 2. 3.1) were executed 
using the auto-edited profiles as produced by ProSoft. For 
Case-2 water comparison, K was computed with and with- 
out editing, ignoring near surface tilts for one example. 
All files (light and dark) were processed first using the NUL 
dark correction. Once all the data were converted to engi- 
neering units, the LOG dark correction was used to subtract 
the DARK casts from the LIGHT. Table 23 presents a sum- 


mary of the ProSoft settings in the prosoft.cfg file for 
each of the instrument systems. 


Table 23. The ProSoft configuration parameters 
used for processing the standard DARR-00 (40-cast) 
data set. 


Parameter 

MN 

WP 

LN 

so 

PROJLILDIST 

0.160 

0.055 

1.779 

0.183 

PRO_EU_DIST 

1.021 

0.045 

1.787 

0.181 

PRES.TAR2 

0.258 

0.030 

0.569 

0.177 

PRESJBIN.INT 

1 

1 

1 

1 

NUM_K_BINS 

5 

5 

5 

5 
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Appendices 

A. The DARR-00 Science Team 

B. The DARR-00 Database 

Appendix A 

The DARR-00 Science Team 

The DARR-00 science team is presented alphabetically. 

Jean-Frangois Berthon 
JRC/SAI/ME T.P. 272 
1-21020 Ispra (VA) 

ITALY 

Voice: 39-0-332-789-934 
Fax: 39-0-332-789-034 

Net: j ean-f rancois . berthon® jrc . it 

Davide D'Alimonte 
JRC/SAI/ME T.P. 272 
1-21020 Ispra (VA) 

ITALY 

Voice: 39-0-332-785-727 
Fax: 39-0-332-789-034 

Net: davide . cP alimonte® jrc . it 

Stanford Hooker 
NASA/GSFC/Code 970.2 
Greenbelt, MD 20771 
Voice: 301-286-9503 
Fax: 301-286-1775 

Net: stan®ardbeg . gsf c .nasa. gov 

Stephane Maritorena 

ICESS/UCSB 

Santa Barbara, CA 93106 

Voice: 805-893-2544 

Fax: 805-893-2578 

Net: stephane®icess . ucsb . edu 

Scott McLean 
Satlantic, Inc. 

Richmond Terminal, Pier 9 
3481 North Marginal Road 
Halifax, Nova Scotia B3K 5X8 
CANADA 

Voice: 01-902-492-4780 
Fax: 01 902-492-4781 

Net: scott®satlantic . com 

Juri Sildam 
Satlantic, Inc. 

Richmond Terminal, Pier 9 
3481 North Marginal Road 
Halifax, Nova Scotia B3K 5X8 
CANADA 

Voice: 01-902-492-4780 
Fax: 01-902-492-4781 

Net: juri®satlantic . com 

Giuseppe Zibordi 
JRC/SAI/ME T.P. 272 
1-21020 Ispra (VA) 

ITALY 

Voice: 39-0-332-785-902 
Fax: 39-0-332-789-034 

Net: giuseppe . zibordi®jrc . it 


Appendix B 
The DARR-00 Database 

A short summary of the data files used to create the DARR-00 
database is given in Table Bl. 

Glossary 

A/D Analog-to-Digital 
AAOT Acqua Alta Oceanographic Tower 
AMT Atlantic Meridional Transect 
AOP Apparent Optical Property 
ASCII- American Standard Code for Information Inter- 
change 

CoASTS Coastal Atmosphere and Sea Time-Series 
COLORS Coastal Region Long-Term Measurements for Col- 
our Remote Sensing Development and Validation 
CT Conductivity and Temperature (probe) 

DAAC Distributed Active Archive Center 
DARR Data Analysis Round Robin 
DARR-94 The First DARR (July 1994) 

DARR-00 The Second DARR (March 2000) 

DCC Dark Current Correction 
DCM Depth of the Chlorophyll Maximum 
DO Deep Ocean 

EOF End-of-File 
ET Eutrophic 

FF Free-Fall 

GMT Greenwich Mean Time 

GS GSFC and Satlantic (comparison) 

GSFC Goddard Space Flight Center 
GUI Graphical User Interface 

HPLC High Performance Liquid Chromatography 

ICESS (UCSB) Institute for Computational Earth System 
Science 

IDL Interactive Data Language 
IOP Inherent Optical Property 

JCR (Royal Research Ship) James Clark Ross 
JG JRC and GSFC (comparison) 

JGOFS Joint Global Ocean Flux Study 
JRC Joint Research Centre 
JS JRC and Satlantic (comparison) 

KMR K from Multiresolution (wavelet analysis) 

LN LoCNESS 

LoCNESS Low-Cost NASA Environmental Sampling System 

miniNESS miniature NASA Environmental Sampling System 
MLD Mixed Layer Depth 
MN miniNESS 
MT Mesotrophic 

NASA National Aeronautics and Space Administration 
NIST National Institute of Standards and Technology 

OC Ocean Color 

OC2v2 Ocean Color (algorithm) 2 Version 2 
OCI-200 Ocean Color Irradiance series 200 (sensor) 
OCR-200 Ocean Color Radiance series 200 (sensor) 

OLL One-Percent Light Level 
OT Oligotrophic 
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Table Bl. The DARR-00 data set. The data that was part of the 40-cast data set is indicated by “(40)” next to the 
file name; the data that was added to create the 50-cast data set is indicated by “(50)” next to the file name. The 
concentration of total chlorophyll a, Cr a , is given in units of milligrams per cubic meter. 


LoCNESS (FF, DO) 

SeaOPS (WC, DO) 

File Name 

Longitude [°] 

Latitide [°] 

Gr a 

File Name 

Longitude [°] 

Latitude [°] 

Cr a 

A07KD005 

(40) 

-10.7960 

41.5791 

0.28 

A060U025 

(40) 

12.4044 

-21.6534 

1.66 

A07KD020 

(40) 

-17.5773 

14.9881 

0.49 

A060D028 

(40) 

12.4042 

-21.6520 

1.66 

A07KD025 

(40) 

-21.3139 

12.1723 

1.18 

A060D029 

(40) 

12.4042 

-21.6529 

1.66 

A07KD028 

(40) 

-23.7503 

4.2843 

0.16 

A060D030 

(40) 

12.1015 

-21.3989 

1.79 

A07KD030 

(40) 

-23.7499 

4.2864 

0.16 

A060U031 

(40) 

12.1011 

-21.3984 

1.79 

A07KD040 

(40) 

-30.4334 

-12.1615 

0.08 

A060D032 

(40) 

11.9995 

-18.9960 

1.86 

A07KD045 

(40) 

-30.5133 

-12.3743 

0.08 

A060D033 

(40) 

11.3314 

-17.6599 

2.43 

A07KD060 

(40) 

-36.8965 

-23.4888 

0.13 

A060D034 

(40) 

11.3315 

-17.6601 

2.43 

A07KD070 

(40) 

-40.4902 

-27.0146 

0.29 

A060D035 

(40) 

11.0733 

-17.4396 

0.82 

A07KD090 

(40) 

-46.6759 

-32.7510 

0.10 

A060D039 

(40) 

-2.6308 

-5.8635 

0.24 

A06LD039 

(50) 

-19.9964 

25.0854 

0.03 

A070U005 

(50) 

-13.6048 

38.7815 

0.09 

A06LD041 

(50) 

-19.5624 

29.0873 

0.03 

A070D007 

(50) 

-17.4900 

36.3089 

0.06 

A06LD045 

(50) 

-17.1673 

32.6517 

0.03 

A070U010 

(50) 

-19.9454 

30.0624 

0.09 

A07KD009 

(50) 

-17.4967 

36.0950 

0.04 






A07KD047 

(50) 

-32.1382 

-16.2339 

0.04 






A07KD056 

(50) 

-35.1391 

-19.8250 

0.06 






A07KD092 

(50) 

-46.7654 

-32.8267 

0.07 






miniNESS (FF, SC) 

WiSPER (WC, SC) 

P56S1H 

(40) 

12.5083 

45.3139 

0.97 

W56S1A 

(40) 

12.5083 

45.3139 

0.97 

P56S4D 

(40) 

12.5083 

45.3139 

1.54 

W56S4A 

(40) 

12.5083 

45.3139 

1.54 

P56S7G 

(40) 

12.5083 

45.3139 

1.22 

W56S7A 

(40) 

12.5083 

45.3139 

1.22 

P58S4B 

(40) 

12.5083 

45.3139 

0.48 

W58S4A 

(40) 

12.5083 

45.3139 

0.48 

P60S4R 

(40) 

12.5083 

45.3139 

0.83 

W60S4A 

(40) 

12.5083 

45.3139 

0.83 

P60S6C 

(40) 

12.5083 

45.3139 

0.65 

W60S6A 

(40) 

12.5083 

45.3139 

0.65 

P61S5J 

(40) 

12.5083 

45.3139 

0.68 

W61S5A 

(40) 

12.5083 

45.3139 

0.68 

P61S6E 

(40) 

12.5083 

45.3139 

0.87 

W61S6A 

(40) 

12.5083 

45.3139 

0.87 

P62S4F 

(40) 

12.5083 

45.3139 

2.02 

W62S4A 

(40) 

12.5083 

45.3139 

2.02 

P63S3B 

(40) 

12.5083 

45.3139 

0.54 

W63S3A 

(40) 

12.5083 

45.3139 

0.54 


PAR Photosynthetically Available Radiation 

S/N Serial Number 
SC Shallow Coastal 
SDY Sequential Day of the Year 
SeaBASS SeaWiFS Bio-Optical Archive and Storage System 
SeaOPS SeaWiFS Optical Profiling System 
SeaWiFS Sea-viewing Wide Field-of-view Sensor 
SIFS Satlantic Instrument Files Standard 
SIRREX SeaWiFS Intercalibration Round- Robin Experiment 
SO SeaOPS 

SOOP SeaWiFS Ocean Optics Protocols 
SQM SeaWiFS Quality Monitor 

THOR Three-Headed Optical Recorder 

UCSB University of California at Santa Barbara 
UPD Unbiased Percent Difference 

VI Version 1 
V2 Version 2 
V3 Version 3 
V4 Version 4 
V5 Version 5 

WC Winch and Crane 


WiSPER 

WP 


a 

a" 


a 

A 

B 

a 

A 


C 

C 

C490 
C520 
C c ( A) 

Cm,n 

Cr a 

Cl 

C2 

C 3 

C 4 


Wire-Stabilized Profiling Environmental Radiome- 
ter 

WiSPER 


Symbols 

The spectral seawater absorption. 

The processor code for the first data product. 

The processor code for the second data product. 
The band-ratio data product being compared. 

The first band-ratio data product being compared. 
The second band-ratio data product being com- 
pared. 

The speed of light. 

Pigment concentration (chlorophyll a plus phaeopig- 
ments) . 

The pigment concentration contribution at 490 nm. 
The pigment concentration contribution at 520 nm. 
The calibration coefficient. 

A summation coefficient in the KMR method. 

The total chlorophyll a concentration. 

A thermal responsivity constant. 

A thermal responsivity constant. 

A thermal responsivity constant. 

A thermal responsivity constant. 


67 




Results of the Second SeaWiFS Data Analysis Round Robin, March 2000 (DARR-00) 


d The SDY. 

33 (A) The calibrated spectral dark value (in engineering 
units). 

S(zi) A dark (profile) value at depth Zi . 

A)} A binned spectral dark profile. 

D( A) The spectral (time) averaged dark value in digital 
counts (usually measured during a special dark cast 
with the caps on the radiometer). 

33(A) The spectral (time) averaged dark value in engineer- 
ing units. 

D c { A) The spectral dark values in a calibration file (in dig- 
ital counts). 

33 C (A) The calibrated spectral dark values from a calibra- 
tion file (in engineering units). 

D 0 ( A) The spectral dark offset (capped darks). 

D S (A) The corrected spectral shutter darks from a log file. 

dz A depth increment interval (z - zq). 

e(A) A spectral regression coefficient. 

F(f(0~,A) The spectral downward irradiance immediately be- 
low the sea surface. 

Ed{ 0 + , A) The incident solar irradiance immediately above the 
sea surface (direct plus diffuse). 

£^(0 + , A) The spectral mean of a set of Ed( 0 + , A) values. 

Ed(z . A) The spectral downward irradiance profile. 

Ei{ 0 + . A) The diffuse solar irradiance immediately above the 
sea surface. 

E u { 0", A The spectral upward irradiance immediately below 
the sea surface. 

E u (z,\) The spectral upward irradiance profile. 

5 The data product being compared. 

The first data product being compared. 

The second data product being compared. 

f r The integer part of the rates ratio. 

f s The sampling frequency. 

Fo(A, d) The mean extraterrestrial solar irradiance corrected 
for the Earth-sun distance. 

Fo(A) The mean extraterrestrial solar irradiance. 

Fe The energy flux. 

Ft The frame index. 

G A data product computed with the GSFC processor. 

0(A) The calibrated value from the GSFC processor [0 
is replaced by Ed(z ), E u (z), L u (z), or £Y(0 + )]. 

h Plank’s constant. 

i A sequential index. 

//(A) The (spectral) immersion factor. 

j A sequential index. 

J A data product computed with the JRC processor. 

3(A) The calibrated value from the JRC processor [3 is 
replaced by E d (z), E u (z ), L u (z ), or Ed{ 0 + )]. 

k A wavelength index. 

K( A) The spectral diffuse attenuation coefficient. 

(AT(z, A)) A median K( A) profile. 

K'(z m , A) The so-called detail in the KMR method. 

Kd( A) The spectral diffuse attenuation coefficient com- 
puted from the downward irradiance. 

Kd( A) The mean diffuse attenuation coefficient. 

Ks d ( A) The spectral diffuse attenuation coefficient compu- 
ted from the downward irradiance; also denoted 
K d { A). 


Ke u ( A) The spectral diffuse attenuation coefficient compu- 
ted from the upward irradiance; also denoted K u ( A). 

K 0 (A) The spectral diffuse attenuation coefficient compu- 
ted with the GSFC processor [0 is L u , Ed , or E u ]. 

A) The spectral diffuse attenuation coefficient compu- 
ted with the JRC processor [3 is L u , Ed . or E u ]. 

Kl( A) The spectral diffuse attenuation coefficient compu- 
ted from the upwelled radiance. 

Kl u ( A) The spectral diffuse attenuation coefficient compu- 
ted from the upwelled radiance; also denoted AY (A). 

AT 0 (A) The spectral diffuse attenuation coefficient compu- 
ted with the Satlantic processor [0 is L u or Ed]. 

AT (A) The spectral diffuse attenuation coefficient compu- 
ted from the upward irradiance. 

AT, (A) The spectral diffuse attenuation coefficient for pure 
water. 

I A wavelength index. 

£(z, A) A spectral light profile. 

(£m>n(A)} The average minimum spectral light values. 

L u (0”,A) The spectral upwelled radiance immediately below 
the sea surface. 

L u (z, A) The spectral upwelled radiance profile. 

Lw{ A) The spectral water-leaving radiance. 

[Lw(A)] n The spectral normalized water-leaving radiance. 

m A sequential index. 

M The boxcar filter width. 

n The number of standard deviations used in a statis- 
tical filter. 

ra^A) The spectral refractive index of seawater. 

N The number of data frames. 

Nb The number of spectral bands being averaged. 

Ng The average number of points in the extrapolation 
ranges for the GSFC processing. 

Nj The number of instruments (usually Nj =2). 

Nj The average number of points in the extrapolation 
ranges for the JRC processing. 

Nm The number of measurements being compared. 

Q n (0 - , A) The spectral nadir Q-factor immediately below the 
sea surface. 

r(A) The spectral diffuse-to-direct solar irradiance ratio. 

R( A) The spectral irradiance reflectance. 

Ars(A) The spectral remote sensing reflectance. 

S A data product computed with the Satlantic pro- 
cessor. 

0(A) The calibrated value from the Satlantic processor 
[0 is replaced by Ed{z), E u (z), L u (z ), or £Y(0 + )]. 

t The time coordinate. 

T The temperature of the radiance or irradiance spec- 
trograph. 

t x The time associated with a data sample. 

to A reference time (generally chosen to coincide with 
the start of a cast). 

t\ The integration time during a measurement. 

t “2 The integration time during calibration. 

v The velocity of a profiler. 

F(A) The raw (spectral) voltage (in digital counts) mea- 
sured by an instrument. 

x The abscissa coordinate. 
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y The ordinate. 

2 The vertical (depth) coordinate. 

2 o A reference depth. 

z\ The bottom (depth) of the extrapolation interval. 

zt, A depth near the bottom of a cast. 

z x The sample index in terms of the depth. 

Zmax The maximum depth for a vertical layer. 

Zmin The minimum depth for a vertical layer. 

a The Fresnel reflection albedo for irradiance from sun 
and sky. 

St A time interval. 

Sz The vertical bin size (in meters). 

<5(A,£i) The spectral first difference operator for times ti 
and U+ 1 . 

6 (A) The spectral average difference. 

6o(d) The Earth-sun distance correction. 

A Ed A difference in Ed between two depth levels (or in- 
dices). 

AN The average difference in the number of points in 
the extrapolation ranges for the JRC and GSFC 
processors. 

Az The maximum allowed resolution. 


6 The solar zenith angle. 

0 The nadir angle. 

A Wavelength. 

A 1-5 The blue-green wavelengths, nominally 412-555 nm 
(channels 1-5). 

A 3/5 The blue-green wavelength ratio, nominally 490 and 
555 nm (channels 3 and 5, respectively). 

\q -7 The red wavelengths, nominally 665-683 nm (chan- 
nels 6-7). 

A 1-7 All the wavelengths, nominally 412-683 nm (chan- 
nels 1-7). 

A k /i Denotes the wavelengths (by channel number) in- 
volved in a band-ratio data product for wavelengths 
k and l (the data product for wavelength k is divided 
by the data product for wavelength /). 

p( A) The spectral Fresnel surface reflectance. 
p s ( A) The sea surface reflectance for downward irradiance. 

o' The standard deviation. 

ip The sensor tilt with respect to the vertical axis. 

The maximum allowed tilt. 
p x The x-axis tilt. 
ip y The y-axis tilt. 

4 >^ B (A k/i) The UPD between two band-ratio data products 
B^Aj k/l ) and B b (A fc/( ). 

4 > 3/5 The average of the band-ratio results calculated for 
490 and 555 nm (bands three and five, respectively). 

<f>(z) Haar’s scaling function. 

0 mn The so-called smoothing function (used with the 
KMR method). 

Xc(A) A spectral regression coefficient. 

Vg B (A) The average UPD between two data products 
and $ B . 

4/^ B (A) The spectral average of the UPD between two data 
products $ A and 


4 >i -7 The spectral UPD average over all (seven) wave- 
lengths, 412-683 nm (channels 1-7). 

4 / 1 - 5 The spectral UPD average of the (five) wavelengths 
in the blue-green part of the spectrum, 412-555 nm 
(channels 1-5). 

4 / 6-7 The spectral UPD average of the (two) wavelengths 
in the red part of the spectrum, 665-683 nm (chan- 
nels 6-7). 

u;(A) The spectral single scattering albedo for seawater. 
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